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ABSTRACT 

A potentially powerful approach to quantum gravity has been developed over the last few years under 
the name of Causal Dynamical Triangulations. Numerical simulations have given very interesting 
results in the cases of two, three and four spacetime dimension. The aim of this thesis is to give an 
introduction to the subject (Chapter 1), and try to push the analytical understanding of these models 
further. This is done by first studying (Chapter 2) the case of a (l+l)-dimensional spacetime coupled 
to matter, in the form of an Ising model, by means of high- and low-temperature expansions. And 
after (Chapter 3) by studying a specific model in (2+1) dimensions, whose solution and continuum 
limit are presented. 
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Preface 



Since ancient times man has always been fascinated by the quest for the ultimate nature of things; the 
need for a deeper understanding has brought philosophers and scientists to look where eyes cannot 
see, be it the very distant or the very small. 

Meditating on the very small, Democritus, already twenty-four centuries ago formulated an atomic 
hypothesis. Many centuries later, the study of atomic phenomena had a fundamental role in the birth 
of quantum mechanics, which today is the basis of all microscopic physics. In particular, the Standard 
Model, which gives a complete description of the strong, the weak and the electromagnetic interactions, 
is built upon the rules of the quantum theory. 

On the other hand it was by looking at large distances that Newton was able to formulate his clas- 
sical theory of the gravitational force, the fourth and the last of the known fundamental interactions. 
Newton's theory of gravity is nowadays known to bo valid only in the special case of small velocities 
and masses, otherwise being replaced by Einstein's theory of general relativity. The relativistic theory 
of gravity has changed our understanding of the universe, and is a necessary ingredient in cosmology. 

Given the very different scales at which the theories of quantum mechanics and general relativity 
arc typically relevant, it would not be a problem that they have different natures, quantum and 
indctcrministic in the first instance, and classical and deterministic in the second instance. Since the 
ratio of the gravitational to electromagnetic force between two protons is of the order of 10^'^*, gravity 
has never appeared in microscopic physics experiments. However, quantum mechanics has taught us 
that in order to probe short distances we need high energy. Additionally, general relativity has taught 
us that energy gravitates. It should then occur that looking at shorter and shorter scales we would 
concentrate enough energy to make the gravitational field come into play. A simple calculation shows 
how important this would be if we would try to probe a distance of the order of the Planck length 
Lp ~ 1.6 X 10~^^cm, since doing that would concentrate enough energy to create a black hole. If 
we do not want to run into the paradox of a classical and a quantum theory coexisting at the same 
energy scale, close to the Planck scale gravity must be quantized too. 

The need for a quantum theory of gravity was advocated already at the dawn of the quantum 
revolution, but has since then revealed itself as one of the most challenging tasks for theoretical 
physicists, remaining unsettled up to now. 

The reasons for this situation are manifold. Already Heisenberg foresaw troubles because the 
coupling constant of gravity, Newton's constant G^v, has inverse-square energy dimensions Hence 
the effective coupling in a pcrturbative expansion is GnE^, and the perturbative theory is expected 

^Throughout this thesis I will use units h = c = 1, so that Gn = and the Planck units are all defined just in 

terms of Newton's constant, Lp = Tp = Mp^ = ^JU^,. 
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to break down when GnE"^ ~ 1, that is, at the Planck energy. In modern language we say that the 
theory is non-renormalizable. 

Furthermore, and besides the many technical difficulties ultimately related to the highly non-linear 
nature of general relativity, there are fundamental problems raised by the novel ideas that these two 
theories have revealed and which appear in contrast with each other. As already mentioned the 
classical nature of general relativity contrasts with the discovery that the world is indetcrministic. On 
the other hand, in quantum physics the role of spacetime as an inert stage upon which everything 
takes place is in contrast with the beautiful discovery that spacetime is itself a dynamical object. 

Physicists have embraced many different attitudes toward these problems. The most popular is to 
think that either one or the other of the two theories has to be modified. Most commonly it is thought 
that general relativity is just a low-energy effective theory for some new and more fundamental theory. 
It is also believed that the new theory should realize a unification with the other three fundamental 
interactions: the strong, weak and electromagnetic. The most popular candidate for such a theory 
seems to be string theory. Alternatively there are also proposals where quantum physics is derived 
from some new physics at the Planck scale. 

Another attitude is to be more conservative and try new ways of looking at the same theory. The 
point here is that maybe it is just the perturbation theory that is wrong and not the whole idea 
of quantizing general relativity as it stands. In such case a non-perturbative treatment would be 
necessary. Of course going non-perturbative is like opening Pandora's box, since most of the times 
in physics the best we are able to do is a perturbative expansion around one of the few cases we 
can solve exactly. Therefore the main problem in this approach is that of devising new tools for the 
non-perturbative quantization of general relativity. 

The models I describe and investigate in the present thesis take this second point of view as their 
starting point, and are known as Causal Dynamical Triangulations. In order to get a handle on non- 
perturbative quantization these models start from a basic observation, namely, that computability 
in general requires a finite number of degrees of freedom. This is something which is experienced in 
everyday life, for example whenever we look at the computer screen where all images are represented 
as a finite collection of pixels. In the case of quantum gravity the degrees of freedom are in the 
configurations of geometry at each spacetime point, and are therefore infinite. The proposal is to 
sample these degrees of freedom by replacing the continuous spacetime with a coarse-grained one 
built with identical triangles. Such a triangulation has to be dynamical, since otherwise it would not 
encode any degree of freedom. Finally causality is imposed, in a precise sense, as a physical ingredient. 

These models have been studied intensely in the last decade and have shown to lead to very 
interesting results. The triangulation of spacetime does what it is meant to do, in the sense that it 
produces a model which can be studied on a computer. Of course at some point one would like to 
make these triangulations finer and finer in order to recover the continuum theory. The computer 
simulations indicate that a continuum limit exists and that it has at least some of the good properties 
we expect it to have. The drawback to this approach is that we have basically no analytical method to 
investigate the properties of such models in 4-dimensional spacetime. Computer simulations are our 
only tool, and they are limited by the computer power. This is certainly an unsatisfactory situation. 

The status quo is radically different in two spacetime dimensions, where the model can be solved 
exactly in a variety of ways and the continuum limit taken analytically. This is maybe no surprise 
because in two dimensions things are a lot easier and we have many techniques for studying discrete 
systems. Nevertheless, exactly solvable models are not so abundant and it is a lucky case that 
causal dynamical triangulations are one of those. A possible objection to such models is that in 
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two dimensions there is no classical general relativity to start with, and that the quantum theory is 
renormalizable because Newton's constant is dimensionless. On the positive side, the two-dimensional 
model can be considered as an interesting playground for the development of tools and ideas, with 
the hope to generalize them at a later stage to the higher-dimensional case. 

Unfortunately any attempt to extend the analytical methods developed in two dimensions to higher 
dimensions has up to now failed. Not even in three spacetime dimensions are we able to solve the full 
model nor to compute any observables. 

The aim of this thesis, besides the non-marginal one of making the author deserve the PhD degree, 
is to try to take a small step beyond this impasse. 

After reviewing in Chapter 1 in more detail what has been briefly mentioned in this preface, in 
particular giving an introduction to the methods and the results of causal dynamical triangulations, 
we will proceed in the next two chapters to illustrate in detail two original pieces of work which both 
represent a small step away from the special case of pure gravity in two dimensions. 

In Chapter 2 we tackle the problem of coupling matter to causal dynamical triangulations in two 
dimensions. To fix a starting point we choose one of the simplest, and certainly the most notorious, 
matter fields on a lattice, the Ising model. We show how an old method of investigation for spin 
systems on a lattice, the high- and low-temperature expansion, can be applied on the dynamical 
lattice represented by the triangulation. The application is not trivial as it requires an ensemble 
average over the lattices for the graphs appearing in the expansion, but we will develop a scheme that 
allows the computation of such averages. The method is quite general and can in principle be used 
for other kinds of matter models. 

The work presented in Chapter 3 goes into the realm of higher dimensions. We introduce a 
particular model of three-dimensional causal dynamical triangulations and show that it can be solved 
in an asymptotic limit and a continuum quantum Hamiltonian is found. This is an important step 
in the understanding of these kinds of models, because this is the first case in which, for a dimension 
larger than two, a continuum limit has been obtained analytically, although many things remain to 
be understood and some approximations still removed. 

Together, the works presented in Chapter 2 and 3 constitute an advance in the analytical under- 
standing of causal dynamical triangulations and have the appealing feature of indicating directions 
for further development. 
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Chapter 1 



Introduction to Causal Dynamical 
Triangulations 



In this Chapter I will review the Causal Dynamical Triangulations (CDT) approach to quantum gravity. First 
I will briefly review the problems we face in quantizing gravity, then give a short introduction to the non- 
causal predecessor of CDT. In Sec. 1.3 I introduce the causal construction for the (l+l)-dimensional case, and 
proceed to review what is known about the higher-dimensional cases in the next two sections. 

Several excellent reviews exist on quantum gravity in general (see for example [2], or [3] for an historical 
overview) and on CDT in particular (see [4] and [5]). Here I will try to review all the necessary background in 
the attempt to make this thesis as self-contained as possible. In particular the first section sets the background 
of facts and ideas against which the research on non-porturbativc quantum gravity takes place, while the other 
sections recall the technical aspects needed for the two final chapters. 

1.1 The non-pert urbative quantum gravity program 

With the term quantum gravity I will refer to the quantum theory, to be constructed in one way or 
another, based on the classical Einstein-Hilbert action 



whore G and A are respectively the Newton's and the cosmological constants, M is the spacetime 
manifold, g^i, the spacetime metric, g its determinant and R the associated Ricci scalar curvature. 

One of the most powerful ways to quantize a classical field theory is the path integral formalism. 
For the case of gravity we would formally write something like 



where the integration should be over the (diffeomorphism-equivalence classes [g] of) metrics on space- 
time. Such a functional integration is rather formal and one still has to make sense of it in some 
way. 

Usually in quantum field theory we would add a source term in the action, something like the 
/ d'^xj{x)4>{x) for the scalar field, and use the path integral as the generating function of the Wightman 
functions, i.e. the vacuum expectation values of products of the field. Knowledge of the generating 




(1.1) 




(1.2) 
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function would then mean full knowledge of our theory. In general one can only hope to be able 
to compute it for the free theory, and switch to the interacting one by studying the perturbative 
expansion. In writing down the perturbative expansion for a Wightman function the path integral 
formalism turns out to be extremely useful. One can think of the path integral as a book-keeping device 
that tells us what to compute at every order. We would then have to actually do these computations, 
and usually would get divergent quantities. If we are lucky the theory is renormalizable and these 
divergences can be removed order by order by a suitable redefinition of a finite number of parameters. 
If the theory is not renormalizable, at every order there will be new parameters to readjust and the 
predictive power would be limited. For this reason wc consider non-renormalizable theories as non- 
fundamental, although they can still be of use as effective low-energy theories (as, for example, the 
Fermi theory of weak interaction). 

In applying the idea of perturbative quantization to gravity, one has to overcome a large number 
of technical difiiculties. Efforts towards the perturbative quantization of gravity have led to impor- 
tant techniques, ranging from the introduction of the "ghost" fields [1] to the development of the 
background field method [6], which have subsequently found a fundamental role in the quantization 
of Yang-Mills theories. The greatest difficulty is certainly that the coupling constant has the dimen- 
sion of inverse-square mass, which makes the theory non-renormalizable. By a dimensional counting 
argument we expect to have worse and worse divergences at higher orders of perturbation theory, 
and so new counterterms will be needed at each order. One could still hope that by some miracle 
the divergences could be reabsorbed, i.e. physical quantities be made finite just by field redefinitions, 
without having to introduce new parameters. This is exactly what happens at one-loop, as 't Hooft 
and Veltman [7] have shown. However the miracle did not last long, since Goroff and Sagnotti [8] 
showed that at two loops divergences appear that require a new counterterm, which cannot be re- 
absorbed in any way. Furthermore the dream of a finite perturbative quantum gravity was already 
broken in the presence of matter even at one loop [7] . 

The problem of non-renormalizable divergences has produced two main lines of thinking: one is 
to think of some other theory, which has general relativity only as a special limit and try to quantize 
the new theory (supergravity, strings, etc.); the other is to think that maybe the quantum theory of 
gravity still exists in a non-perturbative sense (asymptotic safety, canonical loop quantum gravity, 
etc.). The CDT model fits in this last perspective. 

As argued by Weinberg [9] , there is a well-defined scenario under which a non-perturbative theory 
of gravity could still make sense. Namely, if the renormalization group fiow of a perturbatively non- 
renormalizable theory would admit a non-Gaussian ultraviolet fixed point, it could still be considered 
a meaningful theory, in principle with the same predictive power as a renormalizable one. Such a 
scenario goes under the name of "asymptotic safety" , referring to the fact that such a theory would 
be safe from divergences as the cut-off is removed. 

Not many examples of such theories are known. The most relevant for our discussion is the case 
of gravity in 2 -|- e dimensions, for small e [9] . The big question is whether gravity is asymptotically 
safe also at e = 2. There is no answer to it yet, although some positive indications have come out 
in recent years (see [10] for a review), but definitely this idea is at the heart of all non-perturbative 
approaches to quantum gravity. 

Another great stimulus in the direction of non-perturbative quantization came from (2+l)-dimensional 
gravity. At a first look the theory seems to have all the problems of its higher-dimensional analogue, 
in particular since Newton's constant has dimension of length. Despite that, as Witten first showed 
in [11], the non-perturbative theory not only makes sense but is actually an exactly soluble system. 
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Of course a more careful look shows that the two theories have a fundamental difference, with the 
lower-dimensional one having no local degrees of freedom, once the constraints have been enforced. 
Because of this the quantum theory can be made meaningful in several ways, for example, by working 
in reduced phase space or imposing the constraints after quantization (sec; [12] for a review). There 
are still several difficulties, for example, with higher topologies or with the path integral over metrics, 
but for sure it is encouraging. 

A lot of effort has been spent over the years to build a meaningful non-perturbative theory of 
quantum gravity in four spacetime dimensions, but unfortunately without conclusive results yet. One 
reason for this situation is that already in "easier" theories a non-perturbative study of the path 
integral is an extremely difficult task. 

In general the path integral has no rigorous mathematical definition. For the non-rclativistic 
free quantum particle, the analogue of (1.2) can be made meaningful by rotating time to imaginary 
time. The quantum theory then becomes a stochastic theory and the integral has a measure-theoretic 
definition in terms of Wiener measure, whose properties are known. For example, its support is given 
by everywhere non-diffcrcntiablc paths. The Wiener measure is usually obtained by discrctizing the 
time line, which turns the path integral into a finite product of ordinary integrals, and then by taking 
the continuum limit. 

In quantum field theory the situation is mathematically less safe. For the Wick-rotated free scalar 
theory it is known that the support of the theory is on distributions. For the interacting theory, no 
rigorous result is known, but one expects troubles because the fields will be at least distributions, 
whose products are ill defined. This is connected to the problem of ultraviolet divergences in the 
perturbative approach, which in turn is solved by renormalization if the theory is renormalizable. If 
we want to address non-perturbative questions, like the confinement in QCD, or if - like for the case of 
gravity - the theory is non-renormalizable and we want to explore the asymptotic safety hypothesis, 
we have to find a way to define the path integral. What is usually done is to discretize the spacetime 
on a lattice in analogy with the time variable of the non-relativistic case. The Wick-rotated theory 
on the lattice is a well-defined object which we can put on a computer to make calculations. It also 
has the form of a partition function for a statistical mechanical system, making it possible to borrow 
methods from statistical mechanics, like the renormalization group. The quantum field theory should 
be defined as the continuum limit. This limit, where the short-distance cut-off a is sent to zero, should 
be such as to keep the physical mass finite. In the statistical mechanical language this means we have 
to go to the regime of infinite correlation length ^ = {mphysO')~^ , i-e. we have to tune the theory to 
the critical surface. 

Clearly the task of constructing a non-perturbative quantization of gravity along these lines is 
extremely challenging, and even at a quick look already presents a considerable number of difficulties. 
I list some of the most relevant obstacles one encounters in trying to construct a non-perturbative 
theory of gravity: 

(i) no notion of preferred time is present in general relativity, so it is not possible to perform a 
standard Wick rotation; 

(ii) we would like to have a regularization that preserves the diffeomorphism invariance of general 
relativity; this may be not a must, as also ordinary lattice field theory breaks Lorentz symmetry, 
but then one should recover it in the continuum; 

(iii) as in gauge theory we expect divergences from the gauge volume (the infinite volume of diffeo- 
morphism group here); a good gauge fixing would be needed and the associated Faddeev-Popov 
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determinant should be computed; alternatively, explicit field coordinates on the physical con- 
figuration space '^D^ff^M)^^ diffeomorphism-equivalence classes of metrics [5^^] (the so-called 
geometries) should be found; 

(iv) in case wc manage to define a Wick rotation or in case wc start with the Euclidean theory from 
the beginning, we still have to face the problem that the action is unbounded from below, due 
to the sign of the kinetic term for the conformal mode. 

Although since the pioneering work of Leutwyler [13] on a lattice path integral formulation for 
gravity many things have been learned and understood, these issues still remain an open problem. 

An old proposal by Hawking to deal with (i) is to start with an Euclidean theory from the beginning 
[14] , from which physical information can directly be recovered, like for example on the ground state 
of the theory. The hope was also that a solution to the problem of rotating back to Minkowskian 
signature would be found once the Euclidean theory is solved. 

The problem (iv) was also addressed by Hawking and others in the context of formal Euclidean 
path integrals with the proposal to analytically continue the conformal factor to imaginary value [15]. 
A more promising and less ad-hoc solution seems instead the idea that the cure would be naturally 
given by a non-perturbative treatment of the path integral. In this case, non-trivial measure contri- 
butions arising from Faddeev-Popov determinants can in principle cancel the conformal divergence, 
as indicated by some formal calculations [17, 18], inspired by the generally covariant method used by 
Polyakov for the bosonic string [16]. 

Several discretized models have been proposed over time, see [19] for an overview. One of the 
most studied has been probably the quantum Regge calculus, because it is close to the spirit of 
general relativity, coordinate independent, and can be used classically to approximate any geometry. 
A path integral based on it presents several problems. In order to avoid some of these the Dynamical 
Triangulation model was proposed in the '80s, and eventually led to the CDT model at the end of the 
'90s. 

The CDT model is a very concrete and promising way to face the above-mentioned problems, 
with a coordinate-independent discrete path integral over geometries admitting a Wick rotation. I 
will introduce it in this chapter, after having recalled the necessary ingredients borrowed from its 
predecessors. 



1.2 The Dynamical Triangulations approach 

It is well known that Feynman's path integral for quantum mechanics can be defined as a limit for 
N going to infinity of N ordinary integrals for the position at time = k-^ (for /s = 1, A^) of a 
particle during its path from an initial to a final point in a total time t. At finite N only a finite 
number of points, those at time tk, in the path are specified and the path at intermediate times can 
be thought of as straight. This is then called a piecewise linear path. The limit as A'' —> 00 is the 
continuum limit. 

The idea of dynamical triangulations for the quantization of gravity is simply to generalize this 
procedure to higher dimensions. 

The analogue for spacetime geometries of the piecewise linear paths of the Feynman path integral 
is given by piecewise linear manifolds, which I am going to recall in the following subsection. 
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Figure 1.1: A piecewise linear path appearing in the definition of the Feynman path integral. 

1.2.1 Preliminciries: simplicial manifolds and Regge calculus 

A piecewise linear n-dimensional manifold is a collection of n-dimensional (flat) poljd;opes (higher- 
dimensional analogue of polygons (n = 2) and polyhedra {n = 3)) glued together along their (n — 1)- 
dimensional faces in such a way as to preserve the topological dimension. In general it will not be 
possible to embed the object obtained in this way isometrically in M", meaning that there will be 
curvature defects. 

To be more specific and to simplify issues, piecewise linear manifolds are generally taken to be 
simplicial manifolds, i.e. the constituent polytopes are chosen to be simplices-'^. 

A simplex is the higher-dimensional generalization of a triangle, a polytope with the minimal 
number of faces in any given dimension. For example, in three dimensions it is a tetrahedron. In 
general, a n-dimensional simplex C7„ can be defined by its embedding in flat space as the convex hull 
of n -|- 1 afRnely independent points (meaning that no more than m -|- 1 of these points are in the same 
m-dimensional plane) of IR*^, for d>n. Denoting these points by vi, ...,Vn+i, the simplex is given by 

{n+l n+1 ^ 

xgW^, d>n\x = ^ XiVi ; Ai > 0, ^ Ai = 1 I . (1.3) 
i=l j=l J 

A subsimplex of C7„ is a simplex ak (for k <n) whose vertices are a subset of the vertices vt of cr„. 
Subsimplices of dimension n — 1 are called faces, while those of dimension n — 2 are called bones or 
hinges. Regardless of n we will call (Ti also edges or links. 

A collection of simplices glued along their subsimplices is called a simplicial complex. Finally, a 
simplicial manifold is a simplicial complex in which the neighbourhood of any vertex, i.e. the set of 
simplices sharing that same vertex, is homeomorphic to B„, an n-dimensional ball in R". We will 
construct simplicial manifolds by gluing together no more than two simplices along the same face 
(faces l)(4()ii,2,iiig to only ouc simplex will belong to the boundary of the simplicial manifold). Even 

'^Onc advantage of simplices is that their geometry is completely determined by the specification of their edge lengths. 
Appealing to "universality" , this specific choice of building block should not affect the continuum results. There is very 
good evidence of universality in two-dimensional DT models, but the property needs to be verified for individual models. 
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with this restriction, gluing together simpHces at random can give rise to singular points, with a 
neighbourhood not homeomorphic to -B„. For example, in three dimensions such a neighbourhood 
may have a higher-genus surface as its boundary. We will not be interested in this issue in what 
follows, since our models by construction are given in terms of simplicial manifolds. 
It is useful to introduce the so-called^ f-vector of the simplicial manifold M, 

f{M) = {No,N,,...,Nn} , (1.4) 

where Nk is the number of /c-dimensional (sub-)simplices in M. If the manifold has a boundary DM 
we also introduce the /-vector of the boundary, 

f{dM) = {N^,N^,...,Nt,} . (1.5) 

The elements of the /-vectors are not all independent, but have to satisfy topological relations. The 
most famous one is Euler's relation 

n 

^{-irN,=x{M) , (1.6) 

fe=0 

where xi^) is the Euler characteristic of the manifold, which for n = 2 is x(M) = 2 — 2g — b, with 
g being the number of handles (the genus) and b being the number of boundaries, and for n odd is 

given by x(M) = x{dM)/2. 

Besides Euler's relation also the Dehn-Sommerville relations must be satisfied (see [21] for details), 
namely, 

Nl = (1 - {-iT-^)N, + i-iT-'-' K X J • (1-7) 

j=fe+l \ / 

Because Euler's relation is independent of the Dehn-Sommerville relations only in even dimension, we 
find that among f{M) and f{dM) only n clement arc free if n is even, and n+1 if n is odd (respectively 
n/2 and (n + l)/2 if there is no boundary, since then only [(n + l)/2] of the Dehn-Sommerville relations 
are linearly independent). 

We can fix the metric properties of a simplicial manifold by fixing its edge lengths, and remembering 
that the inside of the simplices is flat, it is trivial to see how this fixes distances between points. 
Furthermore, the notion of curvature can be implemented too, as first shown by Regge in [22], and 
later made more rigorous in [23] . 

The idea is to replace the continuum notion of curvature in the following way: 



/ 



(Px^R^2 J2 ^hSh, (1.8) 

J M 

where is the deficit angle at h, 



^ hinges h 



Sh = 27T- 0,„^h (1.9) 

CTnDh 

(the sum is over all n-dimensional simplices cr„ sharing the hinge h, and 9ar,>h is the angle between 
the two faces of (j„ sharing h). Fig. 1.2 illustrates the two-dimensional case. 

The formula (1.8) can be understood by recalling that intrinsic curvature (as encoded in the 
Ricmann tensor) is associated with the variation a vector undergoes when parallel transported around 



^The name derives from the fact that in mathematics its elements are usually denoted /j, rather than Ni as we are 
used to in physics. 
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Fi gure 1.2: An example of a positive deficit angle. 



an infinitesimal closed loop spanned by two vectors. One can then easily realize that the curvature is 
zero inside a simplex (since we have taken it to be flat from the beginning) and that it is also zero at 
all points lying inside the faces (as Fig. 1.2 illustrates at the faces there is only extrinsic curvature). 

The curvature is therefore concentrated at the hinges, and one may think of (1.8) as saying that 
the Ricci scalar curvature is given by a linear combination of delta functions with support on the 
hinges. This fact should not constitute a worry, it is reminiscent of the distributional character of 
fields in quantum field theory. 

As the title in [22] suggests the original motivation behind the introduction of simplicial manifolds 
in physics was to have a coordinate-independent approximation scheme for classical General Relativity. 
With the notion of curvature at hand, one can introduce the discrete analogue of the Einstein-Hilbert 
action (1.1), called the Regge action^. 



vary it with respect to the edge lengths, and obtain a simplicial analogue of the Einstein equations 
which is explicitly independent of coordinates. 

Although this discrote formulation of gravity has had a number of classical applications, it received 
a boost of interest in the eighties with the proposal of using it to define a path integral quantization 
of gravity [26, 27]. 

The prescription for defining the path integral (1.2) from Regge calculus is to take for M a piecewise 
linear manifold with fixed connectivity, and to integrate over the edge lengths, with the exponential 

of the Regge action as weight. 

•^In case the manifold has a boundary, boundary terms must be added in the action. They are essential in the quantum 
theory in order that the transition amplitudes satisiy the correct composition law and have the correct classical limit, as 
was pointed out in [24]. The full continuum action is fixed by the requirements of additivity {i.e. given two contiguous 
regions with metrics g and g' the action for the union of the two regions must be S[g + g'] = S[g] + S[g']) and of 
equivalence of the classical equation of motions with the stationarity of the action under variation of the metric at 
fixed boundary geometry. It is the standard Einstein-Hilbert action plus the so-called Gibbons-Hawking term for the 
boundary, 



where h is the determinant of the induced metric on the boundary, and K is the associated extrinsic curvature. For the 
case of a simplicial manifold, the corresponding discrete action was given in [25]. It amounts to adding the following 
term to the Regge action: 




(1.12) 




(1.10) 



(1.11) 



hinges h 



where = - 
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Quantum Regge calculus has usually been constructed and studied for Euclidean signature, even if 
a Lorentzian version has also been attempted [28] . It was shown to reproduce the standard continuum 
theory in the weak field limit but understanding its non-perturbative properties remains an open 
problem. 

Many choices are possible for the measure, a common choice being of the type 



and there is no agreement on the correct one. Furthermore it has also been strongly argued that a 
gauge-fixing and related Faddeev-Popov determinant are needed [29] , but in this case (at least in two 

dimensions) the measure becomes non-local, making even computer simulations difficult. 

A weakness of quantum Regge calculus is that no analytic treatment is possible, not even in two 
dimensions. 

1.2.2 Dynamical Triangulations 

One of the potential problems of quantum Regge calculus is the use of a fixed simplicial manifold. 
This may prevent the path integration from exploring entire regions of the space of metrics where the 
metric is very different from those realizable for the chosen connectivity of the triangulated manifold. 
On the other hand, other metrics are overcounted and would need a gauge fixing (think for example 
of the flat two-dimensional triangulation and how many ways we have to (over-)count it by varying 
the edge lengths in such a way that the vertices remain in the flat plane). 

A possible solution to the first problem is to include also a sum over simplicial manifolds. The 
idea of having variable connectivity and variable edge lengths was revived recently under the name 
of random Regge calculus [30], but apart from not solving the other problems of Regge calculus, it is 
extremely difficult to handle. 

To simultaneously solve the various problems and introduce a simplification, it was proposed in the 
mid-eighties to fix completely the edge lengths to be all equal and only keep the sum over simplicial 
manifolds [31]. This is the approach we call Dynamical Triangulations (DT), and it is thoroughly 
reviewed in the book [20]. 

Contrary to Regge calculus, this formalism is not well suited for classical approximations, but 
is specifically designed as a functional integral over geometries. It is expected that the ensemble of 
dynamical triangulations is more evenly distributed on the space of all geometries than that of the 
quantTim Regge calculus on a given simplicial manifold [32]. 

Another advantage of dynamical triangulations is that it introduces a natural UV cut-off a, the 
edge length common to all simplices. This provides a regularization of the path integral and allows 
for a clear definition of the continuum limit procedure, which I will explain after having given some 
more details. 

In this approach the integral over geometries is turned into a discrete sum over simplicial manifolds, 



n 



edges e 




(1.13) 




r 



(1.14) 



Furthermore, due to the restriction to identical building blocks, the action simplifies a lot, basically 
reducing to a counting of global variables. The volumes and internal angles of the simplices can only 
assume certain discrete values, and it is an easy exercise to check that the Regge action for Euclidean 



1.2. THE DYNAMICAL TRIANGULATIONS APPROACH 



13 



gravity on a closed manifold becomes 

Sr = KnNn - Kn-2N„-2 , (1.15) 

where we have introduced the dimensionless couplings 

1 , , arccos - 
Hn = AoV;„ + -n(n + 1) 

2 ISttGo (iig) 



'«n-2 = 



8 Go 

where Aq and Go are the bare cosmological and gravitational constants, and where we have used the 
fact that each simplex has n{n + l)/2 internal dihedral angles equal to arccos^. The path integral 
(1.2) is reduced to a combinatorial object, which takes the form of a generating function for the 
number of simplicial manifolds with given numbers of simplices and hinges: 

Z{Kn-2, K„) = E c^^"^"^^^ = E ^e--^"+— , (1.17) 

where wc have also taken into account C (T) , the order of the automorphism group of the triangulation 
T, to avoid overcounting of manifolds with special symmetries'*. 

Due to the existing relations among the Ni, the action is now bounded above and below for fixed 
volume Nn (to put a rough bound, note that since there are n{n+ l)/2 hinges per simplex, it follows 
that Nn-2 < ^^^^A^n). If we rewrite (1.17) as 

Z{Kn-2, «„) = V e-«"^" V -l-e'^"-^"- = V e-«"^"Z(«„_2, iV„) , (1.18) 

we see that Z(k„_2, A^„) is only a finite sum with respect to Nn-2, and so is well defined. In order 
to check whether Z(«;„_2,k;„) is well defined too, we need to check the convergence with respect to 
the sum over N^. Clearly (1.18) will converge for some values of the coupling constants if and only if 
Z{Kn-27 Nn) grows at most exponentially, like Z(k„_2, Nn) ^ e^^", for come constant c. A necessary 
condition for this to be true is that the topology must be fixed, otherwise Z{Kn-2, K,n) would grow 
factorially^. That this condition is also sufficient is proven only for n = 2, but numerical results 
support it also in higher dimension [33]. The existence of an exponential bound is now generally 
believed to be true. 

At this point it is worth noting once more that the regularized Euclidean version of the path 
integral, defined in (1.17), is a partition function for a model of statistical mechanics. I will from now 
on often refer to it and its model-dependent variations as partition functions. The relation between 
Z{Kn-2, ^n) and Z{Kn-2, i^n) is Seen to be that between a canonical and a grand canonical formulation. 
Because in general relativity the volume of spacetime is not fixed we are usually interested in the grand 
canonical partition function. 



^This factor can be understood in the following way. If we would label all the subsimplices of the triangulation and do 
the sum over labelled triangulations, we would have to fax;tor out the number of ways a triangulation can be labelled. In 
general this number is A'^o!A''i!...7V„!, but for triangulations with some symmetry it is less, precisely No\Ni\...Nn,\/C{T). 
The factor 1/C(T) thus survives in the case where we sum over unlabcUcd triangulations, and therefore can be thought 
of as what remains when factoring out the volume of the difEcomorphism group. 

®In two dimensions, a sum over topologies can still be defined after having first performed the sum for any fixed 
genus, in a kind of resummation of the series, but only in a double-scaling limit where the weight of higher-genus 
sums goes to zero like a specific power of the edge length a. Since this procedure is neither unique nor seems to be 
generalizable to higher dimensions (where topologies are not even classified), topology is usually fixed in models of 
dynamical triangulations. We will also stick to a fixed topology in this thesis. 
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Going back to the issue of convergence, we will now see that it is closely related to the way in which 
we can formulate the continuum limit. For fixed topology, if we assume for the canonical partition 
function a large- volume behaviour of the form 

Z{Kn-2, N,,) ^ iv^^^-^^-^e^^"-^)^" X (1 + 0(1/7V„)) , (1-19) 
the partition function is convergent for k„ > K!^{Kn-2), and near the critical line (k„_2) we get^ 

Z{Kn-2, K„) ~ {Kn - <(/^„-2))'-^^"''-^^ • (1.20) 

Since k„ is conjugate to the number of n-simplices, we can extract the average lattice volume in the 
ensemble of triangulations from the partition function as 

^^^^^_ ainZ(.„_2,.„) _ (^21) 
Prom (1.20) it is then clear that the average number of n-simplices diverges at the critical line like 

/AT \ 7(«^»-2) - 2 

{^n) T , (1-22) 

l^n l^n\^n—2j 

and it is exactly there that wc can send the cut-off a to zero, while keeping the physical volume of the 
universe finite. More precisely, we have to approach the critical line according to a proper scaling, like 

~<(Kn-2) + Aa" , (1.23) 

so that 

(y)=a"(jV„)~ ^^^-^^~^ , (1.24) 

where A is the renormalized cosmological constant, conjugate to the volume in the continuum. 

Tuning k„ to obtain a finite continuum volume may not be enough to arrive at an interesting 
continuum limit. We may still need to tune also k;„_2 to some specific value, for example, in order 
to get a divergent (lattice) correlation length between local variables^. In two dimensions, taking the 
limit of infinite lattice volume turns out to be sufficient, leading automatically to a continuum theory 
of two-dimensional quantum gravity. Of course n = 2 is a special case since the gravitational constant 
couples to the Eulcr characteristic and thus plays no role when the topology is fixed. The situation 
in higher dimension is not completely clear, numerical simulations for n = 4 indicate that different 
phases are separated by first-order transitions [34, 35] so that either a non-trivial continuum limit 
(one with propagating degrees of freedom) is automatically reached with the infinite- volume limit or 
otherwise there is no chance to get one. 

In the infinite-volume limit one finds for n = A two phases, a crumpled one (at K2 < K2) with a very 
large Hausdorff dim,ension^ dn (~ 00) and an elongated phase (at K2 > K2) with the characteristics 



''This expression also defines the susceptibility exponent 7(Kn— 2)- The name comes from the fact that it is the 
exponent of the non-analytic part of the second derivative of tlic partition function, which is usually called susceptibility. 

classical example in this respect is the Ising model, where the limit of infinite lattice volume and zero lattice 
spacing is not sufficient to obtain a continuum field theory with local propagating degrees of freedom. This is only 
achieved by also fine-tuning the temperature to a critical value corresponding to a second-order phase transition. 

*The Hausdorff dimension is defined as the leading power in R with which the volume of a ball B(R) of radius R 
scales with the radius. More precisely, B{R) is defined as the set of points x at geodesic distance d{x, xo) < R from 
a given point xq; taking first the manifold average of the volume of B{R) with respect to the center xq, and then the 
ensemble average over the triangulations we define the Hausdorff dimension dn as 

{V{R)) ~ (1.25) 

in the limit of large R. If {V{R)) grows faster than any power of R, we set dn = 00. 
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of a branched polymer and du ~ 2. Since as just mentioned the transition between the two is only 
first-order nowhere a 4-dimensional spacetime can be found, and the model does not seem to have 
good chances of describing our world. 




Figure 1.3: The phase diagram of 3- and 4-dimensional Euclidean dynamical triangulations. The region of convergence, 
from which the critical line is approached, is above the critical line. 



1.2.3 Dual complex and matrix models 

Before moving to the causal version of dynamical triangulations, which has proven to be better behaved 
and which is the main topic of this thesis, it is useful to recall some of the techniques used in the 
investigations on dynamical triangulations. 

Something that wc will use very often in this thesis is the dual complex of a triangulation, which en- 
codes all the combinatorial information wc need. The dual complex is associated with a n-dimensional 
simplicial manifold via a dual mapping, in the sense that the dual complex of the dual complex is 
again the original simplicial manifold. The dual mapping associates to any fc-dimensional subsimplex 
(Tfe a (n — A;)-dimensional dual polytope (Jn-k whose sub-polytopes are dual to higher- dimensional 
subsimpliccs. In particular, the dual to a n-simplex is a vertex, the dual to a face cr„_i is an edge, 
and the dual to a hinge it„_2 is a polygon. 




Figure 1.4: A portion of a triangulation and the corresponding dual graph. 
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If in a triangulation all the building blocks are of the same kind, but its vertices have all possible 
coordination numbers, the situation in the dual complex is reversed, namely, all kind of polytopes 
appear but all the vertices have coordination n + 1. This turns out to be a useful property in many 
situations. In two dimensions, it allows us to regard the dual of a triangulation as the Feynman 
diagram of a (7!)'^-thcory. This is precisely the link between two-dimensional quantum gravity and 
matrix models, on which a lot of literature exists (see, for example, [36]). 

Matrix models are a powerful method, which enables us to solve exactly the two-dimensional DT 
model for the pure gravity case [37] , and for gravity coupled to matter in the form of the Ising model 
[38], Potts model [39] or hard particles [40]. The partition function for the DT model is identified 
with the generating functional for the connected diagrams of the corresponding matrix model (in the 
large- A'' limit), 

Z{X, ...) = Inl d'^'A d^'i?...e-^^^(^'^'-) , (1.26) 
where for example in the pure case the potential is given by 

V{A) = \a^ - , (1.27) 

with g = e~^. 

A generalization of matrix models for higher-dimensional DT has been constructed, but it has not 
led to any advance because of major problems like the fact that there is no equivalent of the large- A'' 
limit, it is not possible to restrict to manifolds, and there is no reduction to eigenvalues. 

In two dimensions, other methods are possible but it is very difficult to generalize any of them to 
higher dimensions, where the most valuable method remains that of Monte Carlo simulations. 



1.3 Enforcing Causality: the (l+l)-dimensional case 
1.3.1 The model 

As I have argued in the previous section, the naive idea of gluing together simplicial building blocks 
without any restriction in general does not work. To begin with one may not get a manifold at 
all. So wc have to^ restrict to gluings that have the properties of a simplicial manifold. Also this is 
not enough, and if wc want to meaningfully define a partition function we have to fix the topology. 
Even if this makes the model well defined, the statistical ensemble is dominated by very pathological 
geometries with an effective dimension which is too large or too small. One possible reason for this is 
that the class of geometries being summed over is still too broad. 

A clear example of this point is the n = 2 DT model. There the effective dimension djy = 4 can 
be explained in terms of baby universes. A baby universe is a triangulated manifold with a loop of 
minimal length (of the order of three edges) as its boundary, along which it is glued to a mother 
universe (see Fig. 1.5a). When gluing triangles without any restriction apart from the fixed topology, 
it turns out that these baby universes dominate in the partition function. In the continuum limit we 
roughly speaking have a small two-dimensional universe attached to each point of the two-dimensional 
mother universe, giving the effective dimension 4. An artist's impression of this situation is given in 
Fig. 1.5b. 

^The necessity of this point may be argued. In the Group Field Theory approach to quantum gravity the appearance 
of conical singularities does not seem to be considered a major problem [41]. 
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(a) (b) 

Figure 1.5: (a): A portion of a triangulation with a baby universe appearing inside a loop made of three links. It 

should be understood that, according to the DT construction, all the triangles should be of equal size. When trying to 
draw a baby universe on a plane like here we are forced to draw smaller and smaller triangles, generating some kind of 
fra<;tal picture, (b) : A pictorial representation of baby universes attached everywhere on the mother universe. 

These baby universes cannot be removed, for example, by imposing that no more than one triangle 
can be contained inside a three-edge loop of the mother universe, because then baby universes will 
simply appear on four-edge loops and so on. It turns out that they can instead be removed, at least 
in two dimensions, by introducing a new principle in the construction of the geometries included in 
the sum. 

In [42] Ambj0rn and Loll introduced such a principle and thereby defined the model of Causal 
Dynamical Triangulations (CDT). The idea is to restrict the path integral to geometries with a well- 
defined causal structure. 

Although at the classical level causality is considered a necessary condition for a theory to be 
physical (for example, a spacctime with pathologies in the causal structure, like the Godel spacetime, 
is usually considered unphysical), this condition has often been put aside in the construction of a 
quantum theory of gravity, especially in its Euclidean formulation. Nevertheless the importance of 
causality has been stressed by many, in particular by Tcitclboim [43] in the path integral approach, 
and has led to the formulation of new approaches like that of Causal Sets [44] . 

The first step in implementing causality into dynamical triangulations is to work with Lorcntzian 
spacetime geometries from the outset, because otherwise it would make no sense to talk of causality. 
To do that in simplicial geometry amounts to using flat simplices with Lorentzian signature. I illustrate 
this first in two dimensions, postponing the cases n = 3, 4 to the following sections. A triangle with 
Lorentzian signature is roughly speaking a triangle cut out of (l-l-l)-dimensional Minkowski space. 
Regge calculus can be defined on such pseudo-Riemannian structures, as has been known for a long 
time [45]. In a dynamical triangulations approach we want to take these triangles to be all equal to 
each other. A convenient choice is therefore to have two edges of time-like type and a space-like one, 
say, of squared length /| = — aa^ (with a > introduced for convenience) and = c? respectively, and 
where gluing is only allowed between edges of the same type. However, this signature implementation 
is not enough to implement causality, since, for example, time-like lines forming closed loops are still 
present. Neither is it in general enough to change the properties of the model with respect to the 
Euclidean formulation. An example of such a model can be formulated [46] as a matrix model with a 
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matrix A associated to space-like edges and a matrix B associated to time-like edges, and a potential 

(1.28) 

The integration over A can be performed explicitly, because it is just a Gaussian and we obtain 

(1.29) 

which describes a DT model with sqiiarcs instead of triangles (the squares are given by pairs of 
triangles glued along the space- like edge) , and which is known to lie in the same universality class as 
the original Euclidean one with triangles [37]. 

To enforce causality and to select geometries in a different way, the following construction, for the 
two-dimensional case, was introduced in [42]. A triangulation is made of t "strips", a strip being a 
triangulation of topology E x [0, 1] (with S chosen to be either or [0, 1] once and for all) where 
the space-like edges belong either to the initial boundary S x {0} or to the final boundary S x {1}. 
Different strips are glued along these boundaries (with the constraint that the lengths of the glued 
boundaries must match), and the topology of the full triangulation is again S x [0, 1] (or E x 5^ if we 
impose periodicity in time). An example of such a construction is shown in Fig. 1.6, together with 
the dual graph, which is obtained with the same rules as explained in the previous section. 



U = Tr 



2 2 



gAB' 



U = Tv 



Ib' - 3-B^ 
2 2 



Figure 1.6: Example of a triangulated piece of spacetime in the CDT model (left), and the corresponding dual graph 

(right). The picture captures the way triangles are glued together, but does not represent faithfully intrinsic distances, 
and therefore the curvature properties of the simplicial geometry. If we were to respect this property, we would not be 
able to draw the picture in a plane, because the model's triangles are all equal. 

It turns out that this new principle for the selection of geometries is relevant for the critical 
properties of the model. Thanks to it the baby universes are eliminated, and the effective dimension 
of the geometries is = 2 [42]. As I will review in the following section, also in higher dimensions 
it is found that dn = n for n-dimensional CDT [n = 3,4). 

Two features are extremely important in this construction: there is a global time t, with respect 
to which there is a foliation into space plus time, and there is no topology change of the slices with 
time. It is thanks to these features that the model belongs to a different universality class than the 
Euclidean one. As was shown in [42] , if one lifts these constraints by allowing the slices to branch from 

into several S^''s and forcing all of them but one to eventually contract to a point and disappear 
without rejoining with any of the others (in this way no handles are formed and the topology of the 
total manifold is unchanged, see Fig. 1.7), these branchings dominate and the model ends up again in 
the same universality class of the Euclidean DT. 
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Figure 1.7; The branching of a universe and the appearance of a baby universe as seen from a time-slicing point of 
view. 

Having formulated the model with Lorentzian signature, one may wonder whether we end up 

having problems in defining the path integral, because of the oscillatory behaviour of the integrand. 
It turns out that instead in our discrete framework we can easily go back and forth from Lorentzian to 
Euclidean by simultaneously changing the sign of all squared time- like lengths. In order to do that, we 
simply have to change the sign of the parameter a, which amounts to an analytical continuation of it 
in the complex plane. A detailed analysis of the dependence of the geometric variables, and thus of the 
action, on the parameter a goes beyond the scope of this thesis, where we are mainly concerned with 
applying analytical methods after having made the continuation to Euclidean signature. It suffices 
to say that under the analytical continuation to negative a, the weights appearing in the sum over 
geometries transform according to 

e^s^R- ^ e-s--""* . (1.30) 

For details see [47]. 

It should be emphasized that this "Wick rotation" is well defined because we are in a discrete setting 
and because we have introduced a distinction between time- and space-like distances. The issue of 
how to perform the inverse Wick rotation comes back once the continuum limit is performed. In the 
two-dimensional case, it is possible to perform a standard analytic continuation in time T —iT and 
obtain a unitary theory [48] . 

In the rest of this thesis we will only use Euclidean signature. 

1.3.2 The solution 

In two spacetime dimensions, the total curvature is a topological invariant proportional to the Euler 
characteristic, and thus we can omit it from the path integral when working at fixed topology. The 
partition function then takes the simple form 




(1.31) 



where now the triangulations T/Va,* (with fixed number of triangles N2 and of given extension t in 
time) have to be constructed according to the causal rules just explained. 



20 



CHAPTER 1. INTRODUCTION TO CAUSAL DYNAMICAL TRI ANGULATIONS 



The model was first solved in [42], but it can be solved in various other ways too, including 
generalizations, like the inclusion of additional weights or of different building blocks (see [49, 50, 51]). 
This indicates a rich integrability structure and supports the universality hypothesis for these models. 
I will not review them all, but rather will recall the method and results we will need in the developments 

of Chapter 2 and 3. 

In [50] a formula was proven expressing the partition function of the (l+l)-dimensional CDT 
model, in its dual formulation and for open boundary conditions in the space direction, as the inverse 
of that of a hard-dimer model in one dimension: 

^*(- = e--) = ^ (1-32) 

where 

hard diincr config. D 

and where is the weight assigned to each triangle and t is the number of time steps. I rederive the 
proof of the inversion formula in appendix A-II for the three-dimensional case, which we will study in 
Chapter 3. 

This result is not accidental, it is quite general and indeed could already be found in [90] , in a more 
general mathematical context, where the notion of heap of pieces was introduced. Roughly speaking 
a heap of pieces is a partially ordered set whose elements can be seen as occupying columns out of a 
finite set of columns and where direct order relations only exist between elements in the same column 
or in neighbouring ones. It is not difficult to realize that the structure of the (l+l)-dimensional CDT 
model is (in its dual picture) that of a heap of dimers^'^. Stated differently this simply means that in 
the dual picture of Fig. 1.6 we arc free to slide the vertical segments horizontally, with the constraint 
that segments in the same slice or in neighbouring ones can neither cross nor touch each other. 

One can also give to the partition function (1.33) a transfer matrix formulation. This is done 
by introducing a two-dimensional vector space associated to each site, with base given so that (1,0) 
corresponds to the empty state and (0,1) to a dimer state. If we then associate a weight 1 to the 
transition empty-empty, a weight ^/u to the transition dimer-empty or empty-dimcr, and a weight 
to the transition dimcr-dimer, we can write the transfer matrix M around a site. Since the dimer 
partition function is given, for the case of periodic boundary conditions, by TrM* we can rewrite 
(1.32) as: 

Ztiu) = . (1.34) 




Since Tr{M*) = X\ + A*_, where 

A± = (l±\/l-4u)/2 (1.35) 
are the eigenvalues of M (with ^/u replaced by iy^ as in (1.34)), we have for large t 

Z^iu) ~ ^ (1.36) 



^"Note that a dimer is the dual of two triangles glued along their space-like edge, that is, it is the analogue of the 
4- vertex in (1.29). Observe that in order to describe the CDT model in terms of a matrix model, one would have to 
impose on it the requirement that in the generated graphs the vertices form a heap of pieces. It is not yet known how 
to implement such a constraint in a matrix model. 
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which is not analytical at m = 1/4, the critical point. 

The critical point can already be found by analyzing the one-step propagator. We know that in 
(1+1) dimensions'^ 

G{h,h;t = l\X)= Yl T^e-^"" (1-37) 



reduces to 



G(^i,Z2;i= 1|A) =e-^('^+'^) 1, (1-38) 

and this can be obtained by the generating function 

li,l2>0 strip-Ti^j^ 

which in turn can now be written with the use of the inversion formula (and considering x and y as 
the square roots of the weights of the boundary dimers) as'^ 

G(x,y;t=l\u = e~'^^) = 7 r — = -i- -. (1.41) 

'V J [ix) 

We see that if we put x = y = 1 {i.e. considering the partition function for At = 1 with free 
boundary conditions on the slices) we get a singularity in = |. This is exactly the critical point, 
and the continuum limit has to be taken by fine-tuning A to its critical value Ac = In 2 as explained 
in the following subsection. 

1.3.3 The continuum limit 

As I have already argued, in order to remove regularization artifacts and obtain a continuum theory 
we must tune the coupling constants to a critical value as a goes to zero. In two dimensions we only 
have one coupling constant, the cosmological constant, which is conjugate to the volume. Tuning 
it to its critical value achieves the infinite-volume limit and as I will show also leads to a divergent 
correlation length for the dynamical degrees of freedom. 

From (1.36) it is easy to see that in the limit A ^ Ac = In 2 (equivalent to u ^ Uc = 1/4) the 
volume diverges like 

If we scale the cosmological constant like A ~ Ac -I- ^Aa^, with A the renormalized cosmological 
constant, and the number of time steps like t — T/a, keeping T finite (a continuous time interval), we 



Following standard notation, we reserve the letter Z for the partition function, and use the letter G for the propa- 
gator. 

'^^We could of course also obtain the final result by simply counting the number of triangulations of the strip with 
open boundary conditions and obtain 

G{x,y;t = l\X)= Yl x^^y'^e-'('^+^^^ + A = ^— ^ , (1.40) 

which coincides with (1.41). The same is true for the other calculations presented in this section, since most of the 
results in (1-1-1) dimensions can be derived in several different ways. We concentrate on the inversion method because 
it is the one that we will be able to generalize to a particular model of (2-)-l)-dimensional CDT in Chapter 3. 
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get a finite continuous volume^^: 



(F) = lim a2(iV2) ~ ^ 

a^O a/A 



(1.44) 



This expression should be compared with the Euchdean one (V) e ^ T/ K^l^, where the geodesic time 
has to scale anomalously like T = a^/'^t, signaling a fractal dimension (see [20] for details). 

We can see that this infinite-volume limit is sufficient to obtain a continuum theory by evaluating 
the average and correlation for the length of the spatial slices (the only degree of freedom in two 
dimensions). These can be derived with the same inversion method by assigning some weight Zi to 
the dimcr at time ti, that is, by introducing the generating function for triangulations with a slice of 
length l{ti) at time ti, and which we denote by Ft^t^ [u, zi). 

By the inversion formula we have 



;(ti) 



JV2 



from which we find 



Tr 



im)) 



1 iV^^ 

.i^/u J 



'1 
.0 ^1, 



lim a 

a->0 



din Ft^ti{u,z) 



1 iV^' 

^ly/u y 
1 



t-tl 



(1.45) 



9zi 

which does not depend on Ti because of the periodic boundary conditions. 
Similarly, from 

F ill y y "S^ _}_i{ti) lit2) 



(1.46) 



JV2 



Tr 



'1 
.0 ^1, 



' 1 iy^' 

Ay/U J 



'I 0' 

,0 Z2, 



1 iV^' 

iy/u J 



t-(t2-tl)' 



(1.47) 



where we have used the cyclic property of the trace to move the element M *i to the end of the trace 
to make the dependence on t2 — h explicit, we find 

1 d'^Ft^ti,t2{u,zi,Z2)_ 

|zi=Z2 = 



{L{Ti)L{T2)) = lim 



Ff,ti,t2(?i, zi, Z2) dzidz2 
Putting things together, we find for the correlation function of two lengths 



^ I (^1 + e-2(^^-^^)^) . (1.48) 



{L{T,)L{T2))-{m)){L{T2)) 



-2(T2-Ti)VA 



which is the continuum limit of the corresponding correlator, 



{i{t,)i{t2))-{m){i{t2)) 



A+A_e 



-(t2-ti)ln^ 



(A+-A_)2 



(1.49) 



(1.50) 



^^This expression is valid for large T. To be precise, we should take into account also the contribution from A_ («) in 
(1.34), leading to 



{V) 



T 

7x 



-2tVa 



1 - 



1 



-2tVA 



(1.43) 



Note however that for any value of T the correction gives a multiplicative factor between 1 and 1/2 which does not 
change the qualitative result (1-44). However, the correction given by X-{u) can become fundamental when looking at 
other quantities, as is shown by the computation of the correlation function following in the text. 
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from where we identify the lattice correlation length as 

C = (ln^)-i. (1.51) 

This diverges at = 1/4, where \+{uc) = A_(uc), demonstrating the existence of a good continmmr 
limit, and leading to a finite continuum correlation length equal to ^ = lima^o = (2-\/A)^^. 

The existence of correlations between successive slices suggests the presence of a Hamiltonian 
operator for the length of the slices. At the same time, from the result (1.49) we see that the variance 
of the length at T2 — Ti, given by 

^ Vim?) - im))' - ^ , (1-52) 

is of the same order as the expectation value (1.46) of the length itself, which implies that the quantum 
fluctuations will dominate and mask any "semiclassical" behaviour (this is in accordance with the fact 
that we do not expect to see any classical behaviour since in two dimensions there is no classical 
theory of gravity). This phenomenon is nicely illustrated by Fig. 1.8. 



Fi gure 1.81 A typical (l+l)-dimcnsional spacetime as obtained from numerical simulations of the CDT model, at 
volume N2 = 18816 and total time t = 168 [4]. 



1.3.4 Transfer matrix, propagator and Hamiltonian in (d+l)-dimensions 

Because of the global proper time we are dealing with a standard evolution framework. In any 
dimension, the discrete CDT model can be formulated in terms of a transfer matrix (see the following 
section 1.4). This is obtained constructing a vector space W out of the geometrical information 
characterizing the constant-time slices, i.e. we define a 1-to-l correspondence between geometries of 
the slices and state- vectors \g) which we use as basis of W^"^. We can then define the transfer matrix as 

^^Remember that at the discrete level, and for finite volume, the number of different geometries is finite and therefore 
we get a finite-dimensional vector space Wv- If we include the sum over volume, we just get W = Wv, which is 
(countably) infinite. 
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that matrix whose elements in such a basis are given by the one-step propagator with fixed boundary 
data^^, 

, At=l ^ ' 



The transfer matrix can be iterated to obtain the f-step propagator 



G{g„gjm ^ {g,\f%} = ^ cJff'''^^^^ ' ^^"^"^^ 
eventually with periodic boundary conditions in time, giving the partition function 

Zt(A)=Trf*. (1.55) 
If T is a bounded, symmetric and positive operator, we can define a lattice Hamiltonian by 

Hiat = -\lnf. (1.56) 

That the definition (1.53) does actually lead to a bounded and symmetric operator is easy to prove 
[47], while in general (for any dimension) the positivity has only been shown for the case where the 
one-step propagator is replaced by a two-step propagator [47], i.e. it has only been shown explicitly 
that the square of the transfer matrix is positive. However, it is then sufficient to modify the definition 
(1.56) into 

Hit = -Ya^^TK (1.57) 

In general Hiat (or H'l^^) is a very complicated operator which depends on the details of the dis- 
cretization, but luckily wc are not really interested in most of them. Rather, wo arc interested in its 
continuum limit (for which it also should not matter whether we started from Hiat or from H'l^^) 

H = lim Hiat , (1.58) 

a— ^0 

to which we refer as the quantum Hamiltonian of the system (see also [52]). 



1.3.5 Back to (l+l)-dimensions 

In the continuum formulation, pure quantum gravity in (l+l)-dimensions can be thought of as a 
string theory with a zero-dimensional target space. In general, coupling to d scalar fields corresponds 
to having a d-dimcnsional target space. In this case, many results are available from the continuum, 
starting with Polyakov's seminal work [16]. Whenever a comparison is possible, it shows that two- 
dimensional Euclidean DT in the continuum limit agrees with the continuum results of non-critical 
{d ^ 26) string theory in conformal gauge g^^ = e'^S^u (which, as shown by Polyakov, is a quantum 
Liouville theory for the Liouville field (j)). 

What about the continuum limit of (l-l-l)-dimensional CDT? It turns out that also CDT repro- 
duces the continuum results, although in a non-standard gauge. The presence of a time slicing, with 
fixed distance between the slices, suggests comparison with results obtained in the proper-time gauge 

,';o() = 1 ■ ,901 = 510 = , 511 = 7(a;°,a;^) . (1.59) 

-•^^Notc that for lack of letters or fantasy I use the same letter for time, triangulation and transfer matrix, but with 
some variation; t is the discrete time, T the continuous time, T indicates a triangulation, T the transfer matrix and Tjj 
its elements in a given basis. I hope this does not lead to any confusion. 
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Two-dimensional quantum gravity in this gauge was studied in [53] and the results agree with the 
continuum results from CDT (see [42] and [49])^^. The main lesson I want to take from [53] is that 
two-dimensional quantum gravity in proper-time gauge amounts to a problem of quantum mechanics 
rather than quantum field theory. The only degree of freedom surviving such a reduction is the length 
of the slices at constant proper-time, i.e. 

L{x°)= / dxV7(a;°,a;i) , (1.60) 

and the Hamiltonian depends only on this variable (and its conjugate momentum). 

Clearly, CDT fit well into this picture, since having all edges of the same length corresponds to 
having no a; ^-dependence in gu. 

It also brings about a simplification, which allows us to handle the combinatorial problem present 
in the discrete approach. Having to keep track only of a finite number of (discrete) degrees of freedom 
(just one in the specific case), we can rely on the method of generating functions. It is well known in 
mathematics (and physics) that if we want to study a sequence {a„}, it is usually easier to investigate 
its generating function f{x) = J2n C'nX" instead. 

The quantities we want to study in quantum gravity are the propagators 

G(5^.(l),5M^(2),r)= / I?[ff]e«^-«Isl , (1.61) 

which in our discretization arc given by the t-stcp propagators (1.54). In (1+1) dimensions, the states 
\gi) are just labelled by the length, which in lattice units is just a natural number and which we will 
denote by I as before. We thus have 

G{h,l2;t\X) = {h\f%) (1.62) 

for the propagator, and 

G{x,y;t\X)= ^ x'^ y'^ G{hM\t) (1.63) 
ii,;2>o 

for its generating function. Once the generating function is known, we can recover the propagator by 

dx f dy G{x,y;t) 



^ 27ri 27ri x^i+^yh+i ' 

where in the last step we have used Cauchy's formula, and where the contour of integration 7 goes 
around the origin without encircling any singularity of G{x, y; t) (which is analytical in a sufficiently 
small neighbourhood of the origin). 

Propagator and generating function are translated into continuum functions by use of the canonical 
scaling relations 

L T 1 

/ = -, t = — , A = Ac + -Ko? (1.65) 
a a 2 

for the geometric variables and the (dimensionless) cosmological constant, by introducing the boundary 
cosmological constants through^ ^ 

x = e-"^ , y = e-"^ , (1.66) 



l^There are some interesting subtleties, wliicli are probably relevant for a better understanding of the relation between 
Euclidean and Causal DT in two dimensions, but which I will not discuss since it is not the topic of this thesis. 

^'^Note that in (1+1) dimensions the boundary cosmological constants do not need to be renormalized since the 
boundaries have no entropy. 
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and by employing a multiplicative renormalization for the propagator and generating function^* ac- 
cording to 

G(X,Y;T\A)= Urn aG(x,y;t\X) , G(Li, La; TIA) = lim -G(h,l2;t\X) . (1.67) 

o— >0 a— >0 a 

The relations (1.63) and (1.64) become in this limit that of a Laplace and inverse Laplace transform 
respectively: 

/■OO noo 

G{X,Y-T\K)= I dLi dL2e-^^'e-^^^G{Li,L2;T\A) (1.68) 
Jo Jo 

^ ^e^^^e^^=G(X,y;r|A) . (1.69) 

We have already seen how to compute the generating function with the inversion formula, namely, 

1 



G{x,y;t\u = e-^^) 



1 iv^\7i\ 

[ix (1.70) 



A+-A_ 



A,+i - ^/uAt{x + y)- (At+i - At)xy 
where A± are as in (1.35) and A* = — A^. From this we get^^ 



G{X,Y;T\A) = limJ-^G{x,y;t\X) 

a— >0 2 



x/A (1-71) 



Asinh(\/AT) + \/Acosh(\/AT)(X + ¥)+ smh{y/XT)XY 
and from its inverse Laplace transform 

^ ' ^ sinh(\/AT) \smh{VAT)J 



(1.72) 



where Io{z) is a modified Bessel function of the first kind. This expression agrees perfectly with the 
continuum calculation (equation (44) of [53] with m = —1/2, where m + 1 is some kind of winding 
number that in our case is one-half, because we have considered open rather than cylindrical topology 
- see [49] for an interpretation of such winding number in CDT) . 

With the same method we can also compute the Hamiltonian. According to the definition (1.58), all 
we need is the one-step propagator, whose generating function was given in (1.41). Prom T = e~"^'"* 
we see that in the continuum limit we have 

/>oo />oo 

G{X,Y;T = a\A) = dLi dL2e-^^'e-^^^Li\f\L2} = 

.OO ''L , , (1-73) 

= J dLiJ dL2e-^'^'e-^^^Li\(l-aH + 0{a^)j\L2) , 

and hence from the generating function of the transfer matrix can get in the continuum limit the 
Laplace transform of the Hamiltonian (in the "L-representation"). We find from (1.41) 

n 1 /2A — — Y^\ 

G{X, Y;T = a| A) = lim -G{x, y;t = 1\X) = - a (^^^^-^ j + 0{a') , (1.74) 



^*The power of a is fixed by the requirement of obtaining (1.68) and (1.69) from (1.63) and (1.64), given (1.65). If 
this would not lead to a finite propagator, it would mean that (1.65) has to be changed, signaling an anomalous scaling. 
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The factor 2 * is necessary because in (1.74) we have to divide by two in order to get the identity. 
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which upon inverse Laplace transformation yields 

(Li|f IL2) = S{Li -L2)-a (-^2^ - ^ + AL2 ) S{Li - L^) + 0{a^) . (1.75) 



H=-L2^- — + ^L2, (1.76) 



We recognize the Hamiltonian 

_^ d_ 

"dl?2 dL2 

which is self-adjoint (with respect to the measure dL) and bounded from below, as it should be. Its 
eigenfunctions and eigenvalues can be found without too much pain, and are respectively 



ML) = <f4A e-^^Ln{2VAL) , (1.77) 

£:„ = VA(2n + 1) , (1.78) 

where n is a non-negative integer, Ln(x) is the n'th Laguerre polynomial and the eigenfunctions tpn{L) 
are orthonormal. They can be used for a consistency check, showing that 

00 

G{LuL2;T\A) = {L^\e-^^\L2) = e-^^"C(ii)V'n(i^2) (1-79) 

indeed reproduces the result (1-72). 

The model is thus completely solved. From the combinatorics of the discrete triangulation we have 
recovered in the continuum limit the quantum-mechanical proper-time formulation of two-dimensional 
gravity. The Hamiltonian has been found and diagonalized. 



1.3.6 Coupling CDT to matter 

Although highly non-trivial, pure quantum gravity may be considered a rather academic topic, in 
the sense that it is a theory of empty space. In order to make contact with the real world, it is 
important to investigate the properties of the coupled system of gravity and matter. In addition, in a 
regime where matter and geometry interact strongly, the presence of matter can in principle alter the 
results of the pure gravity case. For example, in the context of the perturbative approach, we know 
that the one-loop finiteness result for the pure gravity case is spoiled by the coupling to matter [7]. 
In a non-perturbative context, like that of CDT, the addition of matter degrees of freedom on the 
triangulations can in principle affect the short-distance behaviour. 

Besides understanding how matter affects the geometrical properties of spacetime, it is also impor- 
tant to understand how in turn the behaviour of matter is affected by the fluctuations of spacetime. 
Firstly, this will provide additional cross-checks for a correct classical limit of the gravitational dy- 
namics, and secondly, it will be relevant for potential observable eflfects in nature. 

Matter coupling in quantum Regge calculus and DT has been investigated mostly by numerical 
simulations (exceptions are the already mentioned matrix model methods in two-dimensional DT). I 
will not review here the variety of results that has been found, and will move straight to CDT; the 
interested reader is suggested to start from the review [19]. 

Although a two-dimensional spacetime would not be made more physical by a coupling to matter, 
it is a natural starting point and testing ground for analytical and numerical methods, as usual with 
the hope of extending them to higher dimensions at a later stage. 

Matter coupling in two-dimensional CDT models has been introduced and investigated in [54] and 
[55]. As a prime example of a simple matter model, the Ising model was chosen, probably the most 
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thoroughly investigated matter system on fixed lattices. Both on a regular lattice and on the DT 
ensemble this model has been solved exactly. Unfortunately no exact solution has yet been found for 
the Ising model on CDT. 

In [54] the Ising model on two-dimensional CDT was investigated via Monte-Carlo simulations 
and a high-temperature expansion. These methods provided strong evidence that the model belongs 
to the same universality class as the Ising model on a fixed square lattice, and that the Hausdorff 
dimension of the geometries is dn — 2, and therefore unaffected by the presence of matter. In [55] 
these investigations were extended by coupling eight copies of Ising models to CDT. It was observed 
in Monte-Carlo simulations that even though the matter exponents, and hence the universality class, 
remain unaltered, the geometry has undergone a phase transition, to a phase in which dn = 3. 

In Chapter 2, I will present some new findings that strengthen the results of [54]. The method 
used will be that of high- and low-temperature expansions, whose application to the CDT model will 
be developed in detail. In comparison with the preliminary results of [54], the method will also be 
made more algorithmic. 

1.4 (d+l)-dimensional CDT 

In higher dimensions, CDT arc constructed according to the same principle as in two dimensions, only 
that now the spatial manifold E is higher-dimensional and the triangles are replaced by tetrahedra, 
four-simplices and so on [47] . 

At integer times we have different triangulations of S, constructed of d-dimensional com- 

pletely space-like (Euclidean) simplices ad- The spacetime between two such triangulations at times 
t and t -|- 1 is interpolated by (d -I- l)-dimensional simplices ad+i{d — k, k) having a Ud-k and a Uk as 
subsimplices at time t and i + 1, for A; = 0, 1, d. We distinguish d~\-l types of {d + l)-dimcnsional 
simplices in CDT, to which we will also refer as "of type (d — k, fc)". The three types appearing in 
(2+1) dimensions are shown in Fig. 1.9. 




Fi gure 1.9; The three types of tetrahedral building blocks and their intersections with the ^ + 1/2 plane. 

1.4.1 The combinatorics of the (2+l)-dimensional transfer matrix 

We have seen that in two dimensions it is possible to solve the combinatorial problem for the t-step 
propagator directly. In three dimensions the same task does not seem feasible at the moment. We 
have mentioned that in the Euclidean case this has been tried by a generalization of matrix models, 
but up to now without success. In the CDT case we know that such a matrix model- like description 
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is unlikely to work because of the anisotropy of the triangulations. Furthermore, none of the methods 
used to solve (l+l)-dimensional CDT seem to be generalizable in an obvious way. 

On the other hand, we have also seen how in CDT - once we have obtained a continuum Hamil- 
tonian from the transfer matrix - it is possible to construct the finite-time propagator directly in the 
continuum. This approach seems to have more chances of being soluble because the transfer matrix 
is encoded in the triangulation of the one-step propagator. With a single layer of tetrahedra confined 
between two adjacent slices of constant time, the combinatorial problem is simplified, but still by no 
means simple! 

The nice thing about the one-step propagator of (2+l)-dimensional CDT is that it can be formu- 
lated as a two-dimensional problem "with colours" . In order to do this, we introduce a space-like slice 
at constant half-integer time t + 1/2, i.e. a slice that cuts the (73(3 — k,k) tetrahedra in the middle 
as shown in Fig. 1.9. This slice is by construction a piecewise linear manifold too, but in general 
not a simplicial one. It will consist of squares (coming from the intersection of the time-like faces of 
tetrahedra of type (2,2) with the plane a,tt+^) and of triangles (coming from tetrahedra of type (3,1) 
and (1,3)). We can give a colour-code to the links to distinguish whether they come from triangles of 
type (2,1) (blue links) or of type (1,2) (red Hnks), as in Fig. 1.9. 

Again it is useful to switch to the dual picture. This is done by drawing dual links which connect 
the centre of a triangle or square with the centre of its first neighbours and colouring them with the 
same colour as the link they intersect (see Fig. 1.10). The dual of a blue (red) triangle is a blue (red) 
3-vertex, and the dual of a square is a bi-coloured 4-vertex. 




Figure 1.10: The t + l/2 figures and their dual. 

The dual of the intersection pattern is a graph with blue and red 3-vertices, bi-coloured 4- vertices 
and no bi-coloured 2- vertices. This can also be thought of as the superposition of two three- valent 
graphs of different colours. The blue and red graphs are respectively the dual of the triangulation at 
time t and at time t + 1. 

The next question is which part of the information contained in these two graphs we should keep 
track of and which part we need to sum over in order to obtain the transfer matrix. In principle the 
in- and out-states are labelled by the whole combinatorial information of the blue and red graphs, i.e. 
by the full geometric details of the in- and out-triangulation. If it was really the case that we needed 
to keep track of this complete boundary information, we probably would not have gained much with 
our construction of coloured graphs. We would have to store the information about the two graphs, 
for example, in the form of adjacency matrices, and count in how many inequivalent ways they can be 



30 



CHAPTER 1. INTRODUCTION TO CAUSAL DYNAMICAL TRI ANGULATIONS 



superimposed. It would certainly be a challenging task to relate this complex information to analytic 
data and a quantum Hamiltonian. 

In this situation, the special properties of three-dimensional continuum gravity come to our help. 
Wc know from the canonical analysis of (2+l)-dimensional general relativity that the constraints 
reduce the number of degrees of freedom from infinite to finite (sec, for example, [12]). The remaining 
degrees of freedom correspond to the Teichmiiller parameters of the spatial slices. How this is reflected 
in the quantum theory is less clear. The classical canonical constraint analysis can be implemented 
formally in a phase space path integral, by integrating over the lapse N and the shift Ni variables 
of the ADM decomposition. However, the starting point of the non-perturbative CDT path integral 
is quite different. Firstly, it is the non-perturbative implementation of a configuration space path 
integral, and secondly, it contains by construction no gauge degrees of freedom relating to coordinate 
transformations, and works entirely in terms of geometries, that is, on the quotient space of metrics 
modulo diffeomorphisms. 

How the degree-of-freedom counting of classical, three-dimensional gravity is implemented in a 
configuration space path integral (never mind whether discrete or continuous) is not immediately 
obvious. For example, in a non-perturbative geometric formulation like CDT, it would be interesting 
to see explicitly how most of the degrees of freedom labelling the states of the theory (roughly speaking, 
spatial geometries at a given time) are not propagating, and we are therefore dealing with a theory 
of quantum mechanics rather than quantum field theory. In a formal continuum treatment of the 
path integral, it is clear that this property must be encoded in the non-trivial path integral measure 
obtained after gauge-fixing. A continuum treatment in proper-time gauge {N = 1, Ni = 0), which 
presumably comes closest to the coordinate-free CDT formulation, does indeed strongly suggest that 
the kinetic term for the conformal factor of the spatial geometries (the would-be propagating field 
degree of freedom) is cancelled in the path integral by a Faddeev-Popov determinant in the measure, 
leaving over only a finite number of degrees of freedom [18]. Also the numerical simulations of CDT 
in three [63] and four [65] dimensions show that in the phase where an extended geometry is found, 
there is no dynamical trace of the divergent conformal mode. 

On the basis of these arguments we go back to the transfer matrix of CDT and integrate over 
all but a finite number of variables for the in- and out-state. In particular if we choose a spherical 
topology there are no Teichmiiller parameters, and we will keep track of only one piece of information, 
the area of the slices. We define a new set of states 



where J^iA) is the number of triangulations of given area A (equal to the number of triangles 



In a general theory, this would be just a small part of the full transfer matrix. Only for the special 
case of three-dimensional gravity, on the basis of continuum arguments, we may conjecture that 




(1.80) 



on the slice) and any triangulation with a given total number of triangles A, with the aim of 
computing the transfer matrix in this basis. 



Taa' = {A\f\A'). (1.81) 
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would hold in the large-area limit and we may use (j4|T|A') as our transfer matrix. 
In that case we can again define a generating function for the transfer matrix 

G{x,y;t=l\X,K)= x^'y^^A,\f\A2) , (1.84) 

Ai,A2 

where k is the bare dimensionless Newton constant, having the interpretation of a generating function 
for the number of bi-coloured graphs with a given number of 3-vertices and arbitrary number of 4- 
vertices. Like in two-dimensional DT these graphs can be thought of as Feynman diagrams of a matrix 
model. We can write 

G{x,y;t = l\X,K)= lim In / dAdSe-^^^^'+^'-t^'-f^'-f ^s^^) , (1.85) 
iV— >oo TV J 

where u = xe''^'^~^^^, v = ye'^^'^~^^^ and w = e~'=2«-62A ^j.g ^j^g weights associated to N^i, N13 and 
N22, the number of (3, 1)-, (1, 3)- and (2, 2)-blocks respectively, which come from an evaluation of the 
action for the one-step propagator similar to the one reported in appendix A-I for the case studied in 
Chapter 3. 

In a series of papers, Ambj0rn et al. [57, 58, 59] studied a modified version of CDT with pyramids 
replacing the (3,1)- and (1, 3)-tetrahedra, corresponding to a modified, exactly solvable version of 
the matrix model (1.85) with quartic rather than cubic interactions [60, 61], and analyzed its conse- 
quences. The matrix model generates also many configurations which do not have an interpretation as 
regular slices of a three-dimensional simplicial manifold, but instead correspond to configurations with 
"wormholes" . According to the solution of the matrix model there is a critical line (at u = v, otherwise 
it is a critical surface with a richer phase structure) on which two different phases are realized, which 
in [57] were interpreted in terms of a dominance or otherwise of these wormhole geometries. Taking 
the continuum limit of the CDT model would then correspond to approaching that part of the critical 
line where wormhole geometries are rare. 

One main aim of these investigations was to extract, via the analogue of (1.73), the Hamiltonian 
of the system. Unfortunately it turns out to be highly non-trivial to even identify the leading-order 
behaviour giving rise to the identity part of the transfer matrix [59], and up to now a derivation of 
the Hamiltonian has not been possible. 

In Chapter 3, I will present a detailed analysis of a different CDT model for three-dimensional 
gravity for which some more progress can be made in this direction. 

1.4.2 Numerical results 

A valuable tool in CDT is represented by Monte Carlo simulations. Since exactly solvable models 
are rare, numerical simulations play a vital role in physics. Being able to study quantum gravity via 
simulations presents a great advance, which goes back to Regge calculus and DT. Numerical results for 
CDT for the first time produced tangible evidence for the emergence of extended and "semiclassical" 
geometries in higher-dimensional non-perturbative gravity. 

Monte Carlo simulations of (2-|-l)-dimensional CDT with topolog y S^x [0, 1] [62, 63] showed the 
existence of a single phase, where the typical geometry represents a "lump" of spacctime with a three- 
dimensional scaling behaviour. In Fig. 1.11 a plot of the area of the slices as a function of time is 
reproduced, to illustrate the presence of such an extended lump. 

From an analysis of the fluctuations of the area of successive spatial slices, evidence was gathered 
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Figure 1.11: The "lump" of 3d geometry as produced by Monte-Caxlo simulations of three-dimensional CDT [63]. N2 
(vertical axis) is the number of triangles in the slices of constant time t, i.e. the area of the slices. The plot is reflected 
about the central t-axis to give a pictorial impression of the emerging spacetime. 

for an effective description by an action of the form 



which is exactly the kind of action we would find following, for example, [56]. Numerical results 
support the existence of a well-defined dynamics in terms of a reduced set of degrees of freedom (just 
a single one for the sphere). Whether this has to be thought of as an effective mini-superspace-like 
reduction or as a reduction coming from the constraints of the full theory should be determined in 
some other way, for example, by checking the property (1.82). 

A similar result was obtained more recently for (3-|-l)-dimensional CDT [64], where we expect to 
have all the field degrees of freedom of general relativity. There seem to be three different phases, 
separated by lines of first-order transition. Evidence for the emergence of an extended four-dimensional 
universe has been found in one of these phases, the others being again dominated by degenerate 
geometries. 

By an analysis similar to that in (2-|-l) dimensions, an effective action for the volume of spatial 
slices (within the extended phase) was found to be given by [65] 



Introducing a scale factor a{t) by V{t) = a^(t), one finds a simple mini-superspace action often used 
in quantum cosmology, apart from an overall sign. The fact that such sign differs in the effective 
action obtained in CDT (leading to an action bounded from below) from that obtained by simply 
substituting a mini-superspace ansatz in the action (resulting in an unbounded action) supports the 
idea that a correct inclusion of the contribution of the measure in the path-integral would also cure 
the conformal factor problem in a continuum path integral treatment. 

Finally, 1 would like to mention another very interesting result found in (3+1) dimensions from the 
numerical simulations. As already mentioned, various notions of effective dimension can be measured 




(1.86) 




(1.87) 
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within this approach. By a global scaling of the volume-volume correlator it was found that in the 
extended phase the value of the Hausdorff dimension is compatible^" with dn = 4: [66]. Another 
possible notion of effective dimension is the so-called spectral dimension ds, which I will not describe 
in dcrtail here. While also for this dimension a value compatible with = 4 on large scales was found, 
the value of the spectral dimension turns out to be scale-dependent and decreases monotonically 
from = 4 on large scales to ds ~ 2 on short scales [67]. A similar result has subsequently been 
found within a non-perturbative renormalization group approach [68] . Such a dynamical dimensional 
reduction is very appealing in that it indicates a way in which gravity may cure its own ultraviolet 
divergences by virtue of its asymptotic safety. It is conceivable, at least in principle, to construct a 
non-perturbatively improved pcrturbative expansion where due to the dimensional reduction at short 
distances there would be no such divergences. 

1.5 Summary 

In this chapter I have given an introduction to the CDT approach to quantum gravity, together 
with its motivations, background, and results. Particular attention has been given to what is known 
analytically, since the purpose of the work presented in the following chapters is to elaborate on it 
and to try to push the analytical understanding of the CDT model further. 

We have seen that starting from the idea of constructing a non-perturbative quantum theory of 
gravity in analogy with the lattice formulation of quantum field theory, Regge calculus provides us 
with a natural discretization of general relativity. We reviewed the arguments that led from quantum 
Regge calculus, where the triangulation is fixed and the edge-lengths are the dynamical degrees of 
freedom in the path integral, to dynamical triangulations, where the edge-lengths are fixed and the 
combinatorial structure of the triangulation is being summed over. 

The original Euclidean dynamical triangulation approach has a number of advantages over the 
quantum Regge calciilus, but ultimately leads only to macroscopically pathological geometries, the 
main reason for which it has mostly been abandoned. As I have tried to highlight, the causal version 
of DT offers a physically motivated cure to its troubles. By contrast, CDT models have succeeded in 
recovering the large-scale effective dimension d = n when starting from n dimensional building blocks, 
and have provided strong evidence of the emergence of a semi-classical spacetime from the sum over 
geometries. 

One of the great merits of dynamical triangulations is that they give a definite prescription for 
studying quantum gravity by computer simulations, and by that obtain definite (albeit numerical) 
results. On the other hand, the discrete nature of the model makes it difficult to obtain analj^ical 
results, and the only CDT model so far whose continuum limit has been obtained by purely analj^tical 
methods is the one in (H-1) dimensions without matter coupling. 

I will introduce in Chapter 2 a "semi-analytical" technique that can be used in two dimensions 
to obtain information about the critical behaviour (and hence the continuum limit) in the case of 
gravity-matter coupled systems. 

I have reviewed the solution of the (l-l-l)-dimensional pure- gravity CDT model obtained by use 
of the inversion formula for heaps of pieces. The same technique will be of great value in Chapter 3, 
where it will be applied to the evaluation of the one-step propagator for a special version of (2-1-1)- 
dimensional CDT, in the same spirit as the work on the ABAB matrix model mentioned in the last 
section, demonstrating that a continuum limit can be performed analytically also in that case. 



By numerical simulations it is of course impossible to prove that a certain measured number is an integer. 
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Chapter 2 



Adding matter in two dimensions 



In this chapter I will present work on spin systems coupled to Causal Dynamical Triangulations (CDT), 
more specifically, on the Ising model on a (l+l)-dimensionaJ CDT. It is shown how a method invented to 

analyze the critical behaviour of spin systems on flat lattices can be adapted to the fluctuating ensemble of 
curved spacetimes underlying the CDT approach to quantum gravity. A systematic counting of embedded 
graphs is developed to evaluate the thermodynamic functions of the gravity-matter models in a high- and 
low-temperature expansion. This chapter is based on publications [97] and [98]. 

2.1 Coupling quantum gravity to matter 

If a quantum theory of gravity is to describe properties of the real world, it must tell us if and how 
matter and spacetime interact at extremely high energies. Current, incomplete models of quantum 
gravity usually approach the problem by first trying to construct a quantum theory of spacetime's 
geometric degrees of freedom alone, and then adding matter degrees of freedom in some way. One 
is then interested in whether and how the "pure" gravity theory gets modified and how the dynamic 
aspects of quantum geometry may influence the matter behaviour. A priori, different scenarios at very 
short distance scales are thinkable: the behaviour of quantum geometry may be completely dominant, 
geometric and matter degrees of freedom may become indistinguishable, or the matter- coupled theory 
may have no resemblance with the pure theory at all. 

In the previous chapter we have seen how difficult it is to study the quantum behaviour of pure 
gravity. As one may expect, coupling matter to geometry docs not improve the situation, and under- 
standing how the known behaviour of matter models on flat spacetime is affected by the fluctuations 
of geometry is a challenging problem. Most of the difficulties stem from the fact that in the case of 
quantum gravity, conventional calculational methods have to be adapted to a situation where there is 
no fixed backgroTind spacetime, and instead "geometry" is among the dynamical degrees of freedom. 

In a context whore we try to make sense of quantum gravity by a rcgularization of the path integral 
in terms of discretized geometries, like is the case of CDT, matter fields are replaced by lattice models. 
Typical lattice models are spin systems on a lattice, which have been extensively studied in the case 
of regular (flat) lattices, and for which many techniques exist. In this case, the problem of studying 
matter on a fluctuating geometry is transformed into the problem of studying spin systems on an 
ensemble of random lattices, in particular, the causal triangulaticms. It is then not clear which, if any, 
of the techniques developed for the case of a fixed lattice can be generalized to the case of a fiuctuating 
lattice. 
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In this chapter we will look in detail at an example of such a generalized method, for a particular 
class of quantum systems of gravity and matter. As it should be clear, the gravity part will be played 
by the CDT model, and, as announced in the title of the chapter, we will work in the simplified context 
of two spacetime dimensions, where wc have full control on the pure-gravity model, and where many 
exact results for spin systems on fixed lattices are available. 

We will see how a time-honoured method of estimating the critical behaviour of a lattice spin 
system, namely, by considering finitely many terms of a weak-coupling (i.e. a high-temperature) ex- 
pansion of its thermodynamic functions can be adapted to the case of "quantum-gravitating" lattices. 
The method consists in classifying local spin configurations (which can be represented by diagrams 
consisting of edges or dual edges of a triangulated spacetime), and determining their weight, i.e. the 
probability to find such a configuration in the ensemble of dynamical triangulations. Diagrammatic 
techniques have also been employed recently in another coupled model of geometry and matter in 
three dimensions, in the context of a so-called group field theory, an attempt to generalize matrix 
model methods to describe ensembles of geometries in dimension larger than two (see [41] for a recent 
review) . 

The model wc will be studying in detail is that of an Ising model whose spins live on cither the 
vertices or triangles of a two-dimensional simplicial lattice and interact with their nearest neighbours. 
Summing over all triangulated lattices in the usual CDT ensemble, each with an Ising model on it, 
gives rise to the matter-coupled system whose properties we are trying to explore. An exact solution 
to this Lorentzian model has not yet been found. This has to do with the fact that in terms of the 
randomness of the underlying geometry the model lies in between that on a fixed regular lattice (solved 
by Onsager long time ago [69]) and that on purely Euclidean triangulations (solved by Kazakov [38]), 
but is sufficiently different from either to make the exact solution methods known for these cases 
inapplicable. 

As argued in Chapter 1, it is difficult to introduce additional restrictions concerning the local 

nature of geometry into a matrix model without destroying its simplicity and, more importantly, its 
solubility. As seen there, this does not prevent the pure-gravity model to be exactly solvable by other 
methods, but unfortunately the same cannot be said up to now for the Ising-gravity coupled case. 
Of course there is no obstacle to studying such a system by numerical simulations. This was done in 
[54, 55], where strong evidence was found that the spin model behaves in the presence of gravity (at 
least as far as its critical exponents arc concerned) just like it docs on a fixed, regular lattice. This 
is supported by a "semi-analytical" high-temperature expansion, whose initial results were merely 
quoted in the original work of [54] . In the work I am reporting in this chapter, the details of how to 
perform the relevant graph counting on CDT are developed, both for the standard Ising model on CDT 
and on its dual. The algorithms found have allowed us to compute the diagrammatic contributions 
to the Ising susceptibility to order 6 and order 12 for CDT lattices and their duals, respectively. 

The pure counting results were first announced in [97], where we put forward the conjecture that 
the best way to extract information about the universal properties of a matter system on flat space is 
by coupling it to an ensemble of quantum-fluctuating causal triangulations! The strongest evidence 
for this conjecture so far comes from the fact that an evaluation by straightforward ratio method of a 
rather small number of susceptibility coefficients already leads to results in surprisingly close agreement 
with the exactly known value for the susceptibility. As will become clear below, the evidence for this 
conjecture is still too limited to reach any definite conclusion, and the issue is being studied further. 

Generally speaking, one can look at our method and results also from a purely condensed matter 
point of view, as an example of how the universal properties of spin and matter systems are affected 
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by introducing random elements into its definition. For example, the inclusion of random impurities in 
a lattice model leads to Fisher renormalization [70] of the critical exponents in the so-called annealed 
case, where the disorder forms part of the dynamics. By contrast, in the quenched case the Harris- 
Luck criterion is conjectured to hold [71]. To what extent these ideas extend to the case of random 
connectivity of the lattice which is the one relevant to non-perturbative quantum gravity models - is 
still not very clear, despite a number of results in this direction, involving spin models with (quenched) 
geometric disorder coming from Euclidean dynamical triangulations and Voronoi-Delaunay lattices 
(see, for example, [72, 73]). 

One of the reference points for such investigations is that of the two-dimensional Ising model on 
fixed, regular lattices, which in absence of an external magnetic field can be solved exactly in a variety 
of ways (see [74, 75]). In this case, the critical exponents, which characterize the system's behaviour 
near its critical temperature, are given by the so-called Onsager values a = 0, /? = 0.125 and 7 = 1.75 
for the specific heat, spontaneous magnetization and susceptibility. In sharp contrast, in the case of the 
Ising model on Euclidean dynamical triangulations mentioned earlier the disordering effect is so strong 
that the same critical exponents are altered to a = —1, /? = 0.5 and 7 = 2 [38]. On the other hand, 
despite the annealed geometric randomness present, the corresponding Ising model coupled to causal 
dynamical triangulations seems to share the flat-space exponents, indicating a perhaps surprising 
robustness of the Onsager universality class. ^ 

In the piece of work presented here, we will concentrate for simplicity and definiteness on the 
CDT-Ising system, although our method should also apply straightforwardly to other spin models and 
may inspire similar diagrammatic expansions in other, discrete models of quantum gravity. Sec. 2.2 
introduces the high-temperature expansion on regular and dynamical lattices. In Sec. 2.3, after 
recalling some general definitions and results of graph theory, a set of rules is presented which in 
principle allow us to compute the weight of any susceptibility graph on the fiuctuating CDT lattice 
and its dual, and explain the method with the help of illustrative examples. The finite number of 
terms of the series expansions we obtain in this way are analyzed in Sec. 2.4 with the help of the ratio, 
Dlog Fade, and differential approximants methods. Sec. 2.5 contains some remarks on the extension 
of our method to low-temperature expansions. 



2.2 The Ising model and its high-T expansion 

Given any lattice G of volume N with v vertices and I links (edges), the Ising model on G in the 
presence of an external magnetic field h is defined by the partition function 

Zn{K,H)= 6-^''""^'= E e^^«>^^'^^+^^*'^S (2.1) 

{(T«=±l}«eG {ai=±l}«eG 

where W[{cr}] is the Hamiltonian of the spin system, /3 := l/(fcsT) the inverse temperature, fjj the 

spin variable at vertex i, taking values di ~ ±1, and {ij) denotes nearest neighbours. In standard 
notation, wc will use K = f3J, where J > is the ferromagnetic spin coupling, and H = (imh, where 
m is the spin magnetic moment. The fact that OiGj = ±1 allows us to write 

gK<T,:<T, = n + y o-jcr,-) cosh(ii'), 

^ ' ^ (2.2) 

e""* = {l+Tai)cosh{H), 



The study of fermions coupled to two-dimensional causal dynamical triangulations was initiated in [76] . 
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in terms of newly defined variables u = tanh(iir) and r = tanh(if), and rewrite the partition function 
as 

= cos\i\K) cosW{H) + u ai<jj) + r ai) 

= cosh' (JC) cosh'' (i?) ^ [(1 +M^f7i<Tj -|-w^ ^ (Jiajakai + ...)x (2.3) 

{<T^} (ij) {ij),{kl),j^l 

X (l+T^CTi+T^^CTiCTj + ...)] 
i ijtj 

= T cosh'(ii') cosh" (if) F^^\u, r). 

Because of the sum over spin values ±1 it is clear that every term containing a spin at a given vertex 
to some odd power will give a vanishing contribution. Thus we need to keep only terms with even 
powers of Ci's, which arc terms where each ai belongs cither to an even number of nearest neighbour 
couples OiOj (from the M-terms) or to an odd number of them and to one of the spins coming from 
the r-terms. It is not difficult to convince oneself that each such term is in one-to-one correspondence 
with a graph drawn on the lattice whose length is given by the power of u and whose number of 
vertices with odd valence is given by the power of t. The function F'^^^UjT) introduced in eq. (2.3) 
is then the generating function for the number of such graphs that can be drawn on the lattice, and 
the factor 2" comes from the sum over the spin configurations. The representation we have obtained 
in this manner is a high-temperature expansion, since for infinite temperature T we have u = r = 0. 
We will return to this graphic interpretation in the next section, after having recalled some basic 
definitions from graph theory. In what follows, we will only be interested in the case of vanishing 
external magnetic field. 

All we have said up to now does not require any specific properties of the lattice G, but works 
for any lattice, regular or not. This makes it straightforward in the framework of CDT to couple the 
Ising model to gravity. We simply associate a complete Ising model with each triangulation of the 
ensemble, viewed as a lattice, by putting the spins at the vertices of the triangles, and then perform 
the sum over such triangulations, leading to 

Z{K,H,N) = J2 e-^"^"""- (2.4) 

The situation where the volume is allowed to fluctuate is also of interest in a quantum gravity context 
and described by the grand canonical partition function 

Z{K, H)=J2 e-^^'ZiK, H, N), (2.5) 

N 

where the role of chemical potential is played by the cosmological constant A. As usual in dynamically 
triangulated models, the asymptotic behaviour as function of N of the canonical partition function in 
the infinite- volume limit N oo will determine the critical value of A at which the continuum limit 
can be performed. 

In trying to understand the behaviour of the Ising model coupled to quantum gravity in the form of 
causal dynamical triangulation, we will make crucial use of the known probability distribution P{q) of 
the coordination number q of vertices (the number q of links meeting at a vertex) for the pure-gravity 
case. This is relevant because we will be expanding about the point /3 = 0, at which the spin partition 



2.2. THE ISING MODEL AND ITS HIGH-T EXPANSION 



39 



function (2.3) reduces to a trivial term 2'" and the geometry is therefore that of the pure 2d Lorentzian 
gravity model. The probabihty distribution in the thermodynamic Hmit was derived in [54], resulting 
in 

m = 1^- (2.6) 

The main idea is to use this distribution, together with some information about correlations between 
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Figure 2.1; Comparing the probability distributions P{q) of the vertex coordination numbers q for three different 
types of random triangulations: EucUdean dynamical triangulations [77], Poissonian Delaunay triangulations [78] and 
causal dynamical triangulations [54] . 

distributions at different vertices, to re-express the sum over triangulations in the high-temperature 
expansion of the partition function, (2.4), by an average over coordination numbers. In terms of 
geometric randomness, the pure CDT model lies in between the Voronoi-Delaunay triangulations 
based on Poissonian random distributions of vertices and the planar triangulations underlying the 
approach of Euclidean dynamical triangulations. A comparative plot of the probability distributions 
for the coordination numbers in the three different types of geometry is shown in Fig. 2.1. 

Let us point out a further subtlety which arises from the fact that the underlying lattices are not 
fixed, but fluctuating, and therefore do not have a fixed "shape". In this case it can happen that in 
the graph counting for the high-temperature expansion of some thermodynamic quantity topologically 
non-trivial graphs must be taken into account. By this we mean closed non-contractible graphs which 
wind around the spacetime. This does not invalidate the method in principle, but requires a more 
detailed knowledge of the global properties of the spacetime. In two-dimensional CDT, where the 
spacetime topology is usually chosen to be that of a cylinder (a spatial circle moving in time), this 
would be closed loop graphs which wind one or more times around the spatial , and might be as 
short as a single link. Such pinching configurations certainly exist, but their contribution to the graph 
counting has been shown to be irrelevant in CDT, because it is subleading in N in the thermodynamic 
hmit [54]. 
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2.3 Graph embeddings 

2.3.1 Terminology of graph theory and embeddings 

To prepare the ground for our counting prescription, we will first review some of the relevant termi- 
nology of graph theory and embeddings, along the lines of reference [79] . 

A linear graph is a collection of v vertices and I lines (or edges) connecting pairs of vertices. A 
simple graph is a graph in which two vertices are connected by at most one line and in which single 
lines are not allowed to form closed loops. A graph is said to be connected if there is at least one path 
of lines between any two vertices, and disconnected if for some vertex pair there is no such path. 

The cyclomatic number c of a connected graph g is defined as 

c{g)=l-v + l, (2.7) 

and represents the number of independent cycles in the graph. 

The degree (or valence or coordination number) of a vertex is the number of edges incident on that 
vertex. 

Two graphs are said to be isomorphic if they can be put into one-to-one correspondence in such a 
way that their vertices and edges correspond. They are called homeomorphic if they are isomorphic 
after insertion or suppression of any number of vertices of degree 2 (this operation being defined on 
an edge in a trivial way). 

Homeomorphs are thus graphs with the same topology, in particular, with the same cyclomatic 
number. It is then possible to classify graphs in terms of irreducible graphs. An irreducible graph is 
one in which all vertices of degree 2 have been suppressed. (Note that in general it will no longer be 
a simple graph since it will contain single-line loops.) Typical irreducible graphs whose homeomorphs 
we will encounter are tadpoles, dumbbells, figure eights and 9-graphs (Fig. 2.2). 




(a) (b) (c) (d) 

Figure 2.2: Four kinds of homeomorphs: (a) tadpole, (b) dumbbell, (c) figure eight, and (d) 6-graph. 

In the high-temperature expansion of an Ising model it is useful to think of the lattice with its 
V vertices and / edges as a graph G. The graphs appearing in the expansion of the thermodynamic 
function of interest are then graphs embedded in G. A general definition of embeddings of graphs 
is the following. Let G be a general graph. A graph H is a subgraph of G when all its vertices and 
edges are vertices and edges of G. An embedding^ of a graph g in G is a subgraph G' of G which is 
isomorphic to g. 

The lattice constant {g; G)*^^-* of a graph g on G is defined as the number of different embeddings 
of g in G. In practice it is often useful to work with the lattice constant per site, defined as (<?; G) = 

^To be precise, this defines a weak embedding. By contrast, a strong embedding of g in G is defined as any section 
graph G+ of G which is isomorphic to g, where a section graph is a subgraph of G consisting of a subset A of vertices 
and all the edges which connect pairs of (nearest-neighbour) vertices of A. Since all our calculations will concern weak 
embeddings and the corresponding weak lattice constants, we will from now on drop this specification. 
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With these definitions in hand we can now rewrite the function F'-^\u,t) introduced in (2.3) as 

oo 

FWKr)^l + $:$:A%Vr^«, (2.8) 

where -D; are the lattice constants of graphs of length I and with 2s odd vertices, for which F^^^ {u, r) 
serves as a generating function. It follows from the extensive nature of the free energy of the system 
that in the thermodynamic limit 

fW ^u, t) ""-^ e^®("'") = 1 + NQ{u, r) + 0{N^), (2.9) 

where the function 0(w,r) docs not depend on N. From the usual combinatorial theory of non- 
embedded graphs it is well known that the logarithm of the generating function is the generat- 
ing function for connected graphs. Because of the embedded nature of the graphs in the case at 
hand, the disconnected graphs will still contribute to logF(^)(u, r) with a "repulsion" term of sign 

^ ^^{^connected components) — ! 

In our analysis of the Ising model on CDT we will be looking at a quantum-gravitational average 
of the magnetic susceptibility at zero external field. For the usual Ising model the susceptibility x per 
unit volume is given by 

1 dHnZ^ 1 2 ErOgV 

where in the second step we have substituted in the high-temperature expansion of the partition 
function. By virtue of (2.9) it is clear that the susceptibility is independent of N in the infinite- 
volume limit (thus justifying our notation x(^))) and that moreover the denominator in the last term 
in (2.10) does not contribute at lowest order in N, which is the one relevant to the computation. To 
calculate the susceptibility, it is therefore sufficient to compute the term of order A'' in Di^^u'', as 
is well known. 

The computation in the gravity-coupled case proceeds completely analogously. Performing the 
sum over all triangulations at fixed N, as in eq. (2.4), and then letting N ^ oo, one obtains 



for the susceptibility in presence of the gravitational CDT ensemble, where we have dropped the 
irrelevant constant term and introduced an obvious notation for the ensemble average in the second 
step. Whereas it was fairly straightforward to verify the cancellations of terms of higher order in N 
between numerator and denominator in the pure Ising case in eq. (2.10) by an explicit calculation, 
the analogous computation does not seem feasible in the gravity-coupled case, although we know from 
the same general arguments that it must be realized here too. Again the terms in the denominator of 
the expressions on the right-hand side of (2.11) do not contribute at lowest order, which - at least for 
this particular observable - implies that its annealed and quenched gravitational averages coincide. 

2.3.2 Useful results on lattice constants 



We now want to identify which graph topologies appear in the high-temperature expansion of the 
magnetic susceptibility. For convenience we can divide our set of graphs into two subsets, graphs 
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with zero cyclomatic number (also called Cayley trees), and graphs with c > 0. Since the graphs we 
encounter in the susceptibility series must have two and only two odd vertices, it follows that the 
first subset will only contain connected non-selfintersecting open chains (in other words, self-avoiding 
walks), while the second will contain everything else. The connected graphs in the second subset are 
graphs with loops and one or no open ends, while the disconnected ones take the form of a union of 
one of the previous two kinds with any number of even graphs (graphs in which every vertex has even 
degree). 

The problem we meet in the high-temperature expansion is the evaluation of the lattice constants 
of such graphs. It turns out that there are two kinds of theorems which can be applied profitably in 
the context of our random triangulated lattices. The first is very general and reduces the calculation 
of the lattice constants for disconnected graphs into that for connected graphs. This makes it possible 
to compute the lattice constants for the second subset (c > 0) in a unified way. The second theorem 
is less general but powerful, since it gives us a recurrence relation for the expansion coefBcients of the 
susceptibility. 

The reducfAon theorem for disconnected graphs states that if gi and gj are two graphs gi ^ gj and 
G is any graph, then 

{gi U gj;G) = {gi; G){gj;G) - J^idi + 9i = 9k}{9k\ G), (2.12) 

k 

where gi U gj stands for the disjoint union of the two graphs, {gi + gj = gk} is the number of possible 
choices of cmbeddings of gi and gj in gk having g^ as their sum graph, and the summation is over all 
graphs gk obtainable in this way. If gi = gj we can compute the right-hand side of (2.12) as if the two 
graphs had different colours and thus obtain 2{gi LI gf, G) . U gi or gj are themselves non-connected 
graphs we can iterate the theorem. In this way the lattice constant of the disconnected union of n 
connected graphs can be expressed as a polynomial of order n in the lattice constants of connected 
graphs (see Fig. 2.3 for some simple examples). 



a) 



b) 




Figure 2.3; Two examples of the reduction theorem. The graph G has been dropped in the notation. 

The counting theorem for the susceptibility series was first given by Sykes [80] in 1961, but only 

later proved by Nagle and Temperley [81] as a special case of a more general theorem of graph 
combinatorics. Closer examination reveals that the proof makes no reference to the structure of the 
lattice, but only to its coordination number q, which will turn out to be useful when considering 
the lattice dual to the triangulation. The theorem states that the high-temperature series for the 
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susceptibility can be written as 



X(m) = 1 + (1 - (Ju)-^ [qu{l 



(ju) - 2(1 - u'')Si + 8(1 + uf{S2 + Si)] , 



(2.13) 



where q = a + 1 is the coordination number, and the functions 5, are defined as 



Si 



l>0 



S2 



^ ^wdi{w) 



u 



(2.14) 



l>0 \_w>0 



Ss 




l>0 r<r' 



The coefficient di is the sum of lattice constants of even graphs with I lines, while di{v) is restricted to 
even graphs with I lines and vertices of degree 2gi, 2q2, characterized by the weight w — J2i nin " 
l)/2. The function mi{r,r') is the sum of lattice constants of graphs with exactly two odd vertices 
of degrees 2r + 1 and 2r' + 1. Substituting the left-hand side of formula (2.13) by the expansion 
x{u) = 1 + J2n>o ^nu", one finds a recursive relation for the coeSicients a„, namely. 



As a result, at every new order I in our calculation the only new lattice constants we have to evaluate 
are those for graphs with no vertex of degree 1, reducing the computational effort considerably. 

This formula can in principle be applied to any non-regular lattice whose vertices are all of the same 
degree. For a fixed non-regular lattice this observation would usually not be of much help, because 
one would still need to compute at each order the new lattice constants depending on the complete, 
detailed information of the lattice geometry. In our case this difficulty is not present, since we are 
summing over triangulations and can simply substitute in (2.15) the lattice constants averaged over 
the triangulations, without the need of keeping track of the lattice geometry for individual lattices. 

2.3.3 Counting graphs on the CDT trianguleir lattice 

Next, we will analyze the evaluation procedure for lattice constants in the quantum-gravitational CDT 

modcJ. Computations here are made difficult by the randomness of the coordination number. At the 
outset it is not even obvious that any of the known recurrence relations can be used. Our task is to 
list all possible diagrams, and for each count the number of ways it can be embedded in the lattice. 
On a regular lattice such an operation is tedious but straightforward: starting from a fixed vertex, we 
trace out all possible sequences of links of a given length, say, which do not self- intersect. The discrete 
symmetries of the regular lattice simplify this task greatly. In the simplest case, the counting has to 
be done only for a single initial vertex, with all other graphs obtained subsequently by translational 
symmetry. The analogous counting on the dynamical CDT lattices is complicated by the fact that 
the lattice neighbourhoods of different vertices will in general look different. We will deal with this 



ai - 2aai-i + a^ai-2 = 2{l - 2)c?;_2 - 2ldi+ 




(2.15) 



.w>0 w>0 w>0 
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difficulty by setting up an algorithm to count the embedding of a given diagram on the ensemble of 
CDT lattices, making use of the known probability distribution of the vertex coordination numbers. 

Before presenting this algorithm we need some more notation. CDT lattices have two kinds of 
links, time-like and space-like^. In keeping with the usual representation of two-dimensional CDT 
lattices (c.f. Fig. 1.6, left-hand side), wc will draw time- and space-like links in embeddings of graphs 
as (diagonally) upward-pointing and horizontal lines, respectively (see Fig. 2.4). Furthermore we 
need to keep track of the relative up-down or right- left orientation of consecutive links, to distinguish 
between graphs like (c) and (d) in Fig. 2.4 (which give different contributions to the lattice constant). 
In addition, the counting for a graph like (b) will be identical to that for its mirror images under left- 
right and up-down reflections. These mirrored graphs will be counted by multiplying the embedding 
constant of the graph by an appropriate symmetry factor. A symmetry factor 1 is assigned to graphs 




(a) (b) (c) (d) 

Figure 2.4: The possible typologies of embedding, up to reflection symmetry, of a length-2 open chain. The angle 
which the time-like links form with the vertical has no relevance for the graph counting. 

symmetric under both up-down and left-right reflections (for example, graphs (a) and (d)), a factor 2 to 
graphs with only one of the two symmetries (like graph (c)) or symmetric with respect to a composition 
of the two (like graph (c) of Fig. 2.5), and a factor 4 to graphs without reflection symmetry (like graph 
(b) of Fig. 2.4 or (a) and (6) of Fig. 2.5). 




(a) (b) (c) 

Figure 2.5: Three of the different typologies of possible embeddings of a length-3 open chain. 

The lattice constant (51; T) of a graph g embedded in a graph corresponding to a triangulation T 
is computed as 

{g;T) = J2sii9;T)i, (2.16) 

i 

where (g; T)i is the number of ways a particular typology of embedding {g)i can be realized on T, Si 
is the relevant symmetry factor, and the sum extends over all possible typologies of embeddings of 
the graph g. In the context of quantum gravity, we are interested in the average of (2.16) over all 



^Even after performing the Wick rotation to Euclidean signature, these two link types remain combinatorially 
distinguishable because of the special way in which the original Lorentzian simplicial geometries were constructed. 
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triangulations, that is, 

(g) = (2.17) 

which will take into account the probability distribution of the vertex coordination. 

As mentioned earlier, all calculations will be performed in the thermodynamic limit where the 
cosmological constant is tuned to its critical value Ac = In 2. In this limit, the probability of having k 
incoming time-like links at a particular lattice vertex is given by 

P{k) = ^, (2.18) 

with an identical probability for having k outgoing time-like links at a vertex. This probability 
distribution can be derived easily in the following way. At a vortex there is always at least one 
outgoing link, which can be identified with the right vertical edge of the upward-pointing triangle to 
the left of the vertex (see Fig. 2.6). To the right of such edge there will be always another triangle, 
which in the critical limit will be with probability 1/2 either upward- or downward-pointing. If it is 
upward-pointing, the edge in common to the two triangles is the only vertical link coming out of the 
vertex in question, such a situation having probability 1/2. If it is a downward-pointing one, we can 
continue to its right and repeat the observation until a upward-pointing triangle is found at the k-ih 
step, giving a vertex with k vertical outgoing links, such a situation having probability 1/2^^. 




Figure 2.6: A strip in a triangulation of the CDT ensemble, and, in shadow, the triangles contributing to the weight 
at the vertx i. 

Moreover, at one and the same vertex the two probabilities (incoming and outgoing) are indepen- 
dent of each other. In fact, not just for a single vertex are these two probabilities independent, but 

the same is true for all vertices lying in the same space-like, horizontal slice. By contrast, the outgoing 
probability at a vertex will in general condition the incoming probability at another vertex on the 
subsequent horizontal slice. By construction, the probability of having a space-like link to the left and 
to the right of a given vertex is equal to 1. 

Armed with this information we are now ready to compute any of the {g; G)i in (2.16). We do not 
have a general coimting formula in closc^d form, but will formulate a number of rules and tools which 
will enable us to do the counting recursively. We will start with some illustrative examples. Examples 
3 and 4 will serve as our elementary building blocks in more complicated constructions. 

Example 1. The simplest example is that of a length- 1 graph ci, consisting of a single horizontal or 
vertical link. The number of horizontal embeddings is the number of horizontal links on a CDT lattice, 
which by virtue of the Euler relation is given by half the number of triangles, i.e. 7V/2. If we give the 
lattice constant per vertex, as is customary in the literature, it is therefore {ci)h, = 1- By a similar 
topological argument we can immediately deduce the number 2 for the lattice constant of vertical 
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embeddings, but it is instructive to compute it from the probability distribution (2.18) instead. If at a 
vertex we have k outgoing vertical links, there are precisely k ways to embed the one-link graph such 
that it emanates in upward direction from this vertex. Since every vertical embedding is accounted 
for in this way, by assigning the relevant probability and summing over k we obtain 

oo , 

fe=i 

in agreement with the earlier argument. The total contribution to the lattice constant at order 1 is 
therefore {ci)h + {ci)v = 3. 

Example 2. Next, let C2 be a length-2 open chain, and consider the embedding typology C2,6 shown 
in Fig. 2.4(6). The result follows straightforwardly from example 1: there is one horizontal link per 
vertex, and there are in the ensemble on average two ways of attaching an outgoing vertical line to its 
right vertex, yielding {c2)b = 2. For the embedding of Fig. 2.4(a) there is nothing to compute; we have 
(c2)a = 1- The embedding of Fig. 2A{d) is the first configuration we encounter which extends over 
two triangulated strips. Because of the independence of the probabilities (2.18) in successive strips 
we obtain 

oo , 

(^2)d = (E^)'=4- (2-20) 

Example 3. Now consider the embedding C2,c depicted in Fig. 2.4(c). If the vertex on top has k 
incoming links, there are exactly fc(fc — l)/2 ways to realize the embedding, leading to 

(^^)c = f:^|^ = 2. (2.21) 

Note that the same "inverted- V" embedding on a flat triangular lattice T^eg would give (c2; T^e^c = 1 
instead. Putting everything together and multiplying by the appropriate symmetry factors we get 
(c2)a + 4(^2)5 + 2(c2)c + (c2)d = 17 as a total contribution at order 2. 

In view of more complicated graphs, an alternative and more convenient way of doing the calcu- 
lation in example 3 is by focussing on the probability distributions of the vertices lying on the lower 
space-like slice of the triangulated strip (see also Fig. 2.7(a)). The probability of having n space-like 
links in between the two vertices at the ends of the inverted V is given by the probability that there 
are n — 1 vertices in between with precisely one outgoing link each, which is l/2"~^. Thus we obtain 
again 

oo ^ 

(c2;G)e = 2^— ^=4. (2.22) 

Example 4. Consider an open chain of length greater than 2 embedded in such a way that only the 
first and last edges lie along time-like links, with a horizontal chain of m space-like links in between. 
Both time-like links are supposed to lie in the same strip of the triangulation, as illustrated in Fig. 
2.7(6). We are going to prove that the lattice constant of such a configuration is given by 

The result 4 — 2^"™ could be obtained easily by (i) considering the number of ways in which each 
of the two vertices at the top left and right can have an incoming link (giving a total of 2 x 2 = 4 
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Figure 2.7: Two configurations contributing to the calculations of examples 3 and 4. The labels m and n count the 
numbers of space-like links as indicated. 

possibilities) and (ii) subtracting from that the probabihty for the two time-Hke Unks touching each 
other in their lower vertex, which is 2^"™ (the probability that the m — 1 intervening vertices have 

only one incoming link each). However, we rather want to do the counting in a way that keeps track 
of the number of links separating the two vertices at the bottom, for reasons that will become clear 
soon. In other words, we want to determine the lattice constant for a closed polygon consisting of two 
vertical edges, m edges along the top and n along the bottom. Prom a combinatorial point of view this 
can be rephrased as the problem of counting the number of different ways in which n triangles and m 
upside-down triangles can be arranged to form a strip, with the well-known binomial result {"^^")- 
To get the lattice constant for our random lattice we still have to include a probabilistic factor ^ for 
each vertical link in the polygon interior, leading to the factors 2^~"*~" in formula (2.23). 

The strips decomposition. The embedding typology {g)i of any connected graph g extends over a 
well-defined number of strips in the CDT triangulation, and each part of (g)j belongs to a definite strip 
or horizontal slice. Our first step in the counting procedure is to decompose the embedded graphs into 
strip contributions, thus breaking the graph into pieces. The example of Fig. 2.8 will help us explain 
the procedure. It shows a particular embedding of C4, how it is decomposed into two pieces, and how 
its lattice constant can be computed accordingly. All we have to do is multiply the probability of C2,c 
(example 2) with m links at the bottom with the probability of 0^-1-2,6 (example 3), and then sum 
over m, resulting in 

°° 1 on 
m=l 

We can think of this operation as a product of an (infinite-dimensional) vector and matrix, 

1 



i^i(m) = 



m—l ' 



1 m + n^ 
K2(m,n) = - — ; T- 



(2.25) 



2" 

followed by taking the trace in order to impose the open boundary condition, that is, 

(C4). = Tr{Ki ■ K2). (2.26) 

Fig. 2.8 is a simple case with only two vertical lines in each strip. In general there will be more 
lines, like in Fig. 2.9. However, there are only two possibilities for pairs of neighbouring vertical links: 
they can be either disjoint or have exactly one vertex in common, so that the probability associated 
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n 



Figure 2.8: Example of a strip decomposition. The embedding on the left is split into two strips, after which the strip 
contributions for given m are multiplied pairwise and summed over. 

with them in a particular embedding in a particular triangulation will be given by either K2{m, n) or 
Ki(n). We can then associate a probability to a given pattern of vertical lines in a strip with fixed 
distance between them (in terms of the horizontal links in the in- and out-slice), which is the product 
of the corresponding K\- or ii'2-probabilities, times a factor 5 for each shared vertical line between 
two neighbouring /fj-patterns in the same strip (see Fig. 2.9). In the end we have to glue the strips 
back together again, with conditions such as to get the desired graph, and sum over the allowed values 
of horizontal lengths. As an illustrative example, we obtain for Fig. 2.9 the multiple sum 

00 00 00 m2 — 1 00 ^ ^ ^ 

E E[^2(«i+ii,mi)-/i:2(/2,m2)][ifi(mi)-K2(m2-i2,n)-Jri(i2)], (2.27) 

ji = lmi=lm2=2 J2=l n=\ 

where we have grouped the factors corresponding to the same strip in square brackets. As should be 
clear from this example, the gluing and the range of the summations for the various graph embeddings 
will have to be taken care of on a case-by-case basis. 




Figure 2.9: Example of a more complicated strip decomposition where eax;h strip consists of a sequence of elements 
of the kind illustrated in Fig. 2.7. The sum refers to the total sum in the expression (2.27) which has to be performed 
after the proper weights have been assigned to each element of the graph. 

Using these techniques we have repeated and independently confirmed the calculation of the first 
five orders given in [54] and then extended it to order 6. Table 2.1 gives the lattice constants (times 2) 
per number of vertices (which we will call the susceptibility coefficients) of the open chains, the other 
graphs, and of their grand total. The effort involved in this extension is by no means minor. There are 
387 order-6 open chains that have to be counted individually, with the identifications introduced in 
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this section, and for many of them the computation requires a careful use of the strip decomposition 
and a subsequent evaluation of the sums involved. 



n 


open 


closed + disconnected 


total (a„) 


1 


6 





6 


2 


34 





34 


3 


174 





174 


4 


859i 


-12 


847i 


5 


4152| 


-162| 


3990 


6 


19800^ 


-1416i 


18384^1 



Table 2.1: Susceptibility coefficients for the CDT lattice up to order 6. The terminology 'open', 'closed' and 'discon- 
nected' refers to the topology of the graphs. 



2.3.4 Counting graphs on the dual CDT lattice 

On the dual lattice (an example of which is depicted in Fig. 1.6, right-hand side) we can apply formula 
(2.15), which fov q = 3 reduces to 

ai - 4aa;-i + 4aj_2 = 2{l - 2)di-2 - 2ldi + 

2.28 

+ 8 [mi{l, 1) + 2mj_i(l, 1) + m(_2(l, 1)] . 

Because of the low coordination mimbcr the only even graphs that can appear are connected or 
disconnected polygons, and the only contributions to mi are from 6'-graphs and dumbbells with or 
without disconnected polygonal components. For the disconnected graphs we can use the reduction 
theorem. It should be noted that in the overlap decomposition new topologies of graphs appear, 
for example, graphs with more than two odd vertices, still rendering the higher-order computations 
non-trivial. 

It is not difficult to derive the probability distribution relevant for the dual lattice from the one of 
the original lattice. The probability of having two horizontal links at a vertex, one to the right and 
one to the left, is again 1. The probability for an incoming or outgoing vertical link at a vertex of 
a given chain of horizontal links is given by the probability of having a triangle or an upside-down 
triangle in the corresponding strip of the original triangulation. Since the two cases are mutually 
exclusive, the in- and out-probabilities are not independent and each takes the value 1/2, which may 
also be thought of as the weight associated with a triangle in the original lattice. With these rules 
the lattice constant of an open chain of length 1 is 3/2: 1 for the horizontal link (associated to its 
left vertex, say, in order to avoid overcounting), plus 1/2 for the vertical link. Alternatively, this easy 
example can be compiitcd by use of the Eulcr relation. For higher-order graphs it is often convenient 
to go back to the original lattice and count - with the appropriate weights - the configurations which 
can be associated to the dual graph under consideration. 
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Example 1. Let us consider again a closed polygon with only two vertical links, with m links on the 
top and n on the bottom, but this time on the dual CDT lattice (see Fig. 2.10). A careful evaluation 
of the probabilities involved yields 




{Cm,nr = 1^ (.11.1 (2-29) 

for its lattice constant, where k is the number of vertical links of the lattice which lie inside the 
polygon, and (•)* denotes the average lattice constant per vertex of the dual graph. The binomials in 
(2.29) count the number of possible ways to arrange the k internal links, and the remaining probability 
- expressed as a power of 1/2 - turns out not to depend on k. 




Figure 2.10: An embedding of the polygon considered in example 1. 

As already mentioned, the in-and out-probabilities at a vertex of the dual lattice are not indepen- 
dent, which means that we cannot apply the strip decomposition of the original lattice directly, but 
instead have to proceed more carefully. 

Example 2. As an example of this, consider a polygon embedding that extends over two strips (Fig. 
2.11), a type of configuration that appears from order 8 onward. Its lattice constant is given by 

mi-l / ^ _ -|^\ m2-l / ^ _ -|^\ 



fcl=0 \ ^ / fc2=0 

min{ki ,ni+n2— 4) 

X E 

i=0 



fci - I + ^2 \ / ni + n2 - 2 \ 1 

II ^ I 2'"i+'"2+ni+n2+fci— i+fc2+2 ■ 



In eq. (2.30), mi and m2 count the links on the top and bottom horizontal lines, and ni and n2 count 
the two sets of contiguous links on the central horizontal line. The numbers of internal vertical links 
in the upper and lower strip are denoted by ki and k2, while i counts how many links out of the fci 
end on one of the two intermediate horizontal lines. The logic behind the various counting factors 
appearing under the sums is very similar to that of the previous example. 

We have found the method illustrated by examples 1 and 2 the most compact one to deal with the 
counting of graph embeddings on the dual CDT lattice. It has enabled us to push our "counting by 
hand" to order 12, the results of which arc given in Table 2.2. Because of the fixed, low coordination 
number of the dual vertices, there are far fewer graphs at a given order than there are on the original 
CDT lattice. 
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Fi gurc 2.11: An embedding of the polygon considered in example 2. 
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12 


23 


42+1 


78+i 


142+1 
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461+i 


820+i 


1446+i 


2532+ii 



Table 2.2: Susceptibility coefficients for the dual CDT lattice up to order 12. 



2.4 Series analysis 

2.4.1 Review of different methods 

Having computed the two high-temperature series to some order, we will now turn to analyzing them, 
and try to extract information on the critical properties of the underlying gravity-matter systems. 
The standard analytic methods arc well-known from the case of regular lattices and reviewed in [82] . 
We will recall them briefly here in order to be self-contained. 

The first and most straightforward tool is that of the ratio method, which works as follows. As- 
suming a simple behaviour of the susceptibility of the form 

x{u)r^A{u)(l--y'+B{u) (2.31) 

near the critical point Uc, with analytic functions A and B, its series expansion x{u) = 1 + X^„>o ^n'"" 
should yield (to first order in 1/n) 

r„.^ = 1(1 + 1^). (2.32) 

One can then generate sequences of estimates of the critical parameters from sequences of point pairs 
{r„,r„_i} [82], namely, 

n(2 - n)r„ + (n - l)Vn-i 

7n = 7 r\ ' (2.33) 

nr„ - (n - l)r„_i 

Uc,n = r^-^T ■ (2.34) 

nrn - [n - l)rn-i 

In general, these will converge very slowly. In addition, if the series has an oscillatory behaviour 
the method can yield alternating over- and underestimates. In this case a fit of a whole sequence 
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{fmim •••) '"maa;}) Computed with the help of (2.32), may be better suited than the sequence of estimates 
to find the straight Hne masked by the oscillations {rmin is chosen properly to exclude large deviations 
from (2.32) at small n, and rmax is the ratio computed for n = Umax)- 

Oscillations and other irregularities in the expansion are eaused by a more complicated behaviour 
of the thermodynamic function, for example, the presence of other singularities in the complex plane 
close to the physical singularity, which in unfortunate cases may even lie closer to the origin than the 
singularity of interest. In case the behaviour near the physical singularity is like 

x{u)r^A{u)(l--y\ (2.35) 

with A{u) a meromorphic function with singularities close to Uc, a method that should give better 
results than the ratio method is that of the so-called Pade approximants. The method consists in 
the approximation of a function, known through its series expansion to order A/", by a ratio of two 
polynomials of order C and M., subject to the condition C -\- M. < M , 

^% ^ ^ ^yj;'""' , = Fu{x) + 0(x^+^+i) (2.36) 

with Fj^{x) = Yl'k=o '^kX^- By common usage, the notation [L/M] indicates the order of the polyno- 
mials used. 

In our specific case one takes as the function F{x) to be approximated the derivative of the 
logarithm of x(^^)) so that Uc can be recognized as a pole and 7 as the associated residue in 

i.logx(w) = ^— (l + 0(«e-w)). (2.37) 
au Uc — u 

This is also referred to as the Dlog Pade method. Usually one only looks at the tridiagonal band 
[{M — 1)/A/], [M /M] and [(TV + 1)/T\/] because of the invariance of the diagonal Pade approximants 

under Euler transformations. 

The Fade approximants work well only when there is no additive term B(u). In presenee of such 
a term there is a generalization known as differential approximants - to account for functional 
behaviours of the form 

n _ 

X{u) --)"" +B{u), (2.38) 

with B{u) and ^i, -.-j^n analytical functions and ui, a set of singular points. If we are consid- 

ering the logarithmic derivative of a function f{x) we can rewrite (2.36) as 

Pc{x)f{x) - QM{x)r{x) = 0(x^+^+i). (2.39) 

The idea of the differential approximants method is to generalize this equation according to 

RM.{x)Fj^{x) + QmAx)F'j^{x) - Pc{x)Fu{x) = Sk{x) + o[x'^+'^+^^+M.+i)^ (2.40) 

where Rm2 {x) and Sic{x) arc two other polynomials of order M2 and /C, and where one can substitute 
the differential operator by _D = .t^, forcing the point at the origin to be a regular singular point. In 
the following we will only consider the special case Rm2 = 0, giving rise to the so-called inhomogeneous 
Ist-order differential approximant, denoted by [IC/£;Mi]. With this method, the exponent 7 can be 
evaluated as 

7=^^, (2.41) 

Xc^Ji4^ \Xc) 

where Xc is a simple root^ of Qmi{x), which gives an estimate of the critical point. 



The case of multiple roots can also be considered, but in our case never occurs. 
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2.4.2 The CDT lattice series 

The high-temperature series expansion for the Ising susceptibihty on the CDT model is given by 
Xcdt{u) = 1 + 6m + 34u^ + 174u^ + (847 + ^)m^ + 3990w^ + (18384 +^)u^ + O(w^). (2.42) 

Applying formulas (2.33) and (2.34) we get the sequence of estimates for 7 and Uc reported in Table 2.3. 
The results of linearly fitting the whole sequence of ratios {r2, Vmax} instead are listed in Table 2.4, 



n 




7 


3 


0.2489 


1.819 


4 


0.2424 


1.721 





n.24o!) 


1.7!)1 


6 


0.2438 


1.740 



Table 2.3: Two-point linear extrapolations of the critical point and the critical exponent of the series (2.42) obtained 
using (2.33) and (2.34). 

and a plot of the ratios {r2, rg} versus 1/n together with their fit is shown in Fig. 2.12. Despite the 
shortness of the sequence, the linear fit is of very good quality. It is worth noting that if we remove 



^max 


Uc 


7 


3 


0.2488 


1.820 


4 


0.2462 


1.789 


5 


0.2458 


1.783 


(i 


().24r,4 


1.779 



Table 2.4: Linear fits of the sequences {r2, ■■■,rmax} for the series (2.42) with assumed functional form (2.32). 

the point r2, which of course is subject to the largest deviations from a pure 1/n-behaviour due to 
higher-order terms, and fit the truncated sequence {rs, ...,r6}, we obtain 7 = 1.749 and Uc = 0.2441 
(see Table 2.5). 

It is surprising how fast the ratio sequence seems to converge to the rcsiilt 7 = 1.75 of the exact 
solution for the flat regular case, as compared to the analogous one obtained from the scries expansion 
on the plane triangular lattice (see [83] for the expansion coefficients to order u^^), even though we 
cannot claim the result to be conclusive because of the limited number of terms at our disposal (see 
Fig. 3.72). 

What is the picture when we use one of the alternative methods to evaluate the scries expansions? 
Results from the Dlog Pade and differential approximants methods are somewhat inconclusive, most 
likely because of the small number of estimates we can perform with our 6th-order series (see Table 2.6). 
To get a feel for what may be expected at this order, we report in Table 2.7 the corresponding Dlog 
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r 




0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50 0,55 
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Figure 2.12: Plot of the ratios (2.32) for the case of the CDT lattice (data from Table 2.4). 
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Table 2.5: Linear fits of the sequences {rs, ...,rmax} for the series (2.42) with assumed functional form (2.32). 
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Figure 2.13: A comparative plot of estimates of the critical exponent 7 for the Ising model on the CDT versus that 
on a plane triangular lattice (data taken from [83]), obtained via the ratio method. The filled shapes refer to the CDT 
model and the open ones to the plane triangular lattice. The squares indicate fits of {rn^i„, ■■■,rn} sequences, and the 
circles those of sequences {r„^^^+i, ...,rn} with the first ratio eliminated. 

Fade approximants [(-A/' + i)/-A/l to the series for the plane triangular lattice (see again [83]). The 

estimates for the critical exponent 7 in the latter case are clearly closer to the known exact value, 
whereas for the CDT model, with only two values for the diagonal [J\f/J\f], it is quite impossible to 
extrapolate the behaviour of 7. (Note that we are reporting the estimates for the critical point Uc 
only for completeness; they are or course not expected to coincide for the two models.) We have also 
computed the inhomogeneous first-order differential approximants [jC/{J\f + j);J\f] for the CDT series, 
which was defined just before formula (2.41). Values for the critical exponents are not completely off, 
but there are simply too few of them to make any statement about their convergence behaviour (see 
Table 2.8). 
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Table 2.6: Dlog Pade approximants method applied to the CDT series (2.42). 
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Table 2.7: Dlog Pade approximants method applied to the plane triangular series in [83]. 
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Table 2.8: Inhomogeneous Ist-order diflferential approximants {C/(M + j);M] method for the series (2.42). 
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2.4.3 The dual CDT lattice series 

The high-temperature series expansion for the Ising susceptibility on the dual CDT model is given by 

XcDTd(M) = 1 + 3w + + 12u^ + 23w^ + (42 + ^)u^ + (78 + ^)u^ + (142 + 

+258^8 + (461 + + (820 + l)u^° + (1446 + + (2532 + ^)u^^ + 0{u^^) 

16 8 32 16 

(2.43) 

Using expressions (2.33) and (2.34) we get the sequence of ratio- method estimates for 7 and Uc reported 
in Table 2.9. 



^max 


Uc 


7 


3 


0.5 


1. 


4 


0.6 


1.6 


5 


0.6147 


1.713 


6 


0.5800 


1.390 


7 


0.5842 


1.436 


8 


0.5782 


1.360 


9 


0.6057 
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0.6064 


1.769 


11 


0.6093 


1.820 


12 
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2.033 



Table 2.9: Two-point linear extrapolations of the critical point and the critical exponent of the dual series (2.43) 
obtained using (2.33) and (2.34). 



The relative ratio plot is shown in Fig. 2.14, and shows gentle oscillations around the best linear fit of 
the sequence {ra, ri2}, leading to the estimates 7 = 1.568 and Uc = 0.5956. However, these numbers 
do not tell it all. The strong qualitative differences with the planar honeycomb lattice, the appropriate 
non-random version of the dual CDT lattice, are illustrated in Fig. 2.15. It shows a comparative plot of 
the ratios extracted for the two models (the expansion coefficients for the honeycomb lattice, known 
to order u^"^, can be found in [83]). Obviously, the influence of interfering unphysical singularities 
present on the regular honeycomb lattice is much reduced in the fluctuating dual CDT ensemble. 

On the basis of the rather well-behaved results for the Ising susceptibility on dynamically triangu- 
lated spacetimes we conjectured in [97] that "coupling to quantum gravity" may be an optimal method 
to learn about the critical behaviour of a matter or spin system, in the sense that the randomness of 
the underlying geometry eliminates spurious unphysical singularities (due to the presence of discrete 
symmetries in the case of regular lattices), but is not strong enough to change the universality class 
of the matter system. This is a dynamical version of similar conjectures made in the 1980s, which ad- 
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Fi gure 2.15; A comparative plot of the ratios of the regular honeycomb lattice (data taken from [83]) versus the dual 
CDT model. 
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vocated the use of fixed random lattices to improve the convergence behaviour of lattice gauge theory, 
say (see, for example, [84]), but ultimately did not succeed. 

If it is indeed the case that the singularity structure of the thermodynamic quantities in the com- 
plex plane is simplified, this could explain that the simple ratio method does indeed give the best 
result, at least at the relatively low orders we have been considering. In the case of the Ising model 
coupled to Euclidean dynamical triangulations, an analogous result has already been obtained with 
regard to the locus of the zeros of the partition function in the complex-temperature plane, where 
unphysical singularities of the flat regular lattices have been shown to be largely absent [85]. However, 
as already mentioned, the underlying geometries are too disordered to serve our present purpose, be- 
cause their critical matter behaviour is changed compared to the flat case. 
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Table 2.10: Dlog Pade approximants method for the dual CDT series (2.43). 
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Table 2.11: Dlog Pade approximants method for the honeycomb lattice series in [83]. 

Returning to the evaluation of the dual CDT series, we have used both the Dlog Pade and the 
differential approximants methods, in addition to the ratio method already described. Even with the 
larger number of terms compared to the original CDT case the results are reluctant to show any clear 
sign of convergence. Our results for the Dlog Pade approximant are summarized in Table 2.10. The 
corresponding data for the regular honeycomb lattice [83] , which we are including for comparison in 
Table 2.11, seem to show more consistency at the same order of approximation. The result from using 
the inhomogeneous first-order differential approximants for the dual CDT lattice (see Table 2.12) is 
similarly inconclusive. It would clearly be desirable to understand in more analytic terms the influence 
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Table 2.12: Inhomogeneous 1st order differential approximants [C/N -\- method for dual CDT series (2.43). 
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of the fluctuating geometry of the underlying quantum gravitating lattice on the singularity structure 
of spin models like the one we have been considering, and thus determine which of the approximation 
methods is best suited for extracting their critical behaviour. 

2.5 Comments on the low-temperature expansion 

Having dealt with the high-temperature expansion, it is natural to ask whether similar expansion 
techniques can be applied in a low-temperature expansion of our coupled models of matter and gravity, 
in the same way this is possible for spin systems on regular lattices. The first thing to notice is that 
at the point of zero temperature around which one expands, the spins are all frozen to point into the 
same direction, and effectively play no role. The quantum-gravitational ensemble is therefore again 
characterized by the probability distribution (2.6) for the coordination number of pure gravity, as was 
also the case in the limit of infinite temperature. 

In the low-temperature expansion on the regular lattice one starts from a ground state with all 
spins aligned, and then includes perturbations with s overturned spins, obtaining [86] 

lnZN{K,H) = IK + vH + Y,[s;r;G]zPn', (2.44) 

in terms of the expansion parameters z = e^^^ and = e~^^, where K, H, I and v were all defined 
at the beginning of Sec. 3 above. The part linear in N of the (strong) lattice constant (sec [79] for a 
definition) for a graph with s vertices and r links embedded in G is denoted by [s; r; G], and p is the 
number of lattice links which are incident on any of the vertices of the graph, but do not themselves 
belong to the embedded graph. On a lattice of fixed coordination number q it is easy to prove that 
p = qs — 2r. Graphs with given values of s and r therefore contribute only at a specific order. 

The low-tcmpcrature expansion on CDT lattices via graph counting represents a slight compli- 
cation. It is easy to sec that already by overturning just a single spin the number of links whose 
interaction changes sign as a result depends on the spin's position, so that each site will contribute to 
a different order in the perturbative expansion.^ In the case at hand, this diflaculty can be overcome 
by counting the number of elementary polygons (i.e. faces) with given boundary length p on the dual 
lattice. In the case of several overturned spins we can proceed similarly by associating certain patterns 
on the dual lattice with a given order of z. More precisely, to obtain the coefficient at order z^/U^ we 
need to count polygons on the dual lattice which have total boundary length p and are constructed 
out of s faces. Denoting the number of such patterns by 'P{p, s), we can write the free energy per unit 
volume as 

1 In Zn{K, H) = ^K+ + ^Vip, s)zPii' (2.45) 

and the susceptibility at vanishing magnetic field as 

x{z)=^4s'V{p,s)zP. (2.46) 

s,p 

Although the procedure is slightly more involved, and cannot be seen as a computation of (averaged) 
strong lattice constants, the expansion is still doable using the counting techniques introduced earlier in 
this chapter. Table 2.13 gives the results for 'P{p, s) up to order 10. The susceptibility low-temperature 
scries (2.46) turns out to be very irregular and none of the methods of analysis considered seems to 

^By contrast, for the low-temperature expansion on dual CDT lattices, where the coordination number is fixed to 3, 
the expansion parallels that on flat lattices, with the (strong) lattice constants replaced by their ensemble averages. 
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Table 2.13: Free energy coefficients 'P{p,s) as defined in (2.45) for the CDT lattice. 

give a reasonable indication of the critical exponent, at least not from the relatively few terms we 
have computed. Although this is somewhat disappointing, in view of the fact that low-temperature 
expansions for regular lattices are notoriously ill-behaved, it does not come as a total surprise. 

What helps in the evaluation of the low-temperature series for the Ising model on flat triangular, 
square or honeycomb lattices is the fact that their critical temperatures are known exactly from 
duality arguments and thus can be used as an input in so-called biased approximants to considerably 
improve the estimates of critical exponents. Namely, for H = one can apply the standard duality 
transformation [75] which maps the low-temperature expansion (2.44) of the triangular lattice model 
to the high-temperature expansion (2.3) of its dual and vice versa. This can be combined with the 
star-triangle transformation to obtain the critical point. The analogous deriviation for the simple 
square lattice is even easier, because it is self-dual. Unfortunately, the star-triangle transformation is 
not applicable in the CDT case, which prevents us from making a similar argument in the coupled 
Ising-gravity model. In order to pursue the analysis of the low-temperature series further, we would 
therefore have to rely on the evaluation of higher-order terms in the expansion. 

2.6 Summary and outlook 

In this chapter we have given a concrete example of a method which can be used to extract physical 
properties of a strongly coupled model of gravity and matter. We showed how the method - estimating 
critical matter exponents from a series expansion of suitable thermodynamic functions of the system 
- can be adapted successfully from the case of a fixed, flat lattice to that of a fluctuating ensemble of 
geometries, as is relevant in studies of non-perturbative quantum gravity. For the ensemble of two- 
dimensional causal dynamical triangulations, we have formulated an explicit algorithm for counting 
embedded graphs which allows us to do the counting recursively for increasing order. The method 
can in principle be applied to other matter and spin systems which admit a similar diagrammatic 
expansion in terms of weak or strong lattice constants around infinite or zero temperature, and where 
one has sufficient information about the probability distribution of local geometric lattice properties 
like the vertex coordination number. Even in cases where these lattice properties are not available 
explicitly, the lattice constants could still be extracted from simulating the pure gravity ensemble. 
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As a potential spin-off, we noticed that the explicit expansion for the susceptibility in vanishing 
external field of the Ising model on CDT (up to order 6) or dual CDT lattices (up to order 12) 
indicates a more regular behaviour of this function in the complex-temperature plane than that of 
the Ising model on the corresponding fixed triangular or honeycomb lattice. This may also explain 
why the simple ratio method, applied to CDT lattice results from only six orders, gives an excellent 
approximation to the Onsager susceptibility exponent 7 = 1.75. As we saw in some detail, other 
approximation methods, namely, Dlog Fade and differential approximants, do not produce results 
of a similar quality. We believe that for the case of the dual CDT lattice, the series computed is 
simply too short to yield reliable estimates with any of the approximation methods (even at order 
12, the number of terms contributing is much smaller than the number of terms contributing at 
order 6 on the original CDT lattices). At any rate, if one aims to make an argument of improved 
convergence of matter behaviour on fluctuating lattices (see also [97]), it is plausible that this will 
be achieved optimally with the triangulated CDT geometries whose vertex coordination number can 
vary dynamically from 4 all the way to infinity, rather than with the dual tesselations which have a 
fixed, low coordination number of 3. To investigate the convergence issue in more detail will require 
going to higher order than 6 on CDT lattices, which is not really feasible 'by hand', as we have been 
doing so far, but will require the setting up of a computer algorithm to perform (at least part of) the 
graph counting. 

Our work should be seen as contributing to the study of non-perturbative systems of quantum 
gravity coupled to matter, which is only just beginning. One challenging task will be to establish 
computable criteria characterizing and quantifying the influence of geometry on matter and vice versa, 
in the physically relevant case of four spacetime dimensions, something about which we currently know 
close to nothing. Obvious quantities of interest are critical exponents pertaining to geometry (like, for 
instance, the Hausdorff and spectral dimensions of spacetime already measured for the ground state of 
four-dimensional CDT [64]), and critical matter exponents, like that of the susceptibility investigated 
in the present work. One would like to have a classification of possible universality classes of gravity- 
matter models as a function of the characteristics of the ensemble of quantum-fluctuating geometries. 
A necessary condition for viable quantum gravity models is that at sufliciently large distances and 
for sufficiently weak matter fields the correct classical limits for these critical parameters must be 
recovered. 

Two-dimensional toy models like the one we have been considering can serve as a blueprint for 
what phenomena one might expect to find. In two dimensions, there have been several investigations 
of the eff'ect of random geometric disorder on the critical properties of matter systems, and attempts 
to formulate general criteria for when a particular type of disorder is relevant, i.e. will lead to a 
matter behaviour different from that on fixed, regular lattices. Good examples arc the Harris [87] 
and the Harris-Luck [71] criteria, which tie the relevance of random disorder to the value of the 
specific-heat exponent a and to correlations among the disorder degrees of freedom. However, not 
all models fit the predictions, and open problems remain (see [73] and references therein). Causal 
dynamical triangulations coupled to matter add another class of models inspired by quantum gravity, 
whose randomness lies in between that of Poissonian Voronoi-Delaunay triangulations and the highly 
fractal Euclidean dynamical triangulations, and over which one has some analytic control. It will be 
interesting to see to what degree the robustness of the matter behaviour with respect to the geometric 
fluctuations observed so far [54, 55, 97, 98] will persist for different types of spin and matter systems 
and in higher dimensions. 
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Chapter 3 



Moving to three dimensions 



In this chapter I will present work on a three-dimensional model of CDT. The model I will consider is 
a constrained version of the general one reviewed in the first chapter. Rather than allowing all possible 
triangulations of the spatial slices and of the spacctimc in between two of them, only a particular class of 
triangulations will be allowed, in a similar fashion to the construction leading to two-dimensional CDT from 
Euclidean DT. The restriction will permit an analytic evaluation of the partition function for the one-step 
propagator and of its continuum limit. A quantum Hamiltonian for this model will be derived. This chapter 
is based on [99]. 

3.1 The challenge of higher dimensions 

Although the CDT model has been solved exactly in (l-l-l)-dimensions, it turns out to be very difficult 
to extend this treatment to higher dimensions, for which most results so far have come from numerical 
simulations. An exception to this is the matrix- model formulation of the (2-|-l)-dimensional model 
[57]. An exact solution of the relevant matrix model has been given [60], but its rather complicated 
and implicit form has so far been an obstacle to performing its continuum limit analytically, and 
extracting the quantum Hamiltonian. 

This situation is not particularly surprising, as most of the known solvable statistical models are 
only one- or two-dimensional. Even a simple spin model like the Ising model has been solved only in 
two dimensions, and for zero external field. One can therefore already anticipate^ that an extension of 
analytical methods and results to the case of higher-dimensional CDT models will be a challenging 
task. One possible way to approach this problem is to try to modify the model by making simplifying 
assumptions to improve its solubility, but without changing its physical content in the continuum. 
The main aim of the present work is to study a reduced model of three-dimensional CDT, which is a 
particular case of a class of models introduced in [89] . This model possesses some additional "order" , 
which will enable us to bring a number of technical tools to bear on its solution. As we will see, this 
leads to a non-trivial dynamics for the three-dimensional quantum universe, justifying our original 
"truncated" ansatz. 

This chapter is organized as follows. In Sec. 3.2 we describe the geometric configurations whose 

partition function we are going to study and use as a transfer matrix. The configurations are from an 
ensemble of triangulations of product type, interpolating between successive constant-time slices, with 
free boundary conditions and boundary cosmological constants, conjugate to the areas of the slices. 
In Sec. 3.3 we introduce the first tool that will help us in solving the three-dimensional CDT model, 
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the inversion formula for heaps of pieces. In Sec. 3.4 we reformulate the partition function in terms of 
the transfer matrices for a set of one-dimensional hard bi-coloured dimer models. In Sec. 3.5 we study 
the location of the singularities of the partition function and define a critical surface in the parameter 
space, which is the boundary of the region where the model is convergent and well defined. In Sec. 3.6 
we introduce the second tool we will need for the solution, the replica trick for products of random 
matrices. This enables us to identify a critical point on the critical surface where the mean volume 
goes to infinity, thus making it a potential candidate for taking a non-trivial continuum limit. We can 
give an exact solution to the model at the critical point, and approximate it to the desired order in 
the displacement parameter when moving away from it. In Sec. 3.7 we start exploring the continuum 
properties of the model by deriving expectation values for the volume and curvature of spacctime. In 
Sec. 3.8 we construct the area-to-area transfer matrix by summing over unphysical degrees of freedom, 
and explain the procedure for extracting information about the spatial volume of the universe. In 
Sec. 3.9 we show that for a vanishing bare inverse gravitational constant, one obtains a well-defined 
quantum Hamiltonian acting on the Hilbert space of the area eigenstates. By contrast, a canonical 
scaling ansatz for the gravitational constant does not seem to lead to a meaningful result. In Sec. 3.10 
we consider a possible extension of these calculations, encoding further global metric information in 
the form of an additional Teichmiiller parameter. Finally, our conclusions are presented in Sec. 3.11. 
The four appendices at the end of the thesis collect technical results and discussions needed in the 
main text of this chapter. 

3.2 Introducing the model 
3.2.1 The product type triangulations 

Part of our strategy for trying to solve the non-perturbative three-dimensional model of quantum 
gravity defined by causal dynamical triangulations (CDT) is to identify a suitable (sub-)class of all 
CDT configurations whose superposition can be tackled analytically. The aim of such a reduction is 
to simpliiy the analytic treatment, without eliminating relevant degrees of freedom to such a degree 
that the universality class is changed, compared to that of the full CDT model. ^ To achieve this, 
we propose to work with a set of triangulated, causal geometries which have an additional "order" 
imposed on them. This order is mild in the sense of not implying "isometrics" of the triangulations; 
both the local spacetime curvature and the local curvature of spatial slices can still vary arbitrarily. 
Our model therefore has less order than the hexagon model considered in [88], which has flat spatial 
slices. 

The way in which we will introduce more structure or "order" on triangulated spacetimes - in- 
spired by similar earlier ideas [50, 89] - is that of applying the building principle inherent in causal 
dynamical triangulations twice over. This "causality principle" recognizes that in a metric spacetime 
of Lorentzian signature not all directions are equivalent, but there is a distinguished (class of) time 
direction(s), in line with the existence of light cones and causal relations (none of which are present 
in spaces of purely Euclidean signature). As we have seen in Chapter 1, in CDT this causal structure 
is implemented via a discrete global time slicing. This simply means that each simplicial building 
block of a triangulated spacetime must be contained in exactly one spacetime "sandwich" (the region 

'^This assumes that the full CDT model leads to an essentially unique three-dimensional quantum gravity theory. 

Strictly speaking, rather little is known in dimension three about the universality classes of statistical models of random 
geometry like the one we are using. There are certainly more, but we are only interested in those which possess a 
continuum interpretation in terms of quantum gravity. 
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between discrete proper times t and t + 1). This implies a causal ordering on the simplices, without 
constraining the local curvature degrees of freedom. 

In addition to this physically motivated choice of a distinguished time direction, which is one of the 
key ingredients of the approach of causal dynamical triangulations, we will introduce here a second 
distinguished direction, but this time purely for computational convenience, and under the assumption 
that it will not affect the universal properties of the gravitational model^. This will roughly speaking 
correspond to an additional slicing, but this time in one of the spatial directions. As explained in [89], 
the resulting structure can be thought of as a staggered fibration in a piecewise flat setting. 

To illustrate the main idea behind this type of triangulation let us start from the easier (1+1)- 
dimcnsional CDT model. A spacctimc contributing to the sum over geometries is usually described 
as a sequence of triangulated strips glued together, where each strip represents a piece of spacetime 
between proper times t and t + 1, and the lengths of the spatial boundaries of adjacent strips must 
match pairwise. For our purposes, it is useful to think of a strip as being constructed "sideways" 
(see Fig. 3.1): starting from a segment (a time-like link), build a frame on it and then fill the frame 
with some sequence of up- and down-triangles. We will refer to this construction as "building a 
two-dimensional tower over a link" . The entire two-dimensional triangulated spacetime can then be 
regarded as fibration over a one- dimensional chain of links. 




Fi gure 3.1 1 A triangulated strip constructed as a tower over a one- dimensional link. 

More generally, product triangulations [89] are obtained by building towers over (arrays of) higher- 
dimensional simplices, such as triangles. How a three-dimensional tower is built over a triangle is 
illustrated in Fig. 3. 2."^ If the base space is not just a single triangle, but a two-dimensional triangula- 
tion consisting of n triangles, we can construct a three-dimensional product triangulation by erecting 
a tower over each of them, in such a way that the boundary triangulations of the resulting prisms 
again match pairwise. A general (n -I- m)-dimensional product triangulation is a simplicial manifold 
constructed by consistently building (fc + m)-dimensional towers (for all k < n) over the fc-dimensional 
sub-simplices of an n-dimensional base simplicial manifold. 

In the present work we will use this construction in the context of (2-|-l)-dimensional causal 
dynamical triangulations. By definition, these spacetimes can be regarded as product triangulations 
whose base space - like in (1+1) dimensions - is a one-dimensional triangulation in the time direction, 
and the fibres over each link of the base space are the sandwiches between integer values of the 
discrete time t introduced earlier. Now, instead of allowing sandwich geometries which are arbitrary 
triangulations of thickness At = 1, we will impose an additional product structure on them. Namely, 

^An assumption that obviously will have to be borne out by the final result 

'^It should be noted that the figures do not represent faithfully the metric properties of the building blocks, that is, 
their edge lengths. By construction, up to a relative factor between space- and time-like links, the edge lengths in a 
causal dynamical triangulation are all identical. 
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Figure 3.2: A triangulated prism constructed as a tower over a two-dimensional triangle. 

they should themselves have the form of a sequence of up- and down prisms, as depicted in Fig. 3.3 
below. That is, in addition to the slicing of the entire spacetime corresponding to the global discrete 
proper time, each sandwich possesses a discrete slicing in a given spatial direction. The number of 
slices (that is, the number of prisms in a sandwich) is allowed to vary from sandwich to sandwich. 

An equivalent way of characterizing such spacetimes is to regard them as fibrations over a two- 
dimensional base space (with one space- and one time-dimension) which is itself an arbitrary (1-|-1)- 
dimensional CDT, and where over each triangle of the base space we have erected a tower "filled" 
with a sequence of tetrahedra, in such a way that the triangulations on the faces of neighbouring 
towers match and can be glued together consistently. In this way, we may think of our triangulations 
as being (l-|-l-|-l)-dimensional. 

3.2.2 The partition function 

We will concentrate on the dynamics of a "sandwich geometry" , given by the transition amplitude 
from a spatial geometry at time t to one at time t + 1, in other words, the transfer matrix T of the 
causal dynamical triangulation (CDT) model, which in turn contains information about its quantum 
Hamiltonian H by virtue of the relation T = e~°'^ (see [47] for details). 

Since pure gravity in (2+1) dimensions docs not possess any local, propagating degrees of freedom, 
we expect that most details of the spatial geometries will be dynamically irrelevant, leaving only the 
spatial two-volume A{t) and the Teichmiiller parameters (together with their canonically conjugate 
momenta) as physical degrees of freedom. Therefore, rather than calculating the matrix elements 

G(r7i,52,At=l) = (52|T|<?i) (3.1) 

of the transfer matrix from an arbitrary spatial triangulated geometry \gi) at time t to an arbitrary 
\g2) at time t + 1, we will only keep track of the two boundary areas Ai and A2 and evaluate the 
reduced matrix elements 

G(^i,^2,Ai = l) = (^2|T|^i), (3.2) 
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where \Ai) is a normalized linear combination of states \gi) with given area Ai (see [57] for a detailed 
discussion). Note that one is still summing over the same sandwich geometries as in (3.1), so the 
reduced matrix elements still capture the effective dynamics of a// geometric excitations of the sandwich 
geometry.^ Once these matrix elements are known, a next step is to try to extract the continuum 
Hamiltonian operator H of the system from an expansion in the short-distance cutoff a as a ^ 
according to 

{A2\f\Ai) = {A,\e-''"\Ai) = {A,\ (l - aH + Oia^)) \A,). (3.3) 

We will denote by Nij the number of simplices having i vertices on the initial boundary spatial 
slice at time t and j on the final one at t + 1. In this way Ai (resp. A2) is given by N^i (resp. Nis). 
The prescription for evaluating the discrete one-step propagator (3.2) is given by 

G(iV3i,iVi3,At = l)= V ^e-^-«, (3.4) 

-'|JV31."l3 

where the sum is over all sandwich triangulations T with fixed boundary areas A^si and N13, Seh 
the Wick-rotated discrete Einstein-Hilbert action, and Cr the order of the automorphism group of 
T. Taking into account boundary terms in the action to ensure the correct propagator behaviour, we 
find 

SEH^a{Ni3 + N3i)+f5N22+lN (3.5) 

as an explicit expression for the gravitational action, where N is the number of triangle towers in the 
sandwich, and where we have introduced the parameters 

a = —cik + biX 

P = C2k + 62A (3.6) 
7 = csfc 

depending on the dimcnsionless bare cosmological and inverse Newton constants A and k, and on 
positive numerical constants Cj and bi characterizing the geometric construction (see appendix A-I). 

As one might have anticipated from previous investigations of related three-dimensional quantum 
gravity models, the evaluation of (3.4) remains a challenging task, also for our specific choice of 
ensemble of triangulations. Our main aim will be to calculate its discrete Laplace transform 

Z{x,y,At=l) = ^^x^-yiVi3G(^3^^^^3^Ai = 1) 

"""" (3 7) 

= ^^(a;e-")^^^(t/e-")^- E e'^^^^^-^^ 

N31 N13 7]jVi3,JV3i 

which can be thought of either as the generating function of G'(iV3i, iVi3, At = 1) or as the partition 
function of the sandwich geometry with free boundary conditions (thus summing over all values of 
the boundary volumes) and with additional boundary cosmological terms x = e"^'" and y = e"'^""* 
in the action. 

To get back the transfer matrix from Z{x,y,At = 1), one needs to keep x and y distinct and 
variable. However, as shown by other examples (see the end of section 1.3.2 for the two-dimensional 



*The presence of physical degrees of freedom beyond the spatial area depends on the topology of the spatial slices. 
There are none for the case of spherical slices considered in [62, 57], and there is one Teichmuller parameter in the 
present model, where we choose boundary conditions corresponding to cylindrical spatial slices. Work is under way 
to determine whether the calculation of matrix elements presented here is still feasible when the dependence on this 
parameter is kept explicitly. 
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case, and [62] for a numerical illustration in three dimensions), one may be able to extract non-trivial 
information about the phase structure of the model by considering the special values x = y = 1, which 
simplifies the evaluation of (3.7). 

At this stage we will exploit the special product structure of our chosen spacetime geometries to 
split the sum over all sandwich triangulations into two simpler sums. The idea is to first perform 
the sum over all fillings (i.e. triangulations) of towers for a fixed triangulated base strip, and then to 
perform the sum over all possible triangulations of the base strip (see Fig. 3. 3), schematically 

E = E E • (3-8) 

"^sandwich "^base strip "-^towers 





Figure 3.3: A "sandwich geometry" of product type in the (l-l-l-l-l)-dimensional CDT model, built over a given 
triangulated base strip. We split the calculation of its partition function into two parts, calculating first the partition 
function for such a geometry and then summing over all possible sequences of towers. By "towers of type (2,1)" we will 
mean the towers built on a triangle with two vertices at time t and one at time t + 1, and similarly for "towers of type 

(1,2)". Each of the towers is filled up with a sequence of tctrahcdra, as illustrated in Fig. 3.2. There are three different 
possibilities for how a tetrahedron can be oriented inside one of the prisms (sec [89] for more geometric details). 

We are thus naturally led to studying the partition function Zn for sandwich geometries of product 
type with N towers, and investigating its properties in the large- limit. In this sense we are treating 
the second distinguished direction, which defines the slicing within the sandwich geometry, in the 
same way as we treat the time direction in (l+l)-dimensional CDT. It is related to (3.7) via 

Z{x, y,At = l) = Y^ e-'^^ZAr(u, v, w, At = 1). (3.9) 

JV 

Introducing the weights ^/u = xe~", \po = ye~" and w = e~^, we can express Zn as 

ZN{u,v,w,At=l) = J2 E u^v^w'^'', (3.10) 

where the sequences ^jv consist of N prism towers, N — R of them of type (2,1) and R of type (1,2), 
for any R< N. 

3.3 Inversion formula 

The special product form of our sandwich geometries will enable us to make use of an inversion formula, 
relating its partition function to that of a model of one dimension less. As we have seen in Chapter 1, 
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this technique essentially provides a complete solution to the CDT model in dimension (1+1). This is 
not true in the much more involved case of CDT in dimension (1+1+1) we are addressing here, but - 
as we will demonstrate - it will nevertheless allow us to make substantial progress in the evaluation 
of the partition function. 

In order to apply the inversion formula to the case of three-dimensional causal triangulations, 
we will use the fact that the propagator of this model can be characterized in terms of geometric 
data which are "almost" two-dimensional. We have seen this characterization in the first chapter, 
when talking about the ABAB-matrix model. The procedure to be used here in order to obtain the 
two-dimensional dual graph describing the propagator is exactly the same, but the class of graphs we 
obtain is different, reflecting the difference in the class of triangulations. 

In the case of our product triangulations the dual graphs obtained have a special form. Namely, 
part of the dual graph can be drawn as a fixed sequence of straight lines (which we will call "vertical 
lines" , and which are drawn vertically in Fig. 3.4) coming from the prism towers, blue for a tower over 
a type-(2,l) triangle and red for a tower over a type-(l,2) triangle. The sequence of blue and red towers 
depends on the triangulation of the base strip. In addition, the dual graph contains "horizontal links" , 
which start on a vertical line of the same colour and end on the next (to the right, say) vertical line 
of the same colour. The position of the horizontal links encodes the three-dimensional triangulation 
of the original sandwich geometries. 

Looking at the subgraphs of one colour in the dual graph example of Fig. 3.4, one notes that they 
are particular examples of dual graphs for (l+l)-dimensional CDT spacetimes. It is indeed the case 
that the blue graph is precisely the graph dual to the triangulation of the initial two-dimensional 
surface of the sandwich geometry, and the red graph the dual to its final surface. We may therefore 
think of the bi-coloured graph as a particular superposition of these two uni-coloured graphs. Bi- 
coloured graphs of this type are again of the form of heaps of pieces in the sense of [90] , which means 
that we can map them to dimer configurations as in the (l+l)-dimensional situation (c.f. appendix 
A-II), but with the difference that now the one-dimensional dimer model will also be bi-coloured, with 
a fixed sequence of blue and red sites and with the blue (red) dimers linking a blue (red) site with the 
next one of the same colour, as depicted in Fig. 3.4. 





Figure 3.4: From sandwich geometries with product structure to hard dimers: (a) vertical hnes representing towers 
over a triangulated base strip; (b) horizontal links encode how the tetrahedra are connected, the corresponding weights 
can be read off as indicated; (c) projection to a sequence of hard-dimer configurations, as explained in appendix A-II. 
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The upshot of our considerations so far is that for a sandwich geometry associated with a fixed 
base strip S an analogue of formula (1.32) holds, namely, 

where 

^Sfi,n,iu,v,w) = J2 (3.12) 

coloured hard— dimcr config. D^^ 

is the partition function of the coloured dimer problem for a given sequence S of i?i blue and i?2 red 
sites, \D\f) {\D\r) is the number of blue (red) dimers in the configuration D, and \r\D\ is the number 
of crossings between dimers and sites of different colour. We should emphasize that our application 
of the inversion formula relates only to the second sum in (3.8), since we sum over the triangulation 
of the towers for a fixed triangulation of the base strip. This allows us to rewrite (3.10) as 

whose further evaluation will be the subject of the remainder of this chapter. 



3.4 Random matrix formulation 

The coloured-dimcr partition function can be written in terms of transfer matrices, by introducing 
three vectors (1,0,0), (0,1,0) and (0,0,1) corresponding to an empty link, a blue link and a red link 
between neighbouring sites. A transfer matrix can be associated with transitions between these states, 
taking into account that hard dimers are not allowed to touch each other, and that a dimer of a given 
colour connects only sites of the same colour which are first neighbours (among the sites with same 
colour), and may cross sites of the opposite colour in between. The transfer matrices can therefore 
be associated with the sites or vertices of the dimer model, and their explicit form depends on the 
colour label of the site. The transition empty-empty gets a weight 1, the transitions dimer-empty or 
empty-dimer get a weight or y^, depending on the colour, and the crossing of a site by a dimer 
of any colour gets a weight w; all other possibilities have weight 0. In this way we can associate a 
matrix 



with every blue site, and a matrix 
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(3.14) 



(3.15) 



with every red site. Consequently, the partition function Z^'^'^ ^ {u,v,w) can be expressed as a 
product of Ri matrices A and i?2 matrices B, ordered according to the sequence specified by S^n^^n^. 
(Note that the matrices A and B do not commute.) We can now write the partition function (3.11) 
associated with a fixed sequence S^n^^n^ as ^ 



^Our choice of taking the trace in the denominator imphes periodic boundary conditions in one of the spatial 
directions, by gluing together the first and last prism tower of the sandwich. Since we are leaving the other spatial 
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where the matrices A and B are the same as A and B, but with the substitution y/u i\/u, y/v iy/v. 
In order to express the sum over sequences S in eq. (3.10), it is convenient to combine the two matrices 
and define a general transfer matrix 

^ 1 qji^/u (1 - qj)iy/v\ 

Mg. {u, V, w) = qjA +{l-qj)B= qji^ (1 - qj)w , (3.17) 

V(l-9j)i\A^ qjW J 

with qj taking the values 1 (M = ^) or (M = B). The sum over sequences S is then replaced 
by a sum over sequences of elements {0, 1}. In this formulation, the calculation of the partition 
function assumes the form of a problem of products of random matrices [91], and we can rewrite the 
full partition function (3.10) for a fixed number N of prism towers (equivalently, a total number of 
matrices A, B in the sequence) as 

This in turn can be thought of as the average of (3.16) over all possible configurations {9j}jv (with 
any R), with qj = 1 and qj = each having probability p= 1/2. Introducing the notation 

N 

\PN,{<i}{u,v,w)\=Trl[Mg^{u,v,w) (3.19) 

we can finally write 

^"(-'-'"'^^=^^=^"< |P.,,,(U,^)| ^^-^.^-^-^- (^-20) 



3.5 The zeros of the denominators 

An important feature to be noted about the inversion formula is that it maps an infinite series with 
positive coefficients to the inverse of a finite sum with alternating sign coefficients. The latter will 
have a real root corresponding to the radius of convergence of the infinite series. This is evident in 
the (H-l)-dimensional case, formula (1.41). The same is clearly true for the (l+l+l)-dimensional 
case (3.11) too, which will have a two-dimensional locus of zeros in the three-dimensional parameter 
space spanned by u, v and w. However, the location of this locus will depend on the sequence {qj}, 
and we must keep in mind that we still have to sum over all the sequences {qj}- Consequently, 
we will be interested in determining the envelope of all the loci. In the absence of an analytical 
sohition, which appears difficult to come by, we will determine its location in the limit of infinite N 
by studying particular classes of sequences and by numerical computations, in conjunction with a bit 
of (well-motivated) conjecture. 



3.5.1 u = V 

Let us temporarily assume that x = y (and eventually = 1) in the Laplace transform of the one-step 
propagator, (3.7), and therefore u = v. This clearly makes things easier, but should still permit us 

direction open, the spatial slices of our model have the topology of a cylinder. Instead of compactifying one of the 
directions, we could also have left it open, in which case we would have to specify the boundary lengths as boundary 

conditions, or introduce conjugate boundary variables analogous to the x- and {/-variables in (1+1) dimensions, cqs. 
(1.39) and (1.41), and contract the product of matrices with the vectors (1 ixi 1x2) and (1 iy\ iy^)- This would 
complicate the treatment of the partition function considerably. 
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to extract information about the phase structure of the model, as is the case in both two and three 
dimensions (sec [57] and [58] for illustrations of the latter). This is not an implausible proposition in 
the sense that we want to find the critical line in the k-X plane, which concerns the "bulk" behaviour 
of the spacetime and should be insensitive to details of the boundary data, like those encoded in x 
and y. 

In order to determine the zeros of the polynomials \Pn{u,v = Ujw)] = O^i -^g^ ^ — 
u,w), we write this product as ^"i_B"2 ^"m-i^^m fQj. gome sequence of M < iV positive integers 
{ni, ...jHm} which sum up to N, '^j = ^- Note that from u = v follows B = JAJ, where the 
matrix 

/l o\ 

J= 1 (3.21) 

Vo 1 o) 

is a projector {i.e. = Id), so that we can write 

\PN^{gy\ =Tr{A'^'JA'^^J...A'^"-^JA'^"J). (3.22) 

Note furthermore that 







(a")i2 





A" = 


(a")i2 


(a")22 







I 








(3.23) 

^ u;"/ 

where a" is the n-th power of the matrix 




found in (1.34) above, and which has the following properties: 



= —^ r- = (3.25) 

(a")i2 = ix/^(a"-^)ii (3.26) 
(a")22 = (a")ii - (a"-')ii (3.27) 



where A± = (1 ± \/l — 4u)/2 are the eigenvalues of the matrix a {w is the third eigenvalue of the 
matrix A), and Un{x) is the n-th Chebyshev polynomial of the second kind. 

u=0 

For the subcase u = everything is trivial since |P;v(G,0, w)| = 1 for every sequence. (Wc ex- 
clude degenerate sequences with N matrices of one type and none of the other, which would yield 
|PAr(0, 0,w)| = 1 + . At any rate, they would give a negligible contribution in the average.) We 
thus have Zn{0,'w) = 1. 

w=0 

For the subcase w = the trace (3.22) reduces to 

M M ^ 

\pn,{,}\ = ^0 = n(""*)ii ^ n"'^t^"^(^) (3-28) 
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by virtue of (3.25), which imphes that (3.22) becomes the product of M uncoupled partition functions 
of the simple (uncoloured) hard-dimer problem. 

Note first of all that the term in (3.25) does not vanish at u = 0, because its product with 
Un{-2^) gives a polynomial of order [|^] in u whose zeroth order term is 1. Prom mathematical 
handbooks we know that /7„(x) has zeros only in the interval [—1,1], more precisely, at the values 
Xi = cos (^jqij-Ti'), for i= 1, n. It follows that no zeros can occur for m < 1/4. By contrast, for fixed 
u > 1/4 there is always a sequence, for some N, such that l-Pw.fq}! has zeros between 1/4 and u. For 
positive e = u — J, the occurrence of the first root is when — i= < a;i = cos(— Wtt). For small e 

this requires a large n. And this in turn means that (some but not necessarily all of) the sequences 
contributing to the critical behaviour at ty = and u = j possess a large group of consecutive matrices 
of the same kind; particular examples of this are configurations like AB^~^. A predominance of such 
configurations would signal a decoupling of neighbouring slices in the three-dimensional geometry, and 
a geometric degeneracy of the model there. 



u,W7^0 

For both u and w different from zero the general expression for (3.22) becomes highly non-trivial. 
To get a better idea of the distribution of zeros, we used the program Matlab to plot some of the 
roots of \Pn{u,w)\. This is in principle straightforward, but for increasing A'' requires considerable 
computing power, because the number of sequences grows like 2^. A superposition of several such 
plots is shown in Fig. 3.5a.^ From this picture it appears that the singularity-free region for u < 1/4 
found for u; = is also present for w > 0, up to some value of about 0.5, and then starts to shrink 
toward zero. Because of the limit on N in our numerical study, we could not determine whether the 
value u = is reached at w = 1 or rather some small u = e> 0. 

What we can do analytically to put more stringent bounds on the envelope is to study specific 
sequences and find their locus of zeros asymptotically for ^ oo. We find that {u,w) ~ (-jjw) 
continues to be an accumulation point for the zeros of the sequences AB^^^ for all values of w < 1. 
For sequences of the form A^B^ it can easily be shown analytically that a solution o{Tr{A^ B^) = 
for A?" — > oo is given by w = X+{u) = (1 -|- ^/T^^4u)/2. For the alternating sequences [AB)^/'^ we 
find instead a locus given by uw = 4/27, which is tangent to the curve w = X+{u) in the point 
{u,w) = (2/9,2/3), and otherwise lies completely to the right of it (in a u-w-plot like that of Fig. 
3.5a). Many other sequences with regular distribution patterns for A- and matrices whose large- A'' 
limit we have studied seem to have an asymptotic locus of zeros which is a smooth curve lying in 
between w = X^{u) and w = 4/(27u), and passing through their common point (2/9, 2/3). 

Based on these findings, and in the absence of any evidence to the contrary, we conjecture that the 
(m, u')-rcgion free of singularities is defined by the simultaneous conditions u < 1/4 and w < A+(m), 
whose boundary is given by the red curve in Fig. 3.5a. To its right we have plotted the locations 
of zeros, both of random sequences corresponding to some gi's, and of special sequences we have 
been able to treat exactly. Interestingly, this singular curve can be characterized in terms of the 
eigenvalues (A-|-, A_, w) of the matrix A, as the curve along which the largest two of the eigenvalues 
are degenerate. For u> 1/4, the X± are distinct and complex, and thus cannot be equal to w (which 
is real). For u = 1/4, A± are both equal to 1/2 and larger than w, as long as w < 1/2. For w > 1/2, 
the two largest eigenvalues are exactly defined by w being equal to the larger one of the pair A±. This 



Typical A''-values for the sequences ranged between 50 and 200, and the total number of sequences analyzed was on 
the order of 50. 
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Figure 3.5: (a): The critical line, combined from u = 1/4, w = X+{u), together with plots of the zeros of |P]v,{5}| for 
random sequences {q} and for sequences AB^~^ and A~ B~ at different values of N, in the u = v plane, (b): The 
plot of the critical surface in the full parameter space of the (u, v, w). 

Following this line of argument further, we can give yet another characterization of the singular 
line. It can be cheeked that the trace \Pn\ is always real. Denoting the eigenvalues of Pn by Ai, 
i = 1,2,3, their sum must therefore be real. Their product is the determinant of P^r, which from 
(3.22) is easily computed as the product of the determinants, A1A2A3 = X'^X^w'^ = {uw)^ , which is 
always real and positive for u,w > 0. This leaves only three possibilities for the signature/character 

of the A^: (+ + +), (H ) or (+ c c), where c and c denotes complex conjugates. Of course with 

signature (+ + +) the trace can never be zero, and this is precisely the case in the region u < 1/4 

and w < A+(m). From this we cannot go directly to signature (H ) because two of the eigenvalues 

would have to pass through zero, in which case the determinant would become zero, leading to a 
contradiction. The transition occurring at the singular line should be to the region with (+ c c). 
From there, once two of the eigenvalues are complex, their real part can become negative and also the 
trace can become zero. This scenario is confirmed by numerical computations, but a general, algebraic 
proof is at this stage still missing. 

3.5.2 u ^ V 

Using the same mixture of analytical and numerical methods as for the case u — v, and with a com- 
parable computational effort'' we have found a similar picture in the full three-dimensional parameter 
space with u ^ v. The critical surface is determined by the same relations as for u = v above, but 

^amongst other things, by systematically "scanning" the three-dimensional coupling space in {v, ui)-planes for various 
fixed values of u 
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with u substituted by max(M, v). That is, for u < v the critical surface is given by w = 1/4 for w < 1/2 
and w = X+{v) for w > 1/2, while for u > u it is given by u = 1/4 for w < 1/2 and w = A+(m) for 
w > 1/2. Since our extensive numerical and analytical checks have turned up no contradictions to 
this picture, we will in the following assume it to be correct. The critical surface is depicted in Fig. 
3.5b. 



3.6 The large-iV limit 

We will now assume that we are inside the region free of singularities. In this region the partition 

sum is convergent and its large-iV limit well defined. The boundary of this region is the critical line 
(or critical surface for u ^ v). Wc want to characterize the (non-)analytic character of the partition 
function Zjs[{u,v^w^lS.t = 1) on this boundary after the limit for large N has been taken (see the 
discussion in the previous section). Initially the evaluation of this function seems an insurmountable 
task, because according to (3.13) and (3.18) we have to sum over the inverses of the matrix traces 
resulting from the application of the inversion formula. However, it was exactly with this difficulty 
in mind that wc introduced the reformulation of the partition function in terms of random matrix 
products in Sec. 3.4 above. Techniques available for such products of random matrices will help us to 
estimate precisely the limit we are interested in, namely 

L-\{u,v,w) = lim In / , „ . r-r\=— ln2-|- lim —\Q.ZM{u,v,w,At=l), (3.29) 

^ ' ' ^ N^oo N \\Pn{u,v,'w)\ / N^oo N ^^v ' ' ' n \ j 

where, as in (3.20) above, the average has been taken over all random sequences {qj}. (The reason 
for the notation L_i will become clear soon.) The aim of this section is to compute L-\{u,v,w), at 
least in a perturbative expansion around the critical point (which still needs to be identified). 

3.6.1 Generalized Lyapunov exponents 

For a given sequence {qj}, define a real number 7 by 

7=^ln|P;v|. (3.30) 

Since {q^} is random, 7 will be a random variable as well; its probability distribution has, for large 
N , the general form [91] 

7r;v(7) oc e^^(^). (3.31) 

The function 5(7) is called large deviations function; it is a convex function that has a maximum in 
7 = 7, where 7 is defined as 

7= lim ^ln|Piv| (3.32) 

JV— >(x) iV 

and called the maximum Lyapunov characteristic exponent.^ We will choose 5(7) = 0, so that the 
distribution i^Nil) is normalized at the leading order for N ^ 00: 

yd^e^^(^)~e^^(^) = l. (3.33) 



*In general, the Purstenberg theorem [91] guaxcintees that 7 exists with probabihty 1 and is a non-random quantity, 
i.e. 7 = lim^v^oo W^^^ I^JvD- This means that 7 is a self- averaging quantity, and that §(7) must be a function peaked 
around 7. Note that instead |PAr| = e-'^'"' is not in general a self-averaging quantity: it will be shown later in this section 
that (\Pn\) = e^^i and {|Pivp> = where the L„ are 0(1), so that {(|Piv| - (|Pjv|»^> 

order as (IPatI)^, as long as L2 > 2Li. 
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Figure 3.6: Determining the function S{^) numerically for u = v = 0.24 and w = 0.1. The Gaussian 
approximation is shown as a dashed line. The blue continuous line is a fourth-order fit. The straight 
(dashed-dotted) line on the left has slope 1, and is tangent to 5'(7) in 'yli- In this case, the Gaussian 
approximation is not sufficient to compute L_i from (3.36), because in the point 7 — j'Li, S{j) is very 
different from the Gaussian fit. The figure illustrates a "worst-case scenario" in that for the chosen 
parameter values the deviation of the actual curve from the Gaussian approximation is maximal, see 
also Fig. 3.7. See appendix A-III for details on the numerical method. 



We can then write 



|Pjv(w, u, w)| 

where 7* is the solution to the saddle point equation^ 



07 



Similarly, one can define the Legendre transform of S{'^) (or generalized Lyapunov exponent) by^° 

Ln = lim 1 HlP^n = lim 1 ln(e"^^> = ^(7:) + n^:, (3.36) 

N^oQ iV N—*oo iv 

where 7* is the solution of 

^ = -n. (3.37) 

a-y 



^We are assuming that S{-/) is analytic and convex, so that the solution to the saddle point equation exists and is 
unique. Even if we cannot prove this assumption, it is strongly supported by numerical simulations, see Fig. 3.6. The 
same assumption guarantees that L„ is an analytic function of n and is at the basis of the replica method used in the 
next subsection. 

'^'^Obviously if 7rjv{7) = ^(7 — 7), then Ln = n^. 
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Clearly, for integer n, e^-^" = (|Pjv|") is the n-th moment of |Pjv|- If S{j) is Gaussian, 5(7) = 
— ^^2^1^ , we easily find 

7; = 7 + na\ Ln = n^+ ^n^a^. (3.38) 

In this case the knowledge of the first two moments Li and L2 is equivalent to the knowledge of 7 
and (T, and determines the full curve _L„, now regarded as an analytic function for all values of n. If 
5(7) is not exactly Gaussian, a systematic expansion of L„ in powers of n, 

^" = E I""' (3.39) 
fc=i 

can be obtained starting from the expansion of 5'(7) in powers of 7 — 7 and solving eq. (3.37) order 
by order, given that n = 0(7 — 7). The knowledge of the first k integer moments Lfe is equivalent 
to the knowledge of the first k derivatives of 5(7) in 7 = 7 and allows us to reconstruct Z/„ up 
to a given order of approximation. Therefore, one possible way of investigation is to measure S{'y) 
numerically, use a Gaussian fit to estimate 7 and cr^, and eventually compute the corrections coming 
from the cubic, quartic, ... terms of 15(7). An illustrative numerical result for the function 5(7) is 
plotted in Fig. 3.6 for u = t; = 0.24, w = 0.1, where the deviation from Gaussianity is maximal. The 
Gaussian approximation becomes better as the critical point is approached. The insights gleaned from 
this numerical analysis become most powerful when combined with a different tool for computing the 
moments Ln, based on the so-called replica trick, to which we will turn in the next subsection. 



3.6.2 Replica trick 

A more efficient strategy is to compute the integer moments Ln directly with the help of the replica 
trick. For n a positive integer and for any sequence of matrices Mj, j = 1, - ■■ , A'', it is easy to show 
that 

TrH^^J =Trn<", (3-40) 
\ i ) j 

where (8) is the tensor product. Then, if the Mj are independently and identically distributed, 

(IPivD = ( I Tr n ) = Tr ( [] ) = Tr !]«") = Tr (M«")^ ~ (3.41) 



where is the largest eigenvalue of the matrix (M®"), which can be easily evaluated for n small. 
Thus one has, for positive integer n, 

Ln = \nvn. (3.42) 

Knowing the function I/„ for integers < n < A; allows one to compute the first k coefficients Zi, • • • ,lk 

of the n-expansion (3.39) simply by solving a linear system. This yields an approximate expression for 
7 = lim„^o — h a-iid for L_i, which is the quantity we want to compute. Obviously the method 
assumes that L„ is an analytic function of n and would break down if there were a singularity at some 
value of n. 
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Figure 3.7: The numerical result for L_i as a function of u{— v) for w ~ 0.1, 0.5, Wc, 0.9. For 
w > Wc — 2/3, the simplest approximation L_i = — Li (full line) works well. For w = 0.1 and 0.5 
the Gaussian approximation, L^i = L2 — 3Li (continuous line), is needed to approximate the data 
better. For small w and u — > 1/4 it is clear that additional higher-order corrections are necessary, see 
also Fig. 3.6. 



First moment 

The simplest approximation is L„ = lin. The coefficient li is then equal to Li and is simply the 
logarithm of the largest eigenvalue of the matrix (M) = {A + B)/2, that is. 



(3.43) 



i^i{u,v,w) = i (2 + 'w + V(2 - w)2 - 4(u + v)^ 
Li{u, v,w) = In h'i{u, v,w). 
This approximation is equivalent to neglecting the fluctuations of 7 since 

e^^" = (IPatD = e"^^^ = {\PN\r. (3.44) 
The function i>i{u, v, w) has a singularity at 

2 — w 



u + v=^-^ j , (3.45) 

which intersects the critical surface only at the critical point {uc,Vc,Wc) = (2/9,2/9,2/3). This point 
is therefore a good candidate for a singularity of £-1, if it does not cancel at the next order. 

Second moment: Gaussian approximation 

Including the first correction is equivalent to the Gaussian approximation and is obtained by consid- 
ering Ln — hn + ■^hn^ as in (3.38) above. One has 

= 7 - 2ii - ^Lz, h^cj^ ^L2-2Li, (3.46) 
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with Li — In 1^1 as before and L2 = In 1/2, the logarithm of the largest eigenvalue of the 9x9 
matrix (M®^) = (A®^ + S®^)/2. For our searched-for expression this implies L_i = L2 - 3Li. 
We have calculated the eigenvalue 1^2 with the help of Mathematica for u = v, but refrain from 
reporting it here, because it is long and not particularly illuminating. However, the expressions 
simpliiy considerably on the critical line u = v = w — w^, where one has 

to > - 

iy2(w-w'\w) = { r , ^' (3.47) 

' ^ 2-2w + w^ + V4(1-w)2+«;2(2-3m;)2 — ^ 



4 



w< 3. 



What is rather remarkable about this result is that for w > 2/3 we have 1^2 = i^i, i.e. L2 — 2Li. It 
means that the approximation L„ — nLi is exact^^ on the critical line for w > 2/3. On this line we 
then obtain L-i = — Inz/i = — Intw. 

More generally, we can write i_i = —Li + {L2 — 2Li) = —L\ + AZ/2 and study the behaviour of 

AL2 close to the critical line. For later use in our scaling analysis, we consider two linearly independent 
types of perturbation away from a given point on the critical line, labelled by its «;-coordinatc wq. 
The first one is for w = wq — 5 and u = v = w — uP' , thus describing motion along the critical line, for 
which we find 

1^2 



AL2 



for Wo > 2/3, (3.48) 

«i<52 + 0{5^) for Wo = 2/3, 

K2 + Ki5 + 0{6'^) for Wo < 2/3. 



The second one moves away from the critical line inside the diagonal plane u = v, and toward the 
interior of the singularity-free region according to u = Wo — Wq — e, yielding 

{Kie'^ + 0(e') for Wo > 2/3, 

Kge + Kge^/' + 0{^) for wo = 2/3, (3.49) 

K7 + K8e + 0(e^) for Wo < 2/3 

where the Kj in (3.48) and (3.49) are numerical constants. Interestingly, these corrections do not 
change the divergent term in the first derivatives of L_i with respect to u and v at the critical point, 
which govern the presence or otherwise of an infinite-volume limit. Neither, since V2 according to 
(3.47) has no further singularities, can singularities appear at second order at points with wo 7^ 2/3. 
In summary, since all correction terms give only finite contributions to the derivatives, it is clear that 
our first-order analysis already correctly identified (Uci'cW'c) = (2/9,2/9,2/3) as the only point at 
which an infinite- volume limit exists. 

Since AL2 is proportional to e at the critical point, we expect it to contribute to the calculation of 
the Hamiltonian. This would at first seem to necessitate an analytic determination for the eigenvalues 
of (A®^ -|- B®^)/2 for u ^ V, which is currently out of reach. Fortunately, all that is needed is a 
perturbative evaluation after inserting an ansatz for the scaling of u, v and w near the critical point, 
which is a pcrfcHitly fcvisililc task wt^ will perform in Sec. 9 below. 



-"^^To prove this statement one should prove that Vn = ^^"^ ^- Although this ought to be possible, we have 
limited ourselves here to checking, by sampling random values on the critical line above w > 2/3, that it is true for n 
large but finite. 
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Higher moments: beyond the Gaussian approximation 

The precision in determining L_i can be improved furtlicr by computing L3, L4, etc. In Fig. 3.7 
we report our numerical findings for L_i (see appendix A-III), together with the first- and second- 
moment approximation just described. For w > 2/3 the first-order approximation is excellent, and the 
numerical deviations from AL2 for w > 2/3 and any value of u are very small. Around w = 1/2, 
the Gaussian correction starts to be observable. For w < 1/2 and u ~ 1/4, the correction is large, 
c.f. Fig. 3.6, and higher-order corrections have to be taken into account to correctly reproduce the 
numerical result. 

For approaches to the critical point, we have found that ALr, the incremental correction for the 
calculation of L-i which results from adding the r'th one to that of the first r — 1 moments, scales like 
e*"/^ and thus is not expected to contribute to the calculation of the Hamiltonian. We have confirmed 
this by perturbative calculations for L3. 



3.7 Continuum limit, canonical scaling and properties of the 
slices 

As usual in dynamically triangulated models of quantum gravity, our next step will be to search for 
a continuum limit in which the details of the discretization procedure will be "washed out" , and only 

universal, physical properties will remain. A necessary part of this limit is to take our short-distance 
cutoff, the edge length a, to zero, while keeping the physical spacetime volume V finite. This is only 
possible if the number of building blocks in the simplicial manifold (the discrete spacetime volume 
V) simultaneously goes to infinity. More precisely, we have to take the continuum limit following a 
trajectory (parameterized by a) in the coupling-constant space which ends up (for a = 0) in a point 
where the expectation value of the number of building blocks diverges, and which approaches this 
point in such a way that the physical volume stays finite. In other words, we have to renormalize the 
bare cosmological (and possibly other coupling) constants so that relevant physical quantities remain 
finite. 

Our first step, the identification of a critical point with suitable properties, follows from the 
discussion of the previous section: the only point where the derivatives of the partition function-'^^ 
diverge is {uc,Vc,Wc) — (2/9,2/9,2/3). As a second step we will assume a canonical scaling, where 
the critical point is approached according to the canonical dimensions of the corresponding continuum 
coupling constants. To lowest non-trivial order in a, these are by definition 

ko±kc + ^, \nx ^ -Xb,c- a^X, \ny o± -Xb^c - a^Y, A ~ + a^A (3.50) 

for the inverse Newton, the two boundary cosmological and the bulk cosmological constants. In 
addition, a canonical scaling would require T = at and M = aN for the time and spatial extensions. 

Translated into the variables u, v and w (which were defined in the discussion preceding (3.10)), 
the canonical scaling becomes 

U=- g2aci/G-2o=X-2a='6iA y = 1 g2aci/G-2a''Y-2a^biA yj ^ g-acj/G-a^bjA 

9 '9 '3 ■ y ■ ) 

Already without detailed inspection we can anticipate difficulties with the standard interpretation of 
the bulk and boundary cosmological constants as the couplings conjugate to the bulk volume and 



^'^ Remember that the average volume is given by the derivative of the logarithm of the partition function with respect 
to the bare cosmological constant, which translates into a linear combination of u-^, v-^ and w-^, see (3.53). 
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boundary areas. Since k scales with the lowest power of a when approaching the critical point, terms 
proportional to 1/G will be the ones which survive in the lowest-order expressions of the continuum 
limit, unless unexpected cancellations occur. 

We can check this immediately by looking at what this particular choice of scaling implies for 
geometrical quantititcs like the volume of sandwich geometries. As argued in the previous section, the 
second moment in the replica method only contributes to the next-to-leading order, which therefore 
we will need for the Hamiltonian, but not to recover the continuum expression for the volume. This 
means that as a starting point for our computation we can take the simple expression 

Zn ~ e"^^" (3.52) 

as the partition function of a sandwich geometry. Assuming canonical scaling for the three-volume 
and using previous definitions, we have 

(V) = lima3(y), 

a— >0 



which after inserting the scalings (3.51) and taking the limit a ^ yields 

4V(c2 - 2ci)/G 

to lowest order in a. As expected, only the renormalized Newton constant G appears in the continuum 
expression for the volume. In addition, this expression indicates an anomalous scaling of either the 
time variable t or the variable N measuring the linear spatial extension. For example, if we insist on 
the canonical scaling Af = aN for the finite continuum counterpart of N, we find that the volume goes 
to zero like a^/^ instead of a as would be required for a canonical scaling T = at oi the time variable. 

If instead we set fc = identically, corresponding to working with a bare action which only contains 
a cosmological term the sandwich volume scales canonically according to 

(V) = a N — = oAf — , (3.55) 

Since we do not want to touch the canonical scaling of the three-volume, which would affect the 
interpretation of the model in a fundamental way, possible solutions are to either assume a non- 
canonical scaling only for Newton's constant, for example, fc ~ fee + or to insist on canonical 
scaling but choose the constants Cj such that C2 = 2ci and the order-a terms cancel each other. In 
the first case we would obtain 

A.^X + y - (2ci - C2)/G'2 

and in the second 

(V) = a'N ^ '^^ + ^^ , (3.57) 

which are the same, up to a finite rescaling of G. The scaling now is canonical, and the sandwich 
volume is governed by both the boundary cosmological constants and the Newton constant. We will 
see below that related issues appear when we try to derive the continuum Hamiltonian. 

Another quantity one can consider is the total, integrated scalar curvature Tl-tot = / d'^x^/g TZ. Its 
counterpart at the discrete level is exactly the term Rtot multiplying k in the gravitational action. 
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for which we have used the standard Regge prescription. We can express the average total curvature 
as a derivative with respect to k of the partition function, and then take the continuum Umit. The 
analogues of relations (3.53) are given by 

CR-tot) = lim a{Rtot), 

a— »0 

, dlnZN f f d d\ (3-5^) 

(i?*.) = = -csiV + [c, (2u- + 2v-y c,n.- j InZ.. 

Since our spacetime sandwiches in the continuum limit have only infinitesimal thickness and thus 
infinitesimal volume, it is more appropriate to work with the average curvature per volume, 

^=i^. (3.59) 

It is a well-known feature of dynamically triangulated models of gravity that this quantity generically 
diverges like a~^, unless specific cancellations occur. Wc find here not only the same behaviour, but 
even the same algebraic expression to leading order, namely, 

n=^( gl^) , (3.60) 
\2bi+b2j ^ ' 

regardless of whether we assume the scaling fc = 0, the canonical k '^kc^ ajG, or the non-canonical 
fc ~ fee + {a/GY. However, in line with our earlier remarks, if one makes the choice^^ 2ci = C2, 
a diff'erent, less divergent, behaviour is found. Combining this condition and canonical scaling, we 

obtain ^ 

' 3cl/G , fc2-4c,\ l^^^__3 ^3_g^^ 




461 + 262 1261-^62/ V 4 



For = we find instead 



We will encounter somewhat similar scaling relations in the discussion of the quantum Hamiltonian 
in Sec. 3.9 below. 



3.8 Transfer matrix for areas 

One way to construct the propagator for finite times t of a statistical model is by iteration of the 
transfer matrix T, which in our case would take the form {g2\T*\gi)- Although the determination of 
the complete transfer matrix for our three-dimensional model is out of reach, we have already argued 
in Sec. 3.2.2 that most of the detailed dependence of the propagator on the boundary geometries gi will 
be dynamically irrelevant, because of the absence of local degrees of freedom in the three-dimensional 
gravity theory. Following [57], keeping track only of the boundary areas Ai in the one-step propagator 
as we have been doing throughout this work, should be enough^* to determine the finite-time behaviour 
of the propagator in the continuum limit. Introducing the "area states" 

i.)^^i:m, (3.3, 



^^This can be achieved by fixing the finite relative factor r between the time- and space-Hke (squared) edge lengths 
of the triangulations (c.f. appendix A-I). Using eq. (A-7), it corresponds to taking r ~ 0.724. 
^^possibly up to an additional dependence on global Teichmuller parameters 
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where J^{A) is the number of triangulations of a given area A and g\A any triangulation with a given 
total number of triangles A, we obtain that 

Z{x,y,At=l)= x^^y^- ^ = 

^'^-^'^^ (3.64) 

= J2 x^^y^^A2\f\Ai)^U{A,)J^{A2). 

Ai,A2 

Assuming now in line with our earlier reasoning (see also [57]) that 

{A2\f\g\A,) ~ {A2\f\g\^J ^ for ^2,^1^00, (3.65) 
for all pairs g, g' of boundary geometries of the same area, the completeness relation 

dA\A){A\=l (3.66) 

holds in the largc-arca limit and we can use (A2|r|Ai) as our transfer matrix. One remaining problem 
is the appearance of the square-root term M {At)N {AtJ^x) in (3.64), which we will have to deal with 
either before or after performing the inverse Laplace transform of Z{x,y,/S.t = 1). Given that the 
number of triangulations of a given area A scales like 

N{A) ~ A-^e^°^, (3.67) 

with 7 = 1/2 as we will show at the end of this section, one possibility is to apply to Z{x, y, At = 1) 
a fractional derivative operator 

^.-^,-M^.|W=(^)^''(^)^'^(.,,,A. = l). (3.68) 



There are several definitions of fractional derivatives in the literature; what we need here is an operator 
such that D^e^^ = A"e^^, so that applying D'^l^^lly term by term to (3.64) we get rid of the 
entropy factors (the remaining exponential term e^°^ is unproblematic since it only shifts the location 
of the critical point). We need yet another property for our fractional derivative, namely, a "chain 
rule" , since the final expression for Z{x, y, At = 1) will involve some function of x and y. In particular, 
since the final expression is expanded in powers of In a;, it would be nice if the derivative would act 
on powers with the simple rule D"x" cx An explicit representation of a fractional derivative 

operator with the desired properties exists and is reviewed in appendix A-IV. 

Unfortunately, at the relevant value 7 = 1/2, for some terms in the final expression for Z(x, y, At = 
1) the fractional derivatives suffer from convergence problems. These can be circumvented in a some- 
what ad-hoc fashion by integrating by parts or by introducing a regularization in the fractional in- 
tegrals, which should be removed only after the ordinary derivative has been applied. (As explained 
in appendix A-IV, the fractional derivative is actually defined as an ordinary derivative acting on a 
fractional integral.) 

A cleaner alternative to get the correct Hamiltonian is to first work out the inverse Laplace trans- 
form and then identify the contribution coming from the entropy factor and remove it. This method 
on the other hand requires some special care for the following reason. Starting from (3.64), keep- 
ing only the two lowest orders in a (see (3.3)), and doing an inverse Laplace transform will produce 
N{Ai)6{Ai — A2) at the lowest order and some differential operator acting on 5{Ai — A2) at the next 
order. If we divided this operator by M{At) to obtain the Hamiltonian, we would make a mistake 
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since at next-to-leading order A2 ^ Ai + adA, which needs to be taken into account. To do so let us 
rewrite the coefficients of the power series (3.64) as follows: 



{A2\f\Ai),/Jf{A;W^2} = {A2\^M{A) f ^J^^{A) = 

1^1) = 



= {A2\Af{A)f\A^) + {A2\y/Af{A) 



(3.69) 



= J^{A2){A2\{1 - aH)\Ai) - a^/j7(A^{A2\ 



\A,) + 0{a% 



where we have introduced the area operator v4.|j4i) = If we define the auxiliary Hamiltonian 

(^2|T|^i)VWi)W2) =: Af{A2){A2\{l - aH')\A,) + 0{a^), (3.70) 



we can invert (3.70) by using (3.69) to find the real Hamiltonian 



H = H' 



N{A) 



(3.71) 



It turns out that also with this method convergence problems appear in some of the inverse Laplace 
transforms to be performed. Fortunately, since the divergences in the fractional derivative method and 
in the auxiliary Hamiltonian method appear in different terms, they can be regularized consistently 
by requiring the end result in both methods to agree. 

Wc conclude this section by showing that 7 = 1/2 in the scaling relation (3.67) for the number 
of boundary geometries. The entropy of this class of triangulations can be deduced from its leading 
critical behaviour. Summing over all the triangulations with a cosmological weight - as is appropriate 
for evaluating the gravitational partition function in two dimensions - one has 



r 



-\A 



^A-Te-(^-^«)^-(A-Ao)'^-^ 



(3.72) 



Of course, the boundaries of our model are nothing but standard (l-l-l)-dimensional causal dynamical 
triangulations with periodically identified time. For precisely this ensemble we have already established 
earlier the relation (1.34), with whose help we will now be able to find 7 analytically. Summing over 
t to include all triangulations (and including a factor of 1/2*, which is necessary for convergence), we 
simply get 



El '72d 



1 



1 



ln2 + lnA+(u) ^/l - 4u 



(3.73) 



since there is no power-like subleading behaviour in t. Comparing with (3.72), we deduce that for the 
(l+l)-dimensional boundaries we have 7 = 1/2. 



3.9 The Hamiltonian 

The continuum dynamics of the model is encoded in the Hamiltonian operator, which can be derived 
by taking the continuum limit of the transfer matrix, as we explained before. To obtain the area- 
to-area transfer matrix we will first remove the A'^-dependence from the partition function, to avoid 
that the propagator will be dressed with the exponential decay of the one-dimensional boundaries 
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correlation. To do so we normalize^^ Zjq and sum^^ over iV, yielding 

Ziu, w) = V er^{^-^iu,,v,,w.)-L.^{u,v,w)) ^ — i ^ (3_74N 

^ L_i{Uc,Vc,Wc) - L_i{u,v,w) 

(note that L-i{uc,Vc,'Wc) — L-i{u,v,w) > 0). To derive the quantum Hamiltonian in the "X- 
representation" , we substitute the known function L-i{u,v,w) as well as an ansatz for the scaling 
relations into the evolution equation 

7p{x, t + l) = j> ^Z{x, y-\ At = 1)^(2/, t) (3.75) 

for the wave function (remember that x and y were absorbed in u and v), and evaluate its continuum 
limit to order a, namely, 



(l-aHx + 0{a''))i;{X,T) = -a^ —Z{X,-Y;a)7p{Y,T). (3.76) 



In (3.76), ^ is chosen such that the integration contour lies to the right of the singularities of \j:{Y,T) 
and to the left of those of Z{X, —Y; a), around which we have to close it. 

To show that to lowest order the identity operator is reproduced and to extract the Hamiltonian 
at the next order, we clearly need the first two orders in a of L_i{u,v,w) — L_i{uc,Vc,'Wc)- As 
discussed in Sec. 3.6.2 above, to get the order-a^ terms right wo need to use the replica trick up to 
the second moment approximation, i.e. L_i = L2 — 3Li. It docs not matter that we do not have a 
closed analytic expression available, since we will only need the expansion up to order a^. This can be 
found perturbatively around the solution at the critical point (where we can solve for the eigenvalues 
of (^®^ + B®^)/2), with the eigenvalue problem perturbed in accordance with the chosen scaling of 
coupling constants. 

In anticipation of difficulties with the scaling of k we first discuss the results for A; = 0. Given the 
analytic expression of Li and after computing 

we have all we need to calculate the partition function, resulting in 

.^ZiX. -y; A) ^ ^ + (5 + ^ - ^y_^) + 0(„>). (3^78) 

where we have absorbed a finite numerical factor in A. This expression is in perfect agreement with 
our previous discussion, since due to the entropy factor M{A) ^ ^-i/2gAoA fj.Q]-Q ^]-^g boundaries we 
expect the leading term to have an inverse square-root singularity if the transfer matrix is to reduce 



The normalization is needed in order to take the continuum hmit of the partition function at fixed A/". An analogous 
situation is found in (1+1) dimensions (see, for example, [51]). There, if we start with the generating function (1.41) 
for the transfer matrix, we find to leading order twice the identity operator in the continuum limit. Consequently, if we 
iterate (1.41) t times, we will get a multiplicative factor 2*, which has to be removed if we want to obtain the continuum 
T-propagator (the analogue of our partition function at fixed A/"). 

'^^Since we are evaluating the continuum limit, it is perfectly legitimate to sum over N after having performed the 
large- A'^ limit. Performing the continuum limit first and then integrating over the continuous variable associated to A'^ 
would lead to the same result. What is more important is that wc arc not including any power-like corrections to the 
behaviour of Zjv, which could in principle occur. This is justified by the fact that it can be shown with the replica 
method that they are absent for the positive integer n analogues of Zj^. 
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to the identity at lowest order. Using the fractional derivative method before the integration in (3.76) 
or, alternatively, the method of the auxiliary Hamiltonian one finds 

. 3 3 , 1 a 1 1 , , 

for the quantum Hamiltonian in the "^-representation" (with A denoting the (finite) continuum 

arc^a). This looks like a bona- fide Hamiltonian, with a second-order kinetic term and a potential 
term depending on the cosmological constant. In the absence of a dimensionful Newton coupling G, 
the factor multiplying the kinetic term must be present for dimensional reasons. Note also the 
appearance of the factor 3/2 in front of the first-order derivative, ensuring the self-adjointness of the 
Hamiltonian with respect to the trivial d^-measure. By introducing a new variable 

L = aVA (3.80) 

and simultaneously defining new wave functions^^ 

m = ^i^i^h (3.81) 

the Hamiltonian becomes 

with the 1/L-term having disappeared from the potential. Apart from the cosmological term, which 
has the appropriate dimension for a (2+l)-dimensional Hamiltonian, this is exactly the quantum 
Hamiltonian of the (l-l-l)-dimensional CDT model [42] (see also [51] for a similar transformation of 
variables). 

Since by setting fc = we have not included any spacetime-curvature term in the discrete action, it 
is clear that the kinetic terms in (3.79) and (3.82) have their origin in the "entropy" of configurations, 
or, using a continuum language, in the non-trivial path integral measure underlying the dynamically 
triangulated model. This is also underlined by the absence of such a kinetic term from a related 
(l-|-2)-dimensional model considered and solved in [50], which uses the same product structure as our 
model, but works instead with a fixed, flat base manifold. Rather intriguingly, this gravity-inspired 
model (for finite t) can be related through an inversion formula to a problem of hard hexagons on a 
regular triangulation. Its one-step propagator can be seen to resurface in our model as contributing 
just one of the terms in the partition function, namely, .Zjv = l/Tr{AB)'^/^. Solving the model in 
the large- A'' limit does not require the replica trick, and simply leads to Zj^ <^ l/Amax, where Xmax 
is the largest eigenvalue of the matrix AB. It is straightforward to extract the Hamiltonian, which 
only contains a term proportional to the area, and no derivatives, which implies that the area is not 
a dynamical quantity in the continuum limit^*. 



^^With such a change of wave functions wc make sure that the measure is preserved, in the sense that 

j dA i>r(A)i>2(A) = j dL ML)ML)- 



Note that we also require 
which implies 



J dA M'AJH^M'A) = J dL ML)HlML), 
Hl = Ha-^l + VI \ha^l, 



with the A in the Hamiltonian substituted according to (3.80). 

^*The mapping to the hard hexagon model of [50] suggested a fractal dimension df = 12/5 for the (l+2)-dimensional 
siiuplicial complexes. The relation with our results for the associated one-step propagator is currently unclear. We 
conjecture that in our model, where we sum over all base triangulations, the fluctuations of the base serve as a stabilizer 
for the geometry and lead to an effective, non-anomalous dimension df = 3 for finite times. 
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Returning to the analysis of our model, the case with k ^ raises issues similar to those already 
encountered when computing the spacetime volume and curvature. The a/G-term dominates the 
continuum limit, giving rise to 



ZiX, -Y; G, A) = _J£!^ + ^,(6c.-5c. + 8c3)/G-V4X 4y-3(c./G)^ 
^ ' 3^(02 - 2ci)/G (36ci - 18c2)/G ^ 

(3.83) 

which is difficult to interpret within the scheme of the previous section (and to this order does not 
even contain any reference to the cosmological constant A!) Things simplify somewhat once we choose 
C2 = 2ci, in which case we obtain 

Z{X,-Y;G,A) = = + 0{a^). (3.84) 

4C3-C2 + V^-F-|(§) 

We have not bothered to include the next term in the expansion, because already the lowest-order 
expression gives nothing like the desired form J\f{A)S{A — A') upon performing an inverse Laplace 
transform. It is not totally inconceivable that the offensive terms in (3.84) could still be "argued 
away". Note that we have encountered the numerical term (4c3 — C2) previously in the expressions 
(3.61) and (3.62) for the average curvature. In particular, if we had a good reason for why (4c3 — C2) 
could be set to zero (which we currently do not), it would imply that in the case k = the average 
curvature was turned into a finite expression. However, even if (4c3 — C2) could be made to vanish, it 
would still leave us with the problem of having to absorb the shift in the location of the singularity 
from Y = X to Y = X — J (^^Y . In terms of the ^^-representation, the shift can be related to an 
overall multiplicative operator, which replaces the simple transfer matrix e""^-^ by e~"^^e 
This might be related to a redundancy among the renormalized couplings of the model inherent in 
the ansatz (3.50), (3.51), but we have not yet been able to make the relation precise. 



3.10 Going beyond the area representation 

Up to now we have (through the boundary cosmological constants) only taken the area information 
of the two-dimensional boundary geometries into account, but we expect from the canonical analysis 
for cylindrical topology also a Teichmiiller parameter as a dynamical degree of freedom (see, for 
instance, [95]). Although it does not seem possible to define the Teichmiiller parameter for a general 
triangulation exactly, a suitable discrete implementation for the triangulated cylinder should be the 
ratio between the area of the cylinder and the length of the cycle (the boundary of the cylinder) . The 
latter variable in our model is given by the number of matrices A ot B (for the incoming and outgoing 
cylinder respectively) appearing in the traces. The problem of keeping track of the Teichmiiller 
parameter is then reduced to that of keeping track of the number of matrices A ov B (as opposed to 
their sum N = A + B rs wc have been doing up to now). This task is in principle solvable with the 
help of the so-called microcanonical method [96]. 

The microcanonical method is based on a trick similar to the replica one. We consider the ensemble 
of sequences {qj} such that Np = pN variables are equal to 0, withp G [0, 1], i.e. there are Np matrices 
B and Ni-p matrices A. One can define generalized Lyapunov exponents in the microcanonical 
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ensemble according to 



Lnip) = lim ^IndPjvDpjv 

AT— >cx) iV 



1 ( rr n, \ (3-85) 



[pNj 

where the sum is only over sequences which contain pN matrices B. The microcanonical method is 
based on the identity 



{qj}pN \ J=i / {qj}pN i=i 

1 / , Tr (A®" + ^B®")^ 



x=0 



(3.86) 



where C is a contour in the complex plane containing the point z = 0. Calling iyn{z) the largest 
eigenvalue of the matrix (^®" + zB®^)/2, we have 

^ I Tr [] Mq. I ^—(fdz eJV[l"2-nW-pln^] ^ giV[ln2.„(C,)-plnG] ^ (337) 



where the saddle point (p is the solution of 

dlnvniz) 



dlnz 

Expanding the binomial by means of the Stirling formula we have 

I ^ I ~ (,-N[plnp+{l-p)lnil-p)] = ^Na{p) 
[pN 



p . (3.88) 



(3.89) 



and can finally write 

Ln{p) = ln2i/„(Cp) -plnCp - ct{p) . (3.90) 

This shows that the generalized Lyapunov exponents can be evaluated analytically also in the micro- 
canonical ensemble. 

We would now like to use this property to compute the partition function 

Ri,R2 N R<N 

f' dpe-^'^'-P^''-^-P''ZM,p^y2^e-^''' /''rfpe-(/^2-^i)iVgiv(L_i(p)+a(p)) ^ 

(3.91) 



where we have introduced the lengths Ri = N — R = (1 — p)N and R2 = R = pN of the boundaries 
of the incoming and outgoing cylinders, and the associated conjugate variables yiti and 112- We may 
again evaluate the integral by the saddle point method, namely, 

Z(/Xi, At2) - Q-t^iN^-Ni{t,,-^,)p-L_,ip)-a{p)+ln2) ^ ^g gj) 

N 



3.10. GOING BEYOND THE AREA REPRESENTATION 



91 



where p is the solution of the saddle point equation 

Like in the replica trick, we are going to infer L-i{p) from knowledge of the Ln{p)- If we take 
L-i{p) = —Li{p) as lowest-order approximation, we find (using (3.88) and (3.90)) that the saddle 
point equation reduces to 

Cp = e'^^-^^ (^^) • (3.94) 

Since we do not have an analytic expression for ^p, we will make an educated guess at its solution. If 
we replace L-i{p) by i„(p) in the above expressions, we obtain for the saddle point equation (using 
again (3.88) and (3.90)) 

= e-('*^-'*i\ (3.95) 

and for the partition function 

^ ^ ^ Mi-lni/„(e-(''=-''i)) ^ ^ 

We therefore conjecture that the solution of (3.94) is given by 



e 



-H2 



p = , (3.97) 

which reproduces the value (3.95) for Cp- With this conjecture we obtain 

^^"""'^^^ ~ m + ln4^i(e-(.-.0) + 2(M2 - ' M^^^^) ' ^^'^^^ 

If we choose the scaling 

yui ~ /Ltc + aMi , H2—iJ.c- aM2 , (3.99) 

together with the canonical scaling (3.50) for the other variables (with fc = 0, otherwise we run into 
the same problems as before), we find 

Z{X, Y, Ml, Ma) = , + 0{a^) (3.100) 

^ ^ Ml- M2 + VX -Y ^ ^ ^ ^ 

in the continuum limit. That this expression is compatible with the presence of the identity operator 
at lowest order can be seen as follows. Since we are now dealing with a new set of (orthonormal) 
states depending on two labels, 

with Af{A, R) denoting the number of triangulations of fixed area A and length R, we expect to lowest 
order in a that 

J2 (-"l^'-fl' \T\9\a,r} ~ S{A - A')5{R - R')M{A, R) . (3.102) 

3lA,R,g\A' ,R' 

Performing an inverse Laplace transform of (3.100), we find that 

-f 

AfiA,R)c,f^^e^°^ , (3.103) 
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where we have reinserted by hand the exponential term e'^"'^, which is otherwise absorbed in the 
critical value of the coupling constants, and which - as it should - satisfies (see (3.67)) 

dR U{A, R) = N{A) . (3.104) 

/o 

Work is in progress on including the next-order approximation, Z/_i(p) = Z/2(p) — 3Z/i(f)), and deter- 
mining the resulting continuum Hamiltonian. 



/ 

Jo 



3.11 Conclusions 

In this chapter, we have for the first time derived a continuum Hamiltonian from a three-dimensional 
quantum gravity model in terms of Causal Dynamical Triangulations. This was made possible by 
restricting ourselves to a subclass of all triangulations, which possess a product structure in the 
spatial direction, in addition to the usual product structure in the time direction imdcrlying the 
causal nature of the model. This restriction enabled us to apply a number of analytical methods, 
including the inversion formula, a random matrix formulation and the replica trick, to solve the one- 
step propagator with free boundary conditions and then perform its continuum limit. This led to an 
explicit expression for the quantum Hamiltonian (3.79) in >l-space, with wave functions depending on 
the area A of two-dimensional spatial slices. 

Taking into account the special, topological character of three-dimensional pure gravity, we have 
argued that for obtaining the full dynamics of the theory, it is sufficient to keep track of only a finite 
number of boundary data when deriving the model's transfer matrix from the one-step propagator, 
and sum over everything else. This should not be confused with an approach where the remaining 
degrees of freedom (essentially the conformal mode of the metric) are fixed at the outset. In our 
model, all of these are still present, but are summed over in the discrete path integral. For simplicity, 
what we have done in the present piece of work is to keep track of only a single variable, the volume 
(area) of the two-dimensional universe. Because of the cylindrical topology of our spatial slices, we 
expect there to be one additional Teichmiiller parameter, corresponding to the ratio between the area 
of the cylinder and the length of the cycle (the boundary of the cylinder). To include this parameter 
in our model requires to separately keep track of the number of matrices A and B (for the incoming 
and outgoing cylinder) in the traces, instead of only their sum N, as we have done in the main part 
of our work. This task is in principle solvable with the help of the so-called microcanonical method, 
as explained in the previous section. 

The fact that we have derived a non-trivial three-dimensional dynamics does a posteriori justify 
our restriction to a subclass of all three-dimensional Lorentzian triangulations. Since the continuum 
limit we have identified was obtained by setting the bare inverse Newton constant to zero, it is clear 
that the resulting continuum Hamiltonian describes the "collective" effect of the quantum fluctuations 
of we have summed over in the path integral, and which elsewhere we have called the entropy of the 
model. In particular, the net effect of these fluctuations is to generate a kinetic term for the area which 
makes the Hamiltonian bounded below, and therefore has the opposite sign from the corresponding 
term for the "global conformal mode" in the gravitational action. This further supports a mechanism 
already observed elsewhere in CDT models in three and four dimensions, namely, the presence of 
contributions from the path integral measure compensating for the divergence due to the conformal 
mode of the Euclidean gravitational action. 

There is a (small) price we are paying for using the additional constraint of a double product 
structure for the triangulations, which after all was not motivated by physics, but simply by the desire 
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to be able to apply a number of special solution techniques. The application of the inversion formula 
implies that we are evaluating first the limit of infinitely many building blocks in one of the spatial 
directions (the "height" of the towers), and only then the corresponding limit in the complementary 
direction, corresponding to what we have called the large- A'' limit. The effect of taking these two 
limits in sequence, and not simultaneously, is that the w,-w-diagram (Fig. 3.5), which describes the 
behaviour a/ier taking the first of the infinite sums, has only a single point at which an infinite- volume 
limit can be defined. This is in contrast to the corresponding phase diagram for three-dimensional 
CDT, which has an entire critical line Ac(fc) along which the infinite- volume limit can be taken. 

The analytic solution found here has also highlighted a problem with the scaling of Newton's 
constant, which has already been encountered in previous analyses of the phase structiirc of three- 
dimensional CDT models [59]. The issue is that a standard canonical, additive renormalization of 
Newton's constant as in (3.50) does not seem to lead to a good continuum limit, simply because 
the term containing the renormalized Newton constant G has the lowest order in a, and dominates 
everything else. The only way so far in which we have managed to obtain a nontrivial continuum 
limit was by setting fc 0, which implies that the Newton constant docs not play a dynamical role, 
which is of course compatible with the absence of local gravitational excitations in three dimensions. 
We believe a deeper understanding of the relation between certain scalings of the coupling constants 
and that of physically relevant quantities, like the curvature discussed in Sec. 3.7 above would be 
useful, and potentially relevant for a better understanding of the analogous issues in four dimensions. 
- Even if this and other issues still remain open, we think that the work presented here represents an 
important step in the analytical understanding of CDT models in d > 2 dimensions. 
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A-I The discrete action 

We derive in this appendix the exact expression for the action of our model. All the definitions and 
necessary ingredients were already introduced in [47], here we just recall and use them to derive the 
precise expression for our particular case. We start from the Einstein-Hilbert action plus a Gibbons- 
Hawking boundary, that is, 

S = l^ d^x^) (^A - +hj,M '^'^V^^(^)' (^-1) 

and then use Rcggc's prescription for the corresponding quantities on a simplicial manifold [22] to 
obtain the discrete action 

where (T„ is an n-dimensional simplex, Va^ its volume, M is the interior of the simplicial manifold 
and dM its boundary, 9„^^„^ is the dihedral angle of cts at a\, and the subscript b on the coupling 
constants stands for bare. 

In Dynamical Triangulations things simplify considerably because the edge lengths are held fixed. 
In the old, Euclidean models one would have only a single type of building block in any dimension, 
since all the edge lengths were chosen equal. In CDT the Lorentzian nature of the model allows us 
to have two different types of edge lengths, corresponding to time- and space-like directions. Because 
of the way the causal gluing rules are implemented, what used to be a space-like and what used to 
be a time-like link can still be distinguished after having mapped a Lorentzian triangulation to its 
Euclidean counterpart. Let us define the ratio r = where denotes the squared edge length 

for all time-like links, and l1 that for all space-like ones. Its allowed range in the Euclidean signature 
is r > e > 0, where e is a positive constant, which is the lower bound of the triangular inequalities 
and depends on the number of dimensions d (e = 1/4 for d = 2, e = 1/2 for d = 3 and e = 7/12 for 
d = A). An analytic continuation of r in the complex plane to negative values defines the (inverse) 
"Wick rotation" in CDT. Throughout this work, we will stick to positive values for e, corresponding 
to a Euclideanized, real partition function. 

In (2+1) dimensions we have precisely two different kinds of building blocks, tctrahcdra with three 
space-like and three time-like edges, and tetrahedra with two space-like and four time-like edges. We 
will use the notation Nij for the number of (i, j)-simplices of dimension i + j — 1 having i vertices in 
a constant-time slice t and j vertices in the subsequent one at time t+1. The Regge action (A-2) in 
this case becomes 

S =A,a' (^2,2^2,2 + V^3,i(^3,i + ^1,3)) - 7^ f 2^V^iv(f ^ + 7rVfivf/') + 

- [4^2,2^^2,2 + 3^i,3(iVi,3 + iV3,i)] + w{N2,Q + No,2)+ (A-3) 

- [2^^,2^2,2 + 3ei3iNl,3 + A^3,l)] ) , 

where we have distinguished the dihedral angles at time-like and space-like links by the superscripts 
t and s. Relevant volumes and angles are easily found to be 
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1 2r — 1 
cos 3 = - , cos G\^ = , 

^« 4r-3 1 ^ ^ 

cos 02,2 = , cos 02,2 = ■ 

Finally, using Euler's formula and the Dehn-Sonimcrvillc relations for triangulated manifolds with 
boundary we arrive at the expression 

S' = a(iVi3 + 7V3i) + /37V22+7^ (A-6) 

for the discrete action, where a, /3 and 7 depend on the dimcnsionless bare coupling constants A and 
k according to 



/ n ^ r 2r-13^ 1 \, V3r- 1 , , , , 

ct = —TxJr + i-Jr arccos ; — -tt + 3 arccos — k H — A = —c\k + oi A, 

' 4r-l 2 %/3V4r - 1 / 12 



/3 — { ~2-K\fr + 4\/r arccos V 2 arccos 7^ fc H -= — A = C2fc + h^X, "^^ 

V 4r - 1 4r - 1 y 6\/2 

7 = (6TT\/r — tt) k — c^k. 



A-II The inversion formula 

The proof of the inversion formula in (2+1) dimension (for a fixed sequence S^n-^n^ of i?i blue and 
i?2 red towers) proceeds very similarly to that in (1+1) dimensions given in [50]. We first switch to 
the dual picture of a triangulation, and assign a weight u or v to each horizontal edge, depending 
on its colour, and a weight w to each red-blue intersection, as depicted in Fig. 3.4 above. Next, we 




Figure 3.8: Decomposing the dual picture of a sandwich geometry into a sequence of projections onto its one- 
dimensional base line. The projection is obtained by letting the "bottom layer" of the heap of pieces, consisting of its 
lowest-lying edges, drop down to the horizontal line at the base. Each such projection defines a hard-dimer configuration 
on the horizontal line. 

decompose a given configuration into a sequence of projections onto the base, as illustrated in Fig. 3. 8. 
Such a projection gives rise to a hard-dimer configuration, if by a dimer we simply mean an edge 
linking two nearest vertices of the same colour. The "hardness" refers to the fact that any vertex of 

-"^^Note that the boundary of our "sandwiches" of three-dimensional spacetime (the product of a finite cylinder with 
an interval), which has space-like as well as time-like components, is connected and has the topology of a torus. 
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the base can be occupied by at most one dimer, incorporating a mutual avoidance of dimers. Such a 
construction is possible because the dual graph is a heap of pieces. One may think of the vertical lines 
in Fig. 3.8 as tracks along which the horizontal edges can slide up and down, with the only restriction 
of not being allowed to touch or pass each other. Using this decomposition, we can write the partition 
function Zxit) as 



hard dimcr config. C 



(A- 



where the sum extends over all hard-dimer configurations C on the one-dimensional lattice (including 
the empty configuration). The numbers \C\b and \C\r count the blue and red dimers in C, and 
I n C| the crossings between dimers and sites of different colour in the configuration C. For fixed C, 
Z'^'~^\u,v,w) is the restricted partition function involving those configurations having projection C, 
and from which we have factored out the weight y\^^''v^'-^^'^v^^^^ of the projected part. More generally, 
we have the relations 



E 

CDD 



\c\tjcw,„\nc\^(c} 



iu,v,w) 



(A-9) 



valid for any hard-dimer configuration D (eq. (A-8) corresponding to = 0). This expresses the 
fact that by completing any dual geometric configuration under consideration by a given row of 
horizontal edges (corresponding to a hard-dimer configuration D) , one builds each configuration having 
a projection containing D, i.e. having D as a sub-configuration exactly once, see Fig. 3. 9. Let us now 



D: 



)D: 



(a) 



(b) 



Figure 3.9: By completing an arbitrary triangulation of our sandwich geometry with a hard-dimer configuration D 
(a), we build a larger triangulation (b) whose projection C contains D. With this procedure (with fixed D) we build 
all triangulations whose projection C contains D exactly once. 



rewrite (A-9) as 



with 



g{D) = J2fiC)mD) 



g{D) = u 



c 



\DU„\D\r„\nD\ 



(A-10) 

(A-11) 
(A-12) 
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and 

f 1 if C D £>, 

aC,D) = { (A-13) 
I otherwise. 

Thinking of the space of configurations as a vector space in this manner, (A-9) can be thought of as a 
vector-matrix multiphcation, with an upper-triangular matrix ({C, D) and with all diagonal elements 
being equal to 1. Then Q has an inverse matrix /i(-D, C) of the form 

MAC)=^ (A-14) 
I otherwise. 

This property can be verified by noting that 

Y,aC,D)„{D,C') = {^ -^o^ (A-15) 
D I otherwise, 

and 

{1 if C = C", 

V * / 

We can now invert (A-9) as-^° 

JC|,^|C|.^|nc|^^Q ^^^^^^^^^ ^ ^ (A-17) 

Noting that Zs^^^^ {u, v, w) factors out of the sum on the right-hand side, we finally get 

(_y)|CU(_„)|C|.y,|nc|_^(C) (u,v,w) 
V, w) = ^ (_^)|z.|.(_^)|z.|.^|nD| ' (A-18) 

where we have used that \D\ = \D\i, + \D\r, leading to the minus sign in front of both u and v. Picking 
C = 0, we arrive at the fundamental inversion relation (eq. (3.11) above) 

{u, V, rn) = i!Lu,-v,wy ^^'"'^ 

where 

Z^^^,,^{u,v,w)= Yl t^l^l-'t^l^l'-u^l^^l (A-20) 

hard dimer config. D 

denotes the partition function for hard coloured dimers with fugacity u {v) per blue (red) dimer and 
weight w per crossing. 



A-III Numerical computations 

A direct numerical computation of the generalized Lyapunov exponents L„, defined in (3.36) above, 
is possible following the procedure given in [91]. One extracts a large number of values 7 at fixed N, 

■^•^This is the famous Mobius inversion formula of the theory of partially ordered sets, and ^J, is the associated Mobius 
function. 
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Figure 3.10: The numerical results for Ln{N), n = —1, 1, 2, 3, as a function of 1/N for u = 0.24 and 
w = 0.1. The full lines are fits to Ln{N) — i„ + An/N . The clashed lines are the analytical values 
for Li and L2 computed using the replica trick. They differ from the extrapolation by about 0.1%. 



computes the average (e"^^), and plots Ln{N) ^ N ^ ln(e"'*'^) as a function of 1/A^. Usually Ln{N) 
turns out to be linear in 1/N for large TV, 

L„(iV)~L„ + A./iV, (A-21) 

which allows for a good extrapolation of the data to = 00, see Fig. 3.10 for an example. The values 
of 7 are extracted as follows [91]: 

1. construct a random unimodular vector v (e.g., by extracting the components from a Gaussian 
distribution); 

2. extract the numbers qj and apply the matrix Mq^ to v to obtain a vector v' = Yij ^qj^j 

3. compute ~ e^''', which holds to leading order since the modulus of v' is dominated by the 
same largest eigenvalue of Ylj which also dominates Tr Y[j Mq^ . 

This will give us a set oiJ\f values 7^ for 7, from which we can construct a histogram 7r(7) and extract 
the function S{j) at given N, see Fig. 3.6. The coincidence of the curves for different N indicates 
that the asymptotic limit N 00 has been reached. To compute L„, we first have to compute the 
average of 7, 7 = 7ij ^^'^ then use 

By subtracting 7 we can compute the average as long as the quantity N{'-fi — j) ~ ^/N is not too large. 
In practice we are limited to iV < 50, but this is enough to get a very good linear extrapolation to 
A'^ —> 00. The difference between the numerical computation and the exact analytic expressions for Li 
and L2 is of the order of 10"'', see Fig. 3.10. We believe that a similar error affects the computation 
of L_i too. 



1 



«— 1 



(A-22) 
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A-IV Fractional derivatives and inverse Laplace transforms 

Although the notion of fractional calculus may at first appear somewhat extravagant, this is an old and 
well-studied subject in mathematics, with plenty of applications, as testified by the number of books on 
the subject (see, for instance, the list of references in [92, 93]). At the heart of the theory of fractional 
calculus is the definition and study of two operators J" : f{x) — * J^f{x) and : f{x) D'^f{x) 
on a suHiciently large class of functions {f{x)} and for a positive real number a, with the following 
properties. 

1. When a = n is a positive integer, the operator gives the same result as n-fold integration, 
and _D" gives the same result as the usual n-th derivative 

2. The operators of order a = are the identity operator. 

3. J" and D'^ are linear operators. 

4. For any a, (3, the semigroup property holds for J", namely, J"J^f{x) = J"~^^f{x). 

There are many inequivalent definitions of fractional derivative/integral satisfying these properties, 
which we are not going to review here. Instead, we will recall the definition and properties of the one 
that turns out to be relevant for our purposes, the so-called Weyl fractional derivative/integral. (Ref- 
erence [94] has a short review, which also contains a clear explanation for why the various definitions 
are inequivalent.) The Weyl fractional integral is defined as^^ 

JS,o.m = ^-j^ / dtfm - xr-\ (A-23) 

r(a) 

The natural definition of the derivative operator would be to take D'^ ^ = but this is divergent. 
A standard trick is then to define 

Dloofi^) - j]^}^'-''f{x), (A-24) 

where [a] is the integer part of a. The rationale behind this definition is that in this way is the 
left inverse of J^, i.e. D^J!^f{x) = f{x), the same as when a = n is an integer. With this definition 
it can be checked that 

Dlo^e-^^ = (-A)«e-^^ (A-25) 

which is exactly what we need, as explained in the main text. 
According to the analysis in the text, we need to compute 

Z{x, y, At=l)= Dlil^D\ll^Z{x, y, = 1) , (A-26) 

which in the continuum limit we can write as (see Sec. 3.9) 

Z{X, Y) = D^l%^^^D%^^yZ{X^ Y) = I d]IXd]/Xz{X, Y) . (A-27) 

At leading order in a we expect to have a S{Ai ~ A2) in the transfer matrix. Therefore we can 
take just a i-dcrivativc with respect to one of the two argmncnts. From the expression of Z{X, Y) 

■^'^The subscript "00" refers to the extremum of integration. Other choices of extremum would give rise to different 
(and inequivalent) definitions of fractional derivative/integral, which we would denote by Jx,xo > ^",xo ■ 
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given in (3.78) we have at leading order (l.o.) 

^D]llZi„{X,Y) = ^D]ll^^^- lim — [''°dt^=^= 
■ x,oc i-oA , ) ■ x.oo^^^ x,^^dXT{ll2) + (A-28) 

which is proportional to the Laplace transform of the delta function. Note that the limit here has to 
be taken after the derivative, otherwise it would be divergent. This is the kind of convergence problem 
mentioned in Sec. 3.8 (the problem would stay also if we would use two i-dcrivativc instead of one 
^-derivative). This is not really a big problem because what we need is an operator which gives the 
same result as (A-25) when acting on exponential functions. Since for exponential functions the limit 
can be taken both before and after the derivative, we have a freedom in choosing the order of these 
operations. 

For the cosmological term in (3.78), interchanging limit and derivative does not give any problem, 
and we get 



_a 1/4 1/4 A _ 2a 1 
, ^x,oo^y,oo ^ y)3/2 - ^ (X + y)2 ' 



which is proportional to the Laplace transform of aAAi6{Ai — ^2)- 

For the remaining term in (3.78) we have to choose again to take the limit after the derivative. 
What we obtain is 

«ni/4 ni/4 [5 XY \ _ yi^jx^-mxY + Y') 

i^X,oo^Y,oo\^Q + (^X + YyJ - 16(X + y )5/2 • ^"^-"^^^ 

(The constant term gives zero as a result.) The inverse Laplace transform of this term presents another 
divergence. This may be avoided by an integration by parts, but we prefer to take another route and 

use the second method presented in Sec. 3.8. 

The method consists in performing the inverse Laplace transform on Z{X, Y) directly, without 
use of fractional derivatives, and then recognizing the contribution of the entropy in the Hamiltonian. 
Let's see first how this works for the leading order term. We have to compute the integral 

f+ioo+fi . f+ioo+fj. 



YA2 ^ 

ioo+At J-ioo+fi 27ri ^/X -\-Y 



-ioo+n ^TTl J_ 

-ioo + n 2m sJZ 



dY / d^^zA.2^ (A-31) 

ioo+jii ./-ioo+ju 2m VZ 

Because of the square root, the integrand has a branch-cut along the negative Z semi-axis. The 
contour of integration for the Laplace inversion can continued to a Hankel contour around the cut, as 
shown in Fig. 3.11. We have 

HM.A.) = - ^0 (- r 1^.-^ - r = 

V J-^oo+ 2m i^-Z Jo- 2m -%^-Z 



= S{A^ - A,) [ 
Jo 



dZ'_^_z>A,_ 2 _ 2 6iA2-Ai) 



(A-32) 



\/7r ^/A^ 



in which we can recognize the subleading behaviour of (3.67). Multiplying by -x/Ai we obtain exactly 
the same result as from the fractional derivative. 
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Figure 3.11: The Bromwich-Hankel contour used for the inverse Laplace transform. Integrating along this contour 
the function e^^^ /\/^ we obtain zero because no singularities are encircled. Prom this, in the limit R — » oo, e — * 0, we 
obtain the result (A-32). 

For the cosmological term wc have again a convergence problem because instead of the square root 
in the denominator we have a power 3/2 which diverges too fast at zero. This can be removed by 
a formal integration by parts in the previous step, and the final result coincides with the fractional 
derivative one. 

The remaining terms give instead no problems. The constant term gives a term proportional 
to S{Ai)S(A2), which is a non-propagating and non-universal kind of term already familiar in two 
dimensions (see [57] and references therein), and which furthermore gives just zero when multiplied 
by the inverse entropy factors aY^A2^^, in agreement with the fractional derivative result. The last 
term has no branch-cut, and is the usual term encountered in two dimensions (see, for example, [51]). 
After inverse Laplace transform it gives 

H'LS{A2 - A,) = (-^2^ - ^) 5iA, - A,) , (A-33) 

Dividing H'|^^^ by the entropy factor l/\/A2 we find the kinetic part of the auxiliary Hamiltonian 
(3.70): 

Finally, by use of (3.71) (where N'iA) has to be replaced by the subleading term alone of (3.67) because 
the exponential part is absorbed in the critical value of the boundary cosmological constants), we find 
the kinetic term of the final Hamiltonian (3.79). 

We have seen that the convergence problems in the two methods are complementary. The only 
problem encountered with the auxiliary Hamiltonian method is in the cosmological term, which instead 
presents no problems with respect to the fractional derivative. 
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